UNCLASSIFIED 


AD  NUMBER 


AD825616 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  MAY  1967.  Other  requests  shall 
be  referred  to  Air  Force  Cambridge 
Research  Labs.,  Hanscom  Field,  MA. 


AUTHORITY 


AFCRL  ltr,  22  Dec  1971 


THIS  PAGE  IS  UNCLASSIFIED 


AFCRL- 67-0280  (Volume  U 
MAY  1967 

SPECIAL  REPORTS,  HO.  64  (1) 


MICROWAVE  PHYSICS  LABORATORY  PROJECT  4642 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 


L.  C.  HANSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


Proceedings  of  the  Third  Symposium  on  the 
Plasma  Sheath  -  Plasma  Electromagnetics 
of  Hypersonic  Flight 

Volume  L  Radiation  Characteristics  of 
Plasma-Covered  Antennas 

W.  ROTV.AN 
H.  MOORE  * 

R.  PAPA 
J.  LENNON 

Editors 


*  Lowell  Technological  Institute  Research  Foundation,  Lowell,  Mass. 


OFFICE  OF  AEROSPACE  RESEARCH 
United  States  Air  Force 


Abstract 


This  volume  is  a  collection  of  the  unclassified  papers  presented  at  the  Third 
Symposium  on  the  Plasma  Sheath  -  Plasma  Electromagnetics  of  Hypersonic  Flight. 
This  symposium  consisted  of  the  review  of  progress  in  reentry  communication 
studies  during  the  three  year  period  since  the  prior  meeting.  The  program  of  this 
symposium  on  plasma  electromagnetics  of  hypersonic  flight  involves  a  wide  range 
ofvscientific  disciplines,  including  electromagnetics,  aerodynamics,  aerothermo- 
chemistry,  plasma  dynamics,  electronics,  and  high -temperature  phenomena.  The 
papers  were  selected  to  explore  as  many  of  these  facets  of  research,  including  the 
results  of  laboratory,  flight,  and  system  tests,  as  time  permitted. 


Preface 


This  volume  is  a  collection  of  the  unclassified  papers  pre¬ 
sented  at  the  Third  Symposium  on  the  Plasma  Sheath  -  The  Plasma 
Electromagnetics  of  Hypersonic  Flight.  The  three  day  meeting, 
which  examined  the  effects  of  ionized  flow  fields  on  the  transmission, 
re  ception,  and  diffraction  of  radio  waves  from  aerospace  vehicles, 
was  held  in  Boston,  Massachusetts  on  September  21,  22  and  23,  iyb5, 
under  the  sponsorship  of  the  Micrcwave  Physics  Laboratory  of  the 
Air  Force  Cambridge  Research  Laboratories,  Air  Force  Office 
of  Aerospace  Research. 

The  term  Plasma  Sheath  has  been  defined  in  the  two  prior 
symposia  and  through  technical  usage  as  an  envelops  of  ionized 
gas  that  surrounds  a  space  vehicle  as  it  reenters  the  atmosphere 
at  hypersonic  velocities.  Related  topics  which  were  considered 
to  be  within  the  scope  of  this  symposium  include  other  types  of 
ionized  flow  fields,  such  as  the  exhausts  of  chemical,  electrical, 
and  nuclear  rocket  prcpulsion  devices  and  the  wakes  of  reentry 
vehicles.  Also,  the  hypersonic  velocities  which  may  be  attained 
during  launch  and  guided  flight  of  high-performance  missiles  and 
spacecraft-give  rise  to  similar  ionization  effects  as  the  reentry 
plasma  sheath.  These,  ionized  flow  fields  can  all  cause  perturbation 
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or  blackout  o £  space-to-ground  radio  links  and,  therefore,  are  of 
concern  in  the  design  of  aerospace  communications,  telemetry, 
reconnaissance,  ard  electronic  countermeasure  systems.  The 
related  topics  of  radar  cross-section  and  detectability  of  reentry 
bodies  and  their  wakes  were  deliberately  excluded  from  major  con¬ 
sideration  since  these  subjects  are  adequately  treated  in  other 
specialized  meetings. 

A  major  goal  of  the  Third  Plasma  Sheath  Symposium  was  the 
review  of  progress  in  reentry,  communication  studies  during  the 
three  year  period  since  the  prior  meeting.  Several  experimental 
reentry  programs  of  the  military  and  space  agencies  were  initiated 
and  brought  to  peak  activityin-this  interval.  Examples  of  these 
rocket  programs  include  the- Air  Force's  Advanced  Ballistic  Reentry 
Systen-.s  (ABRES)  tests,  flight  studies  of  hypersonic  lifting  reentry 
vehicles  under  the  Air  Force's  Projects  ASSET  and  START  and 
NASA's  Project  FIRE,  and  specific  tests  of  reentry  communication 
techniques  under  the  Air  Force's  Trailblazer  II  flight  program  and 
NASA's  Project  RAM.  These  tests  have  demonstrated  many 
important  principles  of  reentry  communication  systems  design, 
including  the  advantages  of  the  higher  microwave  frequencies  and 
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of  chemical  additives,  and  the  development  of  techniques  for  optimum 
antenna  design  and  location.  Several  other  approaches  for  reducing 
reentry  signal  attenuation,  such  as  techniques  which  use  strong  magnetic 
fields  or  electron  beams,  have  been  studied  extensively  in  the  laboratory, 
but  have  not  yet  been  fully  simulated  or  flight-tested.  The  substantiation 
of  these  techniques  has  been  aided  through  the  development  of  shock  and 
wind  tunnel  facilities  which  permit  realistic  simulation  of  reentry  conditions. 
Theoretical  and  laboratory  programs  have  also  advanced  such  topics  as 
the  radiative  properties  of  plasma-covered  antennas,  nonlinear  inter¬ 
actions  of  plasmas  with  high-power  microwave  signals,  and  non-equilibrium 
flow  field  calculations.  However,  the  requirements  of  ballistic  missile  — ■- 

systems  and  manned  space  flight  programs  demand  an  even  further 
increase  in  reentry  communication  and  related  research. 

The  program  of  this  symposium  on  plasma  electromagnetics  of 
hypersonic  flight  involves  a  wide  range  of  scientific  disciplines,  including 
electromagnetics,  aerodynamics,  aero-thermochemistry,  plasma  dynamics, 
electronics,  and  high-temperature  phenomena.  The  papers-were  selected 
to  explore  as  many  of  these  facets,  of  research,  including  the  ruoiilts 
of  laboratory,  flight,  and  system  tests,  as  time  permitted. 
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The  paper*  presented  in  the  various  sessions  are 
reported  in  the  Symposium  Proceedings.  This  consists  of 
two  separate  sets  of  Air  Force  Cambridge  Research  Laboratories 
reports.  The  classified  papers  from  the  final  session  on  Flight 
Tests  and  Alleviation  Techniques  are  in  a  report  which  has 
already  been  distributed.  Requests  for  additional  copies  of 
this  SECRET  document  by  agencies  of  the  Department  of  Defense, 
their  contractors,  aud  other  Government  Agencies  should  be 
directed  to  the  Defense  Documentation  Center,  Cameron  Station, 
Alexandria,  Virginia,  22314.  In  ordering.  Report  Number 
AFCRL-66-608  (Special  Reports  #53,  Aug  66)  should.be  referenced. 
These  present  reports  contain  the  unclassified  papers  from,  the 
other  sessions.  The  papers  are  not  presented  in  the  same  order 
as  at  the  Symposium,  although  the  papers  in  each  section  treat 
related  topics.  In  addition,  to  assist  the  reader  who  is  interested 
in  the  complete  Proceedings,  the  unclassified  titles  and  abstracts 
of  the  papers  in  the  classified  report  are  included  here. 

We  wish  to  express  our  sincere  thanks  to  those  persons 
who  helped  to  organise  the  symposium:  to  the  members  of  the 
Technical  Paper  Committee,  including  Paul  Huber,  Theodore  Sims, 


and  Robert  Rawhauser,  who  undertook  the  difficult  task  of 
paper  selection;  to  the  Session  Chairmen^Dr.  Ronald  Row, 

Or.  Tetsu  Morita,  Or.  Leopold  B.  Felscn,  Or,  Warren  McBee, 
Dr.  Romayne  F.  Whitmor,  Dr.  Glen  Pippert;-to  Mr.  Charles  E. 
Fills,  Symposium  Director;  and  to  Miss  Alice  Cahill  for 
secretarial  assistance. 

Walter  Rofenan 
Howard  Moore 
Robert  Papa 
John  Lennon 
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I.  A  STUDY  OF  THE  MICROWAVE  PROPAGATION  CHARACTERISTICS  OF  A 
MODEL  ENGINE  EXHAUST  USING  A  GELLED  METALLIZED  LIQUID  FUEL 

R.E.  Compton,  Jr. 

Electronics  Department 
Martin  Company,  Martin-Marietta  Corporation 
Denver  Division 
Denver,  Colorado 

This  report  describes  a  series  of  microwave  attenuation 
measurements  that  were  performed  on  an  11.6%  scale  model  of  an 
improved  Titan  II  engine.  The  model  engine  burned  Alumizine  and 
nitrogen  tetroxide  and  was  fired  into  the  atmosphere  at  an 
altitude  of  6100  feet. 

A  microwave  diagnostic  scheme  was  utilized  to  determine 
the  electron  density  and  collision  frequency  of  the  exhaust  plume 
at  the  exit  plane  and  on  axis  of  the  exhaust.  Measurements  were 
recorded  simultaneously  at  16  and  21  Gc  and  the  resulting  atten¬ 
uations  in  the  order  of  25  db  were  used  to  compute  the  plasma 
properites.  Results  of  the  measurements  show  that  the  free 
electron  density  and  collision  frequency  are  an  order  of  magnitude 
higher  than  those  results  measured  on  a  full  scale.  Titan  II 
sustainer  engine  which  uses  similar  liquid  fuels  without  additives; 

The  addition  of  powdered  aluminum  to  a  liquid  propellant 

engine  creates  a  higher  exi t  temperature  and  a  gas-particle  exhaust 

medium  that  is  ionized  to  a  high  degree.  Test  results  verify  that 
% 

high  performance  propellants  increase  the  severity  of  microwave 
propagation  through  the  exhaust  plume.  The  improved  Titan  II  plume 
characteristics  are  compared  to  other  liquid  and  solid  propellant 
rocket  engines. 
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U.  CO)  ANALYSIS  OF  COMMUNICATIONS  ATTENUATION  DATA  FOR 
BLUNT  AND  SLENDER  REENTRY  VEHICLES  IN  FLIGHT 


R.  H.  Edsall 
P.  E.  Bisbing 
F.  H.  Nuaerich 
General  Electric  Company 
Philadelphia*  Pennsylvania 

The  attenuation  of  communications  signals  transmitted 
fro.i  reentry  vehicles  in  flight  is  assumed  to  be  attributable 
to  vhe  effects  of  ionization  formed  in  the  flow  field  surrounding 
the  vehicle.  Theoretical  analyses  of  hypersonic  flow  are  given, 
including  both  viscous  and  inviscid  flow  regions.  The  effects 
of  chemical  nonequilibrium  and  heat  shield  contaminants  are 
included  in  the  analysis  of  flow  field  ionization.  Nonequili¬ 
brium  effects  are  found  to  predominate  at  high  altitudes.  On 
tha^  other  hand,  the  ionization  in  the  flow  field  surrounding  a 
slender  vehicle  at  low  altitudes  is  very  strongly  affected  by 
the  ablating  heat  shield  material.  The  presence  of  alkali 
metals  in  trace  quantities  can  significantly  enchance  the 
ionization  relative  to  its  corresponding  level  in  pure  air. 

The  attenuation  of  communications  signals  by  the 
reentry  induced  plasma  is  treated  theoretically  by  considering 
the  problem  of  the  interaction  of  plane  electromagnetic  waves 
with  a  one-dimensional  non-uniform  plasma,  rhis  approach  is 
used  in  a  general  parametric  study,  since  it  is  not  subject 
to  the  geometry  of  the  transmitting  system.  Parametric  results 
have  been  obtained  by  numerical  integration  of  Maxwell's 
equations  through  the  plasma  profiles  obtained  from  the  flow 
field  studies,  applying  the  required  boundary  conditions. 


The  results  indicate  in  general  that  blunt  vehicles  are  pre¬ 
dominantly  affected  by  the  volume  electron  density  in  the 
shock  layer;  i.e.,  blackout  occurs  whenever  the  shock  layer 
ionization  is  overdense.  On  the  other  hand,  for  the  slender 
vehicles  the  plasma  may  be  treated  as  a  thin  sheet,  in  which 
the  surface  density  of  electrons  determines  the  propagation 
parameters.  The  latter  result  is  a  consequence  of  the  fact 
that  the  flow  field  ionization  is  confined  to  a  fraction  of 
the  thin  viscous  region  (boundary  layer)  near  the  vehicle 
surface. 

Experimental  data  from  ballistic  missile  flights 
involving  both  slender  and  blunt  reentry  vehicles  are  analyzed. 
In  the  blunt  body  category,  signal  attenuation  and  ionization 
probe  data  from  the  Mark  6  are  analyzed.  Signal  attenuation 
data  for  the  RVX-2A  vehicle  are  also  analyzed.  The  slender 
vehicles  included  for  analysis  of  flight  data  are  REX,  WAC, 

RHV,  and  TVX.  In  general,  it  is  found  that  the  important 
qualitative  and  quantitative  characteristics  of  the  theory 
are  verified  by  the  results  of  the  flight  experiments.  In 
particular,  it  is  found  that  the  use  of  a  non-ablating  conical 
body  reduces  the  signal  attenuation  to  a  negligible  amount, 
for  IC8M  reentry  conditions. 


TTT.  ANALYSIS  OF  PLASMA-SHEATH  ELECTRON  DENSITY 
MEASUREMENTS  AT  ENTRY  VELOCITIES 

William  L.  Grantham 

Experimental  results  are  given  of  a  flight  test 
designed  to  measure  the  electron  density  in  the  flow  field  of  a 
hypersonic  vehicle.  A  comparison  is  made  of  the  flight  measure¬ 
ments  with  nonequilibrium  flow  field  calculations  in  order  to 
assess  the  applicability  of  the  flow  field  theory.  The  microwave 
reflectometer  technique  used  to  make  the  measurements  is  presented 
with  attention  given  to  measurement  accuracy  and  resolution  capa¬ 
bility. 

The  experimental  payload  discussed  was  launched  from 
the  NASA  Wallops  Station,  April  10,  1964,  and  achieved  a  maximum 
velocity  of  18,143  ft/sec  during  the  ascending  portion  of  the 
flight.  A  three-frequency  microwave  reflectometer  system  in  the 
payload  measured  the  plasma-sheath  electron  density  and  standoff 
distance  at  discrete  locations  along  the  body.  A  heat-sink 
type  of  nose  cone  was  used  to  keep  the  plasma  free  of  ablation 
products  during  the  data  period. 
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17.  AN  EXPERIMENTAL  FEASIBILITY  STUDY  OF  INJECTANT  MATERIALS 
TO  ALLEVIATE  MARS  ENTRY  COMMUNICATIONS  BLACKOUT* 

:  Arnold  J.  Kelly 

Jet  Propulsion  Laboratory,  California  Institute  of  Technology 

Pasadena,  California 

A  number  of  Mars  lander  mission  profiles  currently 
under  consideration  involve  atmospheric  braking  and  require  that 
continuous  communication  be  maintained  with  the  spacecraft.  For 
these  missions  to  be  successful  therefore,  requires  that  some 
means  be  provided  to  reduce  the  debilitating  effects  of  shock 
produced  ionization  (expected  to  occur  early  in  the  entry 
maneuver)  upon  the  fixed  frequency  telemetry  signal. 

An  experimental  program  has  been  initiated  to  establish  ' 
the  feasibility  of  using  appropriate  materials,  injected  into 
the  spacecraft  wake,  as  a  means  of  reducing  the  free  electron 
density  level  sufficiently  to  eliminate  its  interference  with 
telemetry.  It  has  been  determined  that  the  flow  conditions 
(pressure,  temperature,  electron  density)  anticipated  behind  the 
bow  shock  of  a  blunt  Martian  entry  body  can  be  accurately  simulated 
by  the  plume  from  an  energetic  (Ducati  type)  arc. 

Four  such  arcs,  operating  simultaneously  in  our  vacuum 
facility,  provide  a  large  (~S0  cm  dia) ,  uniform  plume.  Immersed 
in  this  plume,  its  upper  surface  parallel  and  contiguous  with  the 
centerline  of  the  plume,  is  an  aerodynamically  "clean”  flat  plate, 
through  which  various  electron-capture  materials  are  injected  into 
the  str  'm.  The  electron  density  profile  existing  down  stream  of 


this  plate,  with  and  without  the  presence  of  in j actants,  is 
determined  by  the  use  of  a  swinging  aicrowave  diagnostic  arm. 
This  aicrowave  system,  operating  at  a  frequency  of  24  6CPS, 
wjnitors  electron  densities  in  the  range  of  1  to  7  x  10  cm 
with  a  spatial  resolution  of— 4  ca. 

Electron  density  profiles  data  are  obtained  at 
several  downstream  stations  for  various  injectant  materials. 
From  these  initial  tests  it  is  possible  to  make  estimates  of 
the  relative  efficacy  of  the  species  used. 

The  results  are  classified  ana,  therefore,  the  paper 
should  be  presented  at  the  classified  session. 

*  This  work  presents  the  results  of  one  phase  of  research 
carried  out  at  tho  Jet  Propulsion  Laboratory,  California 
Institute  of  Technology,  under  Contract  No.  NAS7-100,  sponsored 
by  the  National  Aeronautics  and  Space  Administration. 


7.  OEPLETIOH  OF  FREE  ELECTRONS  DUE  TO 
WATER  INJECTION  INTO  HYPERSONIC 
FLOW  FIELDS 

I.E.  Beckwith  and  D.H.  Bushnell 
NASA  -  Langley  Research  Center 

The  results  of  the  RAM  B2  flight  test  of  the  effect  of 
water  injection  on  radio  attenuation  (NASA  TM  X989)  are  analyzed 
for  electron  depletion  mechanisms.  The  mean  concentrations  of 
water  droplets  (or  ice  particles)  in  the  flow  field,  as  required 
in  the  analysis,  are  based  on  correlations  of  maximum  spray 
penetration  obtained  from  wind  tunnel  tests  (NASA  TMX-989). 

The  initial  mean  droplet  diameters  are  computed  from  correlations 
for  cross-current  breakup  of  liquid  jets  in  airstreams  (HACA  TN  40S7). 
A  quasi-one-dimensional  theory  for  the  mean  droplet  motion  and 
evaporation  is  described  briefly  and  used  to  compute  flow  conditions 
and  species  concentrations  downstream  of  the  injection  site.  The 
analysis  is  limited  to  injection  from  the  side  ports  of  the 
RAH  B2  vehicle. 

Both  equilibrium  and  finite  rate  homogeneous  reactions  are 
considered.  In  the  equilibrium  case,  the  only  ionic  reactions  of 
any  importance  were  the  formation  of  N0+  and  GH~*  The  free  elec¬ 
tron  concentrations  were  generally  below  the  critical  value  for 

8  3 

VHF  transmission  (about  7  x  10  electrons  per  cm  )  at  a  distance 
of  one  foot  downstream  from  the  injection  site  as  computed  on  the 
basis  of  equilibrium  reactions. 

Since  equilibrium  conditions  are  unlikely  at  the  low 
densities  encountered,  the  assumption  of  finite  rate  ionic  reactions 


should  be  more  realistic.  The  concentrations  of  neutral  species 
were  assumed  the  ease  as  for  the  equilibrium  calculations  and  the 
enly  ionic  reaction  used  was 

K  ♦  0  «  HO*  ♦  e“ 

with  the  finite  reaction  rate  coefficients  taken  from  the  data  of 
Lin  and  Teare.  Electron  depletion  due  to  this  reaction  was 
insufficient  to  account  for  the  observed  signal  recovery  even  with 
the  considerable  cooling  effect  due  to  water  evaporation. 

This  result  suggests  that  heterogeneous  ionic  reactions 
(that  is,  those  reactions  that  take  place  at  or  near  the  surface 
of  the  droplets)  must  have  been  responsbile  for  the  observed 
effects  of  water  injection.  In  the  analysis  of  these  possible 
heterogeneous  reactions,  the  detailed  mechanisms  involved  were  not 
considered.  The  minimum  water  flow  rates,  injected  during  the  flight 
test,  were  used  to  evaluate  a  factor  which  is  considered  to  be  the 
radius  of  an  effective  cross-section  area  for  electron  depletion 
associated  with  the  individual  droplets.  The  application  of  this 
factor,  as  derived  from  the  RAH  B2  results,  to  predicting 
injection  flow  requirements  under  other  conditions  is  considered. 
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71.  <U)  GEHIHI  REENTRY  COMMUNICATIONS  EXPERIMENT 

L.C.  Schroeder 
NASA  -  Langley 


This  presentation  will  discuss  a  flight  experiaent  in 
which  water  was  injected  frca  a  Seaini  spacecraft  into  the 
reentry  plasna  sheath  to  reduce  the  attenuation  of  RF  coa- 
nunications  and  thereby  increase  the  strength  of  the  signals 
received  at  ground  stations.  Three  different  flow  rates  were 
pulsed  into  the  flow  field  between  the  altitudes  of  270  K  ft. 
and  200  K  ft.  and  effects  on  three  frequencies t  VHF  Voice 
(296.8  BC),  UHF  TM  (230.4  sc)  and  C-’aand  (S69D  esc)  were 
observed.  The  origin ,  deveiopsent  and  description  of  the 
experinent  will  be  reviewed.  A  brief  discussion  of  studies  to 
deteraine  penetration  and  flow  requirements  and  to  predict 
antenna  patterns  in  the  reentry  plasma  will  be  presented. 

Flight  signal  strength  results  will  be  analyzed  as  a  function 
of  water  flow  rate  and  aspect  angle.  Comparison  will  be  made 
of  measured  with  predicted  antenna  patterns  and  flow  rates 
required  to  produce  signal  recovers. 
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VU.  EFFECTS  OF  WATER  DROPLETS  ON  REENTRY  PLASMA  SHEATHS 

Shelby  C.  Kurzius  and  Roy  Ellison 
.AeroChem  Research  Laboratories,  Inc. 

Princeton,  New  Jersey 
a  subsidiary  of  Pfaudler  Permutit  Inc. 

Experimental  studies  of  the  effects  of  water  sprays 
on  suppressing  ion  aui  electron  concentrations  in  seeded  low- 
pressure  combustion  plasmas  have  been  made  to  determine  the 
mechanism  by  which  Witer  sprays  are  effective  in  restoring 
communications  with  blacked-out  reentry  vehicles.  These 
stui  is  were  performed  in  hydrogen-hydrocarbon-air  flames  at 
pressures  of  from  10  to  100  Torr.  Langmuir  and  ion  mass 
spectrometric  probes  were  employed  to  measure  total  ion  and 
electron  densities  ar.d  to  identify  individual  positive  and 
negative  ion  species.  Large  decreases  in  plasma  free  electron 
concentrations  were  observed  upon  introducing  water  sprays  to 
the  plasmas  studied.  These  decreases  have  Lean  shown  not  to 
involve  water  vapor.  Attachment  processes  involving  droplets 
are  similarly  not  involved  in  the  mechanism — nor  are  negative 
ions.  The  results  are  compatible  with  a  rapid  recombination 
process  involving  droplets  as  sticky  third  bodies.  Experimental 
observations  of  the  effect  of  varying  the  mean  droplet  size  and 
droplet  concentration  .on  enhanced  recombination  rates  are 
interpreted  in  terms  of  theoretical  predictions  and  the  implica¬ 
tions  of  these  experiments  to  the  alleviation  of  attenuation 
levels  experienced  in  communications  made  through  reentry  plasma 
sheaths  are  discussed. 
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I.  WAVEGUIDE  AIKCTTAXCE  FOR  RADIATZO*  HOT  rUJSHK  1ATESS 
TBBORT  AMD  EXFERM3T 
by 

Janie  Gale J 3  and  Michael  H.  Hsntsoai 

Applied  Research  Laboratory 

Sylvania  Electronic.  Systems 
A  Division  of  Sylvania  Electric  Products,  Ine. 

Maltbaa,  Massachusetts  02154 

ABSTRACT 

Hcvaured  admittance  of  an  X-band  waveguide  radiating  into  dielectric  and 
plasm  layers  is  ccarparcd  with  variatioaally  caaputed  admittance  figures. 

The  waveguide  adaittesce  measurements  for  polystyrene  sheets  of  varying 
thicknesses  cospared  closely  with  calculations.  The  admittance  is  computed 
first  by  iscuaing  the  principal  \*ve  guide  mode  and  second  by  using  a  super¬ 
position  of  alne  and  shifted  cosine  waves  as  trial  functions  for  aperture 
fields.  Although  the  aperture  fields  are  shewn  to  differ  significantly  in 
the  two  cases,  the  computed  admittance  data  are  nearly  the  same  because  of 
the  stationary  character  of  the  admittance  expression. 

The  admittance  measurements  for  plasma  layers  are  made  during  the 
diffusion  controlled  afterglow  of  a  pulsed  discharge  contained  in  a  bell- 
shaped  vessel.  Thin  Tteflon  sheets  are  used  to  confine  the  plasm  layers  to 
an  approximate  size  ox'  5x5x  3 fk  inches.  The  electron  density  of  the  plasma 
profile  is  measured  by  Langmuir  probes  at  various  times  during  the  afterglow. 
The  plasma  density  was  measured  over  approximately  TO  percent  of  the  thiciaess 
of  the  plasm  layers,  and  was  found  to  'as  nearly  constant  in  this  range  which 
excluded  regions  near  the  boundary.  Rome ri cal  solutions  of  the  equations  for 
aobipolar  diffusion  during  the  afterglow  in  a  rectangular  geometry  show  that 
the  lateral  p-asaa  density  variations  are  negligible  in  the  vicinity  of  the 
waveguide.  The  measured  admittance  and  plasm  density  data  are  shown  to 
agree  with  calculations  made  for  homogeneous  plasma  layars,  if  a  correction 
is  made  for  elevated  Ion  temperatures  far  times  of  approximately  10  to  20  ps 
following  the  discharge.  Several  computational  models  consider  plasma 
stratifications  near  the  boundary  of  the  plasma  layer  but  they  do  not  improve 
the  agreement  between  measurement  and  calculations. 
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waveguide  aehitiarce  for  radiatioi  nrro  plasma  layers 

THEORY  AHD  EXPERBfflfT 
b7 

Janie  Galejs  end  Michael  H.  Mention! 

1.  MMOBOCTTO* 

In  past  slot  antenna  measurements  in  the  presence  of  the  oLasna  sheath 
[Cloutier  and  Bachynski,  156 3 1  there  have  bee a  no  attempts  to  correlate 
measured impedances  or  reflection  coefficients  vith  theoretical  calculations 
because  of  the  lack  of  a  suitable  theoretical  models.  The  analytical  techniques 
for  confuting  the  admittance  of  slots  radiating  Into  layered  plasmas  have  been 
developed' recently  for  isotropic  plasma  that  can  be  represented  as  a  simple 
dielectric  of  a-relative  dielectric  constant  of  less  than  unity  [Galejs  1964, 
1965a  and  b;  Villeneuve ,  1965] •  Ibis  paper  will  attempt  to  correlate  the  above 
theoretical  Impedance  calculations  with  laboratory  measurements.  The  simplest 
measurements  are  possible  with  slots  radiating  into  dielectric  layers  and 
after  reviewing  the  theory  in  section  2,  such  measurements  will  be  described 
in  section  3.  Thu  experimental  setup  used  for  generating  a  plasm  layer  will 
be  discussed  in  section  4.  In  section  5  the  measured  admittance  figures  will 
ha  compared  vith  calculations. 

2,  nnn  of  theory 

Variational  expressions  have  been  derived  for  the  slot  admittance  of 
waveguides  radiating  into  plasm  layers.  The  final  expressions  involve  a 
double  summation  [Galejs  1965a],  when  enclosing  the  plasm  layer  by  a  large 
rectangular  waveguide,  or  an  infinite  double  Integral  for  a  plasm  layer  that 
if,  not  bounded  laterally  [Villeneuve,  1965;  Galejs,  1965b].  Both  formulations 
have  been  shown  to  become  equivalent  vith  Increasing  size  of  the  large  waveguide. 
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Use  slot  admittance  is  computed  as  a  function  of  tbe  aperture  fields  of  the 
waveguide.  The  trial  function  for  thfe  fields  and  waveguide  aperture  is  as  stated 
to  be  a  superposition  of  the  sine  and  shifted  cosine  waves,  For  a  rectangular 
waveguide  aperture  which  is  symmetrical  with  respect  to  the  x  and  y  coordinate 
axis  the  trial  functions  are  of  the  form 

fA(x,y)  *  ain[k(|x|4)]  (1) 

fg(x,y)  »  1  -  C08[k(jx|-/)3 

where  k  is  an  arbitrary  parameter  and  Si  is  the  length  of  the  aperture.  The 
amplitudes  of  the  trial  functions  follow  from  the  variational  formulation  and 
it  is  possible  to  compute  ah  approximation  to  the  field  distribution  in  the 
waveguide  aperture  in  terms  of  (l)  and  (s). 

It  is  also  possible  to  ussme  the  principal  waveguide  mode  as  the  aperture 
fields. 

Impedance  calculations  will  be  reported  using  both  the  two-ter*  trial 
functions  (l)  and  (2)  car  the  principal  waveguide  mode  as  the  aperture  fields. 

3.  XMPEDAHCE  FOR  DIELECTRIC  LAXERS 

Impedance  measurements  are  simplest  to  carry  out  far  dielectric  layers. 

At  the  ssme  time  this  provides  a  non-trivial  check  of  tbe  theoretical  develop¬ 
ments  .  The  presence  of  surface  waves  for  dielectric  layers  of  *p/a0  >  1 
causes  additional  singularities  in  the  integrands  of  the  admittance  expression 
Mnd  tbe  numerical  evaluation  of  the  Integrals  using  a  constant  mesh  size 
approximation  can  be  expected  to  be  less  accurate  than  in  the  plasma  case 
(«p/«Q  <  l)  where  there  are  no  surface  waves. 

The  admittance  of  an  X-band  waveguide  radiating  into  plasma  sheath  of 
varying  thickness  d  Is  shown  in  Figure  1.  Tbe  seasurements  axe  mads  with  a 
waveguide  of  aperture  length  and  width  y^.  In  the  computational  models 


[Galejs,  1965*]  the  dielectric  layer  is  assuMd  to  be  enclosed  laterally  by 

a  large  waveguide  of  inside  dimensions  of  xQ  and  yQ.  The  calculations  for 

this  antenna  geometry  are  aade  first  with  the  principal  waveguide  mode  as  the 

electric  aperture  fields  and  also  for  a  two -ter*  trial  function  using  the 

saae  value  of  k.  (k  «  0.625  kQ,  k0  ■  “  ^ ^0«0)  •  two  sets  of  calculations 

agree  closely  except  for  sheath  thicknesses  of  approximtely  0.5  css.  The 

aperture  field  distribution  which  is  computed  with  a  tvo-ter*  trial  function 

is  shown  in  figure  2.  It  exhibits  significant  deviations  from  the  principal 

waveguide  sods.  The  field  distribution  is  nearly  triangular  for  d  «  0,5  ca, 

where  the  two  sets  of  impedance  computations  shown  in  Figure  1  differ  most. 

The  impedance  calculations  have  been  repeated  also  f  or  k  »  k  in  lie  trial 

o 

functions.  These  impedance  computations  do  not  differ  from  those  with 
k  »  0.625  *0,  but  there  are-minor  differences  in  the  computed  field  distribut¬ 
ions  although  the  same  basic  trend  of  the  fields  are  observable.  There  are 
substantial  deviations  of  the  fields  from  those  of  the  principal  waveguide 
mode  in  the  aperture.  However,  due  to  the  stationary  character  of  the  admitt¬ 
ance  expressions  the  principal  waveguide  mode  can  be  used  as  aperture  fields 
for  impedance  calculations . 

k.  LABORATORT  OramTICW  OF  A  PLASMA  LATER 


The  plasma-covered  slot  antenna  was  simulated  in  a  "bell  Jar  vacuum 
system  os  indicated  in  the  diagram  of  Fig.  3.  The  X-b*nd  waveguide  has  a  Mica 
vmcuw  window  and  is  inserted  through  a  vactn*  seal  at  the  base  of  the  bell  Jar. 
This  system  could  be  pumped  down  to  a  pressure  of  about  10 torr  before  filling 
It  with  belli**  gas  of  tank-grade  quality.  Alirndm*  blocks  of  8  *  3/4  x  3/k  in. 


Argon  was  used  initially,  but  it  required  approximately  three  times  lower 
pressure  for  achieving  a  breakdown  and  hence  it  was  more  susceptible  to 
system  Impurities. 
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vere  placed  on  top  of  the  rover  teflon  sheet  approximately  5  Inches  apart  In 
a  direction  parallel  to  the  short  dimension  of  the  rectangular  waveguide 
opening.  A  plasma  layer  symmetrical  vith  respect  to  a  bisecting  plane  parallel 
to  the  ground  plane  is  required  for  obtaining  a  predictable  embipolar  diffusion 
profile.  To  achieve  this,  a  teflon  sheet  vas  placed  on  top  of  the  electrodes. 
This  sheet  also  served  as  a  support  for  a  double  Imngsmir  probe  during  electron- 
density  measurements.  Figure  4  shows  a  photograph  of  the  bell  Jar  and  of  the 
associated  microwave  circuitry.  A  photograph  of  the  operating  discharge  is 
shown  in  Fig.  5.  A  teflon  strip  surrounds  the  electrodes  in  order  to  further 
coniine  the  discharge  to  the  space  between  electrodes.  The  plasma  layer  is 
generated  in  a  bell  Jar  during  sa  active  2  Msec  discharge  and  it  declines 
thereafter  according  to  the  rate  of  electron  removal  froa  the  afterglow.  Both 
the  electron  density  and  the  slot  admittance  were  measured  independently  during 
the  time  following  the  active  discharge. 

4.2  jangndr  Probe  Measurements 

Electron  and  ion  number  densities  and  electron  temperatures  were 
measured  as  a  function  of  time  using  a  double  Langmuir  probe,  Ihe  probe,  shown 
in  Fig.  6,  consists  of  two  1  sa»  tungsten  wires  mounted  5  ns  apart.  Ihe  sires 
are  headed  with  pyrex  glass  and  have  polished  planar  collecting  areas  flush 
vith  the  glass.  To  provide  rf  shielding  the  headed  wires  are  enclosed  by  n 
5/8"  pyrex  tube,  gold  plated  on  the  inside  surface.  All  the  external  circuitry 
of  the  probe  was  housed  in  a  Faraday  cage.  Ihe  probe,  which  was  removed  during 
antenna  admittance  measurements,  vas  placed  directly  opposite  to  the  waveguide 
aperture  at  various  distances  from  it. 

In  reducing  the  double  probe  data  the  Equivalent  resistance  Method  described 
by  Johnson  and  Halter  [1950 J  vas  utilized,  ihe  number  density  of  the  positive 
ions,  n+,  is  given  as 

1/2 


n+  =  const.  (lp/A)(M/T+) 


(3) 
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where  ■  positive  saturation  current,  A  =  effective  collecting  area, 

M  •  ionic  saas,  and  T+  »  ionic  teapcratures.  The  electron  density  can  toe 
coesputed  only  vhen  a  quasi-neutrality  is  aesused  and  the  ion  temperature  and 
Bass  is  knows.  Our  aeasuresents  shoved  that  the  log  ip  vs  voltage  plot  is 
linear.  In  this  case  the  electron  distribution  is  Maxwellian  and  the  electron 
tenperature,  ?g,  can  be  evaluated  (Johnson  and  Halter,  19fO] .  Also  the  electron- 
atoa  awaentus  transfer  collision  frequency,  v  ,  has  been  evaluated  from  pub- 
liehed  cross -section  data.  The  contribution  to  the  total  collision  frequency 
caused  by  electron-ion  collisions  val  is  saall  because  Te  is  high  for  high 
electron-ion  densities  and  v  T  *1,9  [Mentzoni,  1965]. 

SI  6 

Die  double  Iangsuir  probe  neasureaents  yielded  generally  unambiguous 
probe  characteristics  with  a  surprisingly  flat  saturation  curve  and  sharp 
"knees"  as  shovn  In  Fig.  7  fre*  which  the  values  of  and  vq  were  computed. 

Ihe  scattering  of  experimental  points  shovn  in  the  figure  is  typical  for  nost 
of  the  Beasureents.  The  electron  energy  in  the  early  portion  of  the  after¬ 
glow  was  quite  high  ranging  from  I*  to  about  8  eV  depending  upon  discharge 

current  density.  At  8  eV  and  p  =  .9  torr  the  electron  collision  frequency  v 

a, 

in  heliua  is  1.5  x  109  sec”1  yielding  a  typical  ratio  v^/u  «  .025  for  frequenci-js 
in  the  X-band. 

A  typical  electron  density  profile  is  shovn  in  Fig.  8.  The  electron 
density  is  shovn  to  be  approxinately  constant  for  distances  exceeding  1/8  in. 
froa  the  insulating  layer  confining  the  discharge.  Ho  aeasuresents  were  Bade 
for  closer  distances. 

The  variations  of  the  electron  density  profile  in  the  lateral  direction 
have  been  estiaated  by  solving  the  diffusion  equation  for  the  aabipolar  case 
in  discharge  afterglow  for  the  rectangular  voluse  beeveen  the  electrodes  and 
the  teflon  plates  (Mentzoni,  19f>3) •  The  electron  density  is  nearly  constant 
over  opproxiBii  tely  80  percent  of  she  area  and  the  electron  density  gradient 
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appears  to  be  negligible  in  the  direction  parallel  to  the  aperture  plane  in 
the  vicinity  of  the  slot. 

h.3  Admittance  Measurements 

The  waveguide  admittance  vas  measured  in  the  absence  of  plans*  vith 
standard  slotted  line  techniques.  A  balanced  microwave  reflectaaeter  bridge 
vas  used  for  the  plasma  measurements .  A  typical  reflectaaeter  output  signal 
is  shovn  in  Fig.  9*  lac  bridge  is  balanced  to  produce  zero  output  In  absence 
of  plasma,  and  the  output  approaches  again  zero  as  electron  density  c?  the 
discharge  decays.  . > 

The  discharge  vas  ignited  500  tines  per  second.  Using  a  pulsed 
signal  that  is  triggered  free  a  variable  delay  line  at  a  tins  t  following  the 
discharge,  it  is  possible  to  display  the  reflectaaeter  output  at  this  tine  x. 
After  rebalancing  the  re  flee  tone  ter  vith  the  aid  of  a  precision  attenuator  and 
phase  shifter  the  attenuation  and  phase  shift  settings  along  vith  the  reflection 
coefficients  that  have  been  computed  in  the  absence  of  plana*  give  the  required 
inputs  for  an  impedance  calculation. 

5.  AOCTTASCE  FOR  PLASMA  LAYER 

It  follows  free  the  theory  ci'  a  homogeneous  plasma  layer  that  the  conduct¬ 
ance  of  the  slot  is  decreased  and  that  the  susceptance  becomes  more  inductive 
vith  an  increasing  density  of  the  plaeaa  (increasing  plasma  frequency  u  or 

p 

decreasing  dielectric  constant  .  aese  changes  are  mono  tonic  with  increas¬ 
ing  plasma  density.  In  the  experimental  set-up  shewn  in  Fig.  3  it  is  necessary 
to  have  an  insulating  plate  to  Isolate  the  discharge  from  the  ground  plane  and 
also  a  second  dielectric  plate  is  required  to  limit  the  thickness  of  the 
plassa  layer.  The  effects  of  these  dielectric  layers  are  illustrated  for 
two  sets  of  conditions  in  Fig.  10.  An  increasing  electric  thickness  of  the 
top  plate  tends  to  change  the  shape  of  the  admittance  curve.  Even  sore  marked 


dmge s  have.-  >wn  observed  by  increasing  the  electrical  thickness  of  the 
bottca  plate:  a  1/8  inch  thick  pyrex  glass  sheet  (*d/<0  "  ^.52)  causes  an 
increase  of  both  tb';  conduciaace  and  of  the  magnitude  of  the  susceptance  as 
<^/«o  are  decreased  fro*  1.0  to  0.6.  It  is  obviously  desirable  to  keep  the 
tvodl  electric  plats*  as  thin  as  possible  in  order  to  obtain  oyproxirstely 
the  case  adaittance  behavior  as  for  a  plana  layer  without  the  dielectric 
plates.  Sheets  of  15  and  31  ail  teflon  («d/*o  *  2.08)  vere  celected  as  the 
bottaa  and  top  cover  plates  of  the  plaaaa  layer.  Bss  effects  of  plasaa 
stratification  have  been  examined  for  the  four  gnaae  tries  shown  in  Fig.  11. 
Measured  plasaa  profiles  shoved  an  essentially  constant  plasaa  density  for 
distances  exceeding  l/8  in.  frea  the  dielectric  boundaries,  but  there  verc 
no  aeasureaents  aade  at  closer  distances.  The  plasaa  density  will  taper  off 
near  the  boundaries  and  the  geometries  shown  in  Tig.  11  vere  examined  as 
possible  approximations  to  this  taper.  The  caparisons  include  a  hoaogeneous 
plasaa  layer  of  a  plasaa  layer  of  the  saae  total  thickness  but  vith  a 
0.1  in.  boundary  layer  of  («0  +  «p)/2  or  eQ,  and  a  plasaa  layer  vith  0.05  in. 
boundary  layer  of  «o>  The  corresponding  adaittance  data  are  shovn  in  Fig.  12. 
The  decreasing  plasaa  density  near  the  dielectric  boundaries  increases  the 
conductivity  sllgitly,  but  it  affects  susceptance  acre  significantly.  The 
susceptance  increases  particularly  far  a  larger  discontinuity  near  the 
boundary  (Model  3  of  0.1  inch  boundary  layer  cf  cQ). 

The  aeasurtd  adaittance  date  axe  shovn  in  Fig.  13  together  vith  the 
Modal  1  of  the  homogeneous  plasaa  layer,  but  the  tapered  plasaa  codeia  would 
not  improve  the  agreeaent  vith  aeasureaents.  The  solid  dots  are  all  computed 
vith  T+  *  300°K.  The  two  experimental  points  corresponding  to  the  smallest 
values  of  <p  represent  aeasureaents  sade  10  and  20  aicroseconds  after  initiating 
the  discharge.  If  these  aeasureaents  are  corrected  for  ion  tea?«iatures  cf 
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2000  and  600°K  respectively  the  circles  are  obtained  which  yield  good  agree - 
aent  both  for  G  end  3.  Such  elevated  ion  tesperutures  in  tbs  early  discharge 
afterglow  ore  not  unreasonable. 

There  is  e  fair  agreesaent  between  the  theory  and  the  experts  r.ital  data. 
This  theoretical  t<ork  does  not  account  for  the  reflections  caused  by  discharge 
electrodes,  reflections  free  the  bell  Jar  or  for  the  finite  sine  of  the  ground 
plane.  In  view  of  these  uncertainties  it  is  not  expected  to  achieve  a  better 
agreement  between  aeasuresssnts  and  theory  using  the  present  experimental 


procedure. 
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II.  ADMITTANCE  OF  A  ['USHA-CQVERED  CYLINDRICAL  ANTENNA  * 

6.  Meltz,  P.  J.  Freyheit,  and  C.  D.  Lustig 
Sperry  Rand  Research  Center,  Sudbury,  Hass.,  01776,  U.S.A. 

Abstract 

Based  on  a  variational  formulation,  an  expression  is  derived  for  the 
admittance  of  a  circumferential  gap  in  an  infinite  cylinder  covered  by  a 
set  of  coaxial  dielectric  or  plasma  layers.  The  gap  is  fed  by  the  dominant 
TEM  node  in  a  narrowly  -  spaced,  parallel-plate,  radial  waveguide.  Numerical 
computations  are  presented  for  a  multiple-layer  model  of  a  lossy  plasma' and 
for  a  dielectric  layer  which  supports  surface  waves.  These  computations 
are  compared  with  admittance  measurements  of  a  similarly  excited  metal  cylinder 
surrounded  by  an  annular  plasma.  The  plasma  is  produced  in  a  long  double¬ 
wall  cylindrical  vessel  by  a  pulsed  d-o  discharge  in  argon  at  0.4  Torr.  In 
the  experiment,  the  radial  guide  is  excited  by  a  small  axial  probe  fed  by 
a  conventional  coaxial  air  line.  The  admittance  viewed  in  the  coaxial  line 
is  related  to  the  gap  admittance  by  sn  experimentally  determined  two-port 
network  representing  the  junction  region.  The  measurements  are  in  fair 
agreement  with  calculated  values  based  on  a  three-layer  model  of  the  inhomo¬ 
geneous  plasma  and  the  discharge  tube. 


*The  paper  presented  at  the  Third  Plasma  Sheath  Symposium  by 
G.  Meltz,  P.J.  Freyheit,  andC.D.  Lustig  "The  Admittance  of  a 
Plasma -Covered  Cylindrical  Antenna"  has  been  extensively  modified 
by  the  authors.  Due  to  the  modifications  the  authors  have  requested 
that  this  paper  be  presented  only  in  Abstract,  Interested  readers  are 
referred  to  Radio  Science  Vol.  2  (New  Series),  No.  2,  February  1967, 
pp.  203-224  for  tlv  complete  version  of  the  paper. 
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IHTRODUCTIOH 

The  admittance  characteristics  of  three  types  of  olotted -cylinder 
antennas  are  considered,  nassely,  the  long  axial  slot,  the  finite  axial  slot, 
and  the  gap  antenna  backed  by  a  short-circuited  radial  cavity. 

Moat  of  the  present  effort  is  confined  to  a  study  of  the  long  slot  on 
the  cylinder.  This  was  done  because  the  computational  req-oirejoent.fi  are  not 
too  severe,  yet  the  results  should  he  indicative  of  the  trends  experienced 
by  finite  aperture o  in  a  reentry  environment.  Systematic  computations  of 
the  admittance  per  unit  length  are  made  for  various  coating  conditions. 

Firet,  the  case  of  no  coating  is  considered,  and  results  are  compared 
with  those  for  identical  slots  on  flat  ground  planeB.  Sec  card,  the  special 
case  of  the  plasma  resonant  (  a  •  o>p,  y  ■  0  )  coating  is  analyzed  to  deter¬ 
mine  whether  or  not  a  measurement  of  the  input  admittance  at  resonance  will 
yield  diagnostic  information.  Finally,  the  case  of  a  general  coating  is 
treated  to  ascertain  the  'effects  of  loeses  and  plasma  i  nhoaogenelties  on 
the  admittance. 

L-b6l6 


42 


She  exterior  portion  of  finite  apertures  on  coated  rylindr  are  not 
too  difficult  to  analyse  but  can  lead,  to  ccaputatioc^l  difficulties,  even 
in  the  far  field  vhere  a/ysptotic  er.panaicns  nay  oe  sade  to  simplify  the 
field  expressions.  In  the  near  field,  no  such  expansions  are  allowed; 
therefore,  too  computational  requireaents  of  the  external  adsittance  are 
more  severe.  Furthermore,  the  else  of  the  cylinder  chosen  for  the  present 
analysis  (based  on  adaptability  to  feed  vith  the  Blotted  line)  vao  so  large 
that  additional  difficulties  were  introduced  into  the  computations  due  to 
the  large  argusentc  at  the  Bessel  and  Hanhel  functions  involved. 

The  vary  on  the  large  cylinder  includes  specific  confutations  of  the 
conductance  for  the  uncoated  cylinder,  and  experimental  results  for  both 
coated  and  uncoated  cylinders.  Adsittance  expressions  for  the  coated 
cylinder  are  gi  ven  and  preliminary  calculations  of  the  conductance  are 
discussed.  Erperi'eental  and  theoretical  results  for  radiation  pcttems 
are  also  discussed. 

finally,  the  gap  antenna,  backed  by  a  short  circuited  radial  cavity 
and  fed  by  a  current  eleaent  et  the  periphery  1b  briefly  discussed.  Spe- 
ceflc  computations  of  the  impedance  for  the  case  of  no  coating  are  Bade, 
and  expressions  vith  coating  are  given. 

THE  UCPIEITE  SLOT  OB  A  COATED  CXLQ02R 
"he  geometry  at  the  structure  considered  is  show  in  figure  1.  A  long 
slot  is  cut  Into  an  infinite  conducting  cylinder,  and  is  excited  by  an 
electric  field  vhich  is  unifora  along  the  axle,  end  across  the  slot,  and 
vh'ich  varies  in  tine  as  e+Jmt.  Oie  structure  is  coated  vith  a  dielectric, 
Viose  complex  index  of  refraction,  H  nay  vary  in  the  radial  direction. 
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The  tangential  fields,  Eg  and  Fyj,  at  any  point  exterior  to  the 
cylinder  are  described  by  a  Fourier  aeries  in  aziswth,  i.e.. 


Hs(pj0)  * 

gi«<l 

^  ^(p)e-'3^ 


(3a) 


(Xb) 


Rote  that  equations  (l)  ore  independent  of  the  axial  coor'*-Ate,  and 
that  the  aodal  coefficients  H^p)  and  j^g(p)  depend  only  ca  the 
radial  coordinate.  These  coefficients  are  determined  by  a  straightforward 
application  of  the  boundary  conditions  (appendix  H). 

The  complex  power  per  unit  length  radiated  frc*  the  aperture  is  found 
by  integrating  the  radial  component  of  the  Pqynting  vector  across  the  slot, 
i.e.. 


r*W2  rr°!" 

p  *  /  .  ,  Sp«a  /  &j*(a,jJ)Hj,(a,0)d|5 

The  substitution  of  equations  (l)  into  equation  (2)  gives 
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which  is  nothing  more  than  Parseval's  theorem  in  cylindrical  coordinates.1 


3»  external  or  radiation  adaittance  is  defined  by2 
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■  &ce 
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vbere  V0  Is  the  applied  potential  across  the  slot.  Suppose,  no*  that  the 
slot  is  fed  by  a  parallel-plate  waveguide  excited  in  the  TEH  node.  It  is 
asetsed  that  all  higher  modes  are  negligible  so  that  is  the  only  caa- 
pooent  of  tangential  electric  field  across  the  aperture,  and  is  of  the  fora 


*n-T?  (5) 

If  the  slot  is  sufficiently  thin,  then  2^  a  and  d^  »  e  fi0.  There¬ 
fore  the  following  transfora  pair  exists: 
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Te  .!.[#) 
2**  (^o/2) 


(6b) 


Saving  specified  the  electric  field  at  the  aperture,  the  external  adait- 
taace  aay  be  cc^utea- 
SO  CQU3W 

If  the  antenna  radiates  into  free  space,  the  normalized  external 
admittance  per  unit  length  as  tterived  in  appendix  H  is 


(7) 


Equation  (7)  vac  progrwaed  on  on  electronic  eczsputer  and  results  vere 
obtained  far  tvo  cylinders  haring  different  gaps  whooe  circuaference-to- 
vavelength  ratio  ranged  frcs  8.4  to  4J.4.  The  results  are  given  in 


figure  2.  ?or  reference,  the  admittance  of  a  slot  of  the  seme  width  on  a 
flat  ground  plane  ae  given  by  Harrington^  ie  also  plotted  in  figure  2. 

To  the  accuracy  with  which  one  could  read  Harrington' a  graph,  the  ’admit¬ 
tance  of  the  slot  on  the  cylinder  throughout  the  whole  range  of  C  was 
very  dose  to  that  of  a  slot  on  a  flat  ground  plane. 

RESONANT  CQATCHQ 

It  has  been  shown  (ref.  k)  that  only  the  a  *  0  code  ia  supported  at 
the  plasm  rcoanant  condition,  Jf  •>  0,  (i.e.,  v  ■>  0,  top  »)}  therefore, 
the  radiation  patterns  of  an  infinitely  long  slot  (or  the  equatorial  pat¬ 
terns  of  a  finite  slot)  are  circular,  regardless  of  the  sire  of  the  cylin¬ 
der  and  thickness  of  the  sheath  aa  long  as  each  are  finite.  Since  the 
sudden  defoliation  of  the  radiation  pattern  into  a  circle  suggests  a 
method  of  diagnosing  &  reentry  plasma,  the  admittance  was  also  Investigated 
to  determine  its  behavior  at  plassa  resonance  for  possible  scans  of  diag¬ 
nostics.  At  the  resonant  condition,  thie  admittance  per  unit  length 
(appendix  XX)  is: 


The  results  of  equation  (8)  arc  plotted  in  figure  5  a*  a  function  of  C 
end  V.  Again  two  cylinders  were  chosen  with  C  ranging  between  8.%  and 
45.4.  Note  that  for  V  n  1.00  +  e,  which  corresponds  to  a  vanishing 
plassa  thickness,  the  conductance  decreases  frcn  the  free- space  value  by 
an  order  of  magnitude.  Also,  the  eusceptance  can  change  from  a  largo 
capecitative  to  a  sail  inductive  value  or  zero  when  the  structure  is 


costed  vlth  this  resonant  plasma.  As  the  plasaa  beccesea  thicker,  the 
admittance  per  unit  length  approaches  zero,  This  tendency  for  the  admit¬ 
tance  to  drop  at  resonance  suggests  that  a  measurement  of  admittance  could 
also  be  used  for  diagnostics. 

The  admittance  of  the  slot  has  not  yet  been  computed  as  the  index  of 
refraction  Increases  from  aero  to  ease  sail  finite  value.  However,  based 
upon  calculated  changes  in  the  radiation  patterns  due  to  a  departure 

frees  resonance,4  one  concludes  that  additional  azimuthal  nodes  (a  »  1,  2 ,  3, 
etc.)  are  rapidly  introduced  into  the  field  dependence,  at  least  for  C 
ranging  between  22.7  and  U5.4.  As  such,  one  can  expect  that  these  modes 
will  alter,  perhaps  significantly,  the  input  admittance  free  vhat  it  is  at 
resonance;  the  extent  of  the  alteration  increasing  with  increasing  c. 
However,  for  small  C,  this  alteration  may  still  be  sufficiently  small  to 
cause  a  measurement  of  the  input  admittance  near  and  at  plasma  resonance 
to  be  a  useful  plasma  diagnostic  tool.  Me  remains  to  be  investigated. 
aasmjj  cornua 

The  method  described  by  Svift5  was  used  to  analyze  the  slotted  cyl¬ 
inder  coated  with  a  dielectric  having  a  coupler  index  of  reflection,  in 
this  approach,  the  wave  equation  is  numerically  integrated  through  the 
plasma,  thereby  avoiding  the  computational  problem  of  evaluating  Bcss*l  and 
Hankel  functions  of  complex  arguments,  which  describe  the  functional 
behavior  of  the  fields  within  the  coating.  Furthermore,  if  this  metnod  ie 


4? 


Using  Svift's  notation,  the  no  realised  admittance  for  either  hceo- 
geneous  or  inhomogeneous  coatings  (see  appendix  II  for  an  outline)  is: 


Yei 


.  2  «(a)  V-  *»(*)  + 

Js*5-vm*jv« 


(9) 


vhere  the  prise  indicates  differentiation  with  respect  to  the  radial 
parameter. 

Computations  of  equation  (9)  are  plotted  in  figures  5  and  6  as  a 

function  of  ccoplex  index  of  refraction  for  coating  thickness  corresponding 

to  I  £).  m  0.25  end  a  slot  width  corresponding  to  »  0.25  for 

hr  hr 

cylinders  of  site  C  »  h,  8,  and  12.  It  is  important  to  note  thrt  the 
admittance  is  relatively  insensitive  to  variations  in  C.  One  is  therefore 
teapted  to  conclude  that  the  admittance  of  identical  slots  on  cylinders 
and  flat  ground  planes  are,  for  all  practical  purposes,  the  case*  if  C  sJ  4 
and  if  the  loss  angle  &  of  the  coating  is  between  99°  and  lSo°,  and  if 
the  magnitude  of  H  1. 

Flow  field  analysis  shows  that  the  electron  density  and  collision 
frequency  say  very  considerably  within  the  plsssa  sheath.  A  typical  example 
of  the  distribution  along  a  normal  to  the  vehicle  is  shown  in  figure  7,  and 
was  chosen  as  a  coating  for  a  cylinder  of  physical  radius  k.152  ca  sod  an 
aperture  width  of  1.0l£  caT.  admittance  was  computed  ae  a  function  of 
the  exciting  frequency,  and  the  results  are  given  in  figure  8. 

*A  more  general  conjecture  of  this  type  was  suggested  during  conversa¬ 
tions  with  one  of  the  writers  (OST)  by  V.  Retain  of  APCRL  prior  to  the  time 
these  extensive  computations  were  perfrraed. 

^This  width  corresponds  to  the  width  of  a  standard  X-baad  waveg>-ide. 
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nineteen  sdaittance  points  were  ccsputed  in  the  frequency  interval  of 
10.0  to  11.8  tec,  yet  this  maker  vas  insufficient  to  establish  a  saooth 
curve  bu  cause  of  trial  1  -scale  fluctuations.  nonetheless,  ccese  interesting 
features  axe  revealed.  The  aost  striking  effect  occurs  in  the  region  of 
peak  plassa  frequency  (IO.76  tec),  where  the  conductance  decreases  sharply 
and  susceptance  begins  to  decrease  eonotoaically. 

this  seess  to  be  consistent  vlth  the  results  of  the  resonant  plosna 
coating.  At  frequencies  above  resonance,  the  conductance  rises  rapidly, 
but  the  susceptance  regains  relatively  constant.  As  the  frequency 
increases,  the  admittance  should  approach  the  no  coating  values.  Belov 
resonance,  the  curves  are  fluctuating  too  such  to  suggest  any  general 
con elusions. 

Pot  the  types  of  distributions  shows  in  figure  7,  the  peak  plasm 
frequency  seeas  to  be  a  sharp  dividing  point  for  the  admittance  properties 
of  slots  on  cylinders.  Whether  this  is  true  in  general  resins  to  be  sea''. 

TEE  AXIAL  BLOT  OH  A  COATED  CfUMSa 
The  geoeetiy  is  slovn  in  figure  ?.  A  waveguide,  excited  in  the  TEq^ 
atade  opens  cm  to  a  cyiiidrical  ground  plane,  with  the  long  dissension  of  the 
waveguide  parallel  to  'she  axis  of  the  cylinder.  Higher-order  codec  are 
neglected  In  the  vavepuide,  and  the  slot  is  assisted  to  be  thin  enough  so 
that  over  the  slot  •<  Eff  and  a  d0  =  dr). 

Since  the  apertu--e  is  finite,  and  since  the  tangential  fields  vary 
vlth  t;  the  fields  are  described  by  continuous  mode ft  in  axis  and  discrete 
iodea  in  ariauth,  i.e., 

E-(p,0,s)  -  f~  f  ls*(p) c^'e'^dh  (10a) 

Sa-®  ^  -0» 


E0(p,0,s)  -  Y2  f  t^(p)*"JhV,5^dl 1  (10b) 

And  the  cccplex  power  radiated  by  the  aperture  is 

r*c ->o/2 

?  *2  J  J  .  .  s0*v«j0j*)Hs(a,0,s)a  d0  ds  (11) 
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Using  Parseval's  theorem,  equation  (11)  can  he  rewritten  in  the  following 
fora 

€9  pf* 

*  -  ^ac)2*  J  (32; 

And,  the  external  admittance  iar 


The  external  admit  wace  can  now  he  coeluted  using  equations  (13)  and  (lib) 
In  connection  with  the  solution  of  the  boundary- value  problem,  vhich  gives 

SO  CGASSKO 


Jot  this  case  of  no  coating,  the  normalised  external  conductance  and 


susceptance  are 
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Equation  (15a)  van  computed  over  the  X~hand  range  of  6. >5  to  13,12  kmc 
corresponding  to  a  range  .of  slot  lengths  of  l/2  to  1  wavelength  for  the 
TEqi_  mode,  and  the  results  are  shovn  in  figure  10,  For  this  range  of  fre¬ 
quencies  and  the  13-inch-diameter  cylinder  used  (2a  «*  13  in.)  the  parameter 
C  increases  froa  22. 7  at  k  =  0.>  to  k$.k  at  k  =  1.0. 

Aa  a  partial  check  of  the  computations,  the  width  of  the  slot  at 
k  «  0.5  was  allowed  to  approach  zero  in  order  to  compare  the  results  with 
those  of  Wait  (ref.  5)  for  the  thin  resonant  slot.  The  computations  per¬ 
formed  here  givs  -  0,383  compered  to  Wait's  KcvjWalt  “  C.3&.  This 
partial  check  served  as  a  go-ahead  for  proceeding  with  the  other  cases. 


EXPSRIKESTAL  RESULTS  -  SO  COATIHO 


The  particular  size  of  aperture  and  range  of  cylinder  size  C  consid¬ 
ered  in  most  of  the  coaputations  were  choeen  vith  X-band  experiments  ir 
mind.  The  cylinder  shown  in  figure  11  is  1J  inches  la  diameter  and 
24  inches  long,  and  is  just  large  enough  to  contain  a  Hewlett-Packard 
slotted  line.  The  aperture  size  is  0.4  inch  X  0.9  inch,  which  corresponds 
to  the  inner  dimensions  of  a  standard  RQ-52U  waveguide.  The  short  circuit 
was  realized  by  placing  a  small  brass  plate  curved  to  fit  the  cylinder  sur¬ 
face  and  held  in  place  by  a  strap.  Photographs  of  the  end  view  showing  the 
feed  arrangement  and  the  measurement  setup  iB  shewn  in  figures  11  and  12. 

The  experimental  and  theoretical  results  are  shown  in  figure  13.  The 
measured  and  calculated  values  of  the  input  conductance  differ  by,  at  most, 
5  percent.  The  agreement  la  sufficiently  close  to  conclude  that  the  effect 
of  higher-order  modes  and/or  computational  errors  are  negligible  for  the 
large-sized  uncoated  cylinder  used  here.  Agreement  between  first-order 
theory  and  experiment  can  also  be  expected  to  be  good  for  smaller  cylinders 
with  correspondingly  thinner  slots, 

KaPETCWERTAL  JESUITS  -  FOLTETHTLEHE  CQAKHG 

The  cylinder  used  above  for  the  noncoated  condition  was  coated  with  a 
polyethylene  coating  (representing,  electrical! v,  an  ablative  coating)  of  a 
quarter-inch  thickness,  i.e.,  T  ■  b  -  a  =  0.25  inch,  and  measurements  of 
input  waveguide  admittance  were  made.  The  dielectric  constant  of  this 
coating  material  vac  first  measured  over  the  X-band  frequency  range  to 
ascertain  the  correctness  of  the  published  value  of  fij.  =  2.25.  Pour 
samples  were  cut  out  of  the  polyethylene  stock  sheet  and  their  dielectric 
constant  was  measured  by  both  the  Von-Hippel  method  and  by  directly 
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Manuring  tins  guide  wavelength  in  a  slotted  line  completely  filled  with 
the  sample  material.  Both  methods  of  measurement  gave  an  er  which  was 
within  approximately  ±1  percent  of  the  published  value  of  2.85*  This  is 
within  the  accuracy  of  the  measurements;  thus,  one  can  justifiably  take 
the  dielectric  constant  as  2.25  within  this  accuracy. 

The  polyethylene  sleeve  was  "heat  fitted"  on  the  cylinder  so  as  to 
make  a  snug  fit,  and  it  is  estimated  that  the  accuracy  of  concentricity 
of  the  outer  surface  and  inner  surface  of  the  dielectric  coating  was  within 
10.010  inch.  The  slip  fitting  of  the  dielectric  sleeve  on  the  metal  cylin¬ 
der  is  depicted  in  the  photograph  of  figure  14. 

The  measurements  of  admittance  were  made  in  the  conventional  manner  as 
with  the  condition  of  no  coating.  To  realize  the  short  circuit  condition 
and  to  avoid  the  necessity  of  removing  the  dielectric  sleeve  and  then 
replacing  it  at  each  frequency  used,  the  short  was  first  placed  on  the  non- 
coated  cylinder  at  each  of  the  frequencies  to  be  used  and  the  frequency 
setting  was  accurately  determined  by  means  of  a  frequency  meter  accurate  to 
within  ±2j^  me  (Hewlett-Packard  Model  X55HB).  These  exact  frequencies  were 
then  reused  with  the  short  removed  and  the  dielectric  sleeve  in  position. 

Pints  of  gin  and  bj_n  are  shown  in  figures  15  and  1 6,  respectively. 

Measurements  were  then  repeated  for  four  different  circumferential 
positions  of  tbs  dielectric  sleeve  at  each  of  several  frequencies,  as  shown 
in  table  I  and  figures  15  and  16.  The  dielectric  position l.  were  separated 
by  increments  of  90°  and  designated  by  positions  A,  B,  C,  and  D  as  defined 
in  table  I.  It  is  seen  frea  figures  15  and  1 6  that  the  resultant  admit¬ 
tance  differs  for  each  position  of  the  dielectric  even  though  the  frequency 
was  held  constant.  This  can  he  attributed  to  one  or  both  of  the  following 


reasons!  for  each  position  the  "effective  dielectric  constant'1  of  the 
coating  differs  due  to  the  approximately  ±1  percent  deviation  in  the  cir¬ 
cumferential  variation  of  dielectric  constant  which  exists}  for  each  posi¬ 
tion  the  effective  thickness  of  the  dielectric  coating  differs  due  to  the 
inner  and  outer  radial  variation  of  approximately  ±0,010  inch.  In  either 
case,  such  changes  influence  the  mean  electrical  circumferential  length  of 
the  coating,  C,  here  defined  by: 

C  -  ^  H  -  |<W  +  !)K  (17) 


where  a  is  the  Dean  radius  ^a  -  ft  ^  and  H  is  the  mean  refractive 
index.  It  is  seen  that  the  change  in  C  due  to  changes  in  a  (i.e.,  W) 


and/or  H  is 


dC  »  |(W  +  l)m  +  H  £X 


which  since  W  •>  1  and  ■  AT/a  ,  AT  -  ACb-'O,  Is 


/£  £X  +  H  & 

fa  » 


It  is  recognised  that  the  first  tern  contributing  to  dC  is  the 
change  -a  circumferential  electrical  length  due  to  the  change  in  the 
refractive  index  of  the  coating  whereas  the  second  is  due  to  the  change  in 
thickness  of  the  coating.  How,  here  upper  limit 6  of  and  AF/a  are 
approximately  »  ±0.02  and  AT/a  =  +0.0CA  so  that  ££  *•  C  (0.02^), 

i£T^.e  *>  O.Q2hC~.y  For  the  cylinder  used  at  X-band,  CU«y  *>  therefore, 
the  maxima  change  in  C  can  be  <*•  1.0;  i.e.,  the  mean  circumferen¬ 

tial  electrical  length  of  the  coating  can  approach  the  order  of  a  wave¬ 
length  in  the  coating  cater! al.  Intuitively,  oue  would  expect  that  such  a 
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change  'ould  very  veil  l&d  to  ■*.  large  ehcngt  i.  input  admittance.  It  is 
seen  that  the  sajor  contriiiytit to  CO  is  frox  M.  Thus,  for  large 
cylinder!*  (large  C'a)  small  changee  in  the  refractive  Index  my  he 
expected  to  account  for  the  observed  changes  in  admittance. 

To  rale  meaningful  measurements  of  the  i  .put  admittance  and  to  com¬ 
pare  them  vith  computed  result*,  it  follows  tb*;  stricter  tolerances  on  the 
refractive  index  and  coating  thickness  will  he  nr cost ary  if  one  uses  elec¬ 
trically  large  cylinders.  These  tolerances  vill  he  c  -•  aspondingly 
reduced  (as  saen  from  eq.  (19))  for  a an  tier  cylinders.  It  would  seen 
then  that  for  an  initial  test  of  theory  it  say  be  sore  s^proyrie.te  to  use 
smaller  cylinders. 

The  radiation  patterns  also  seemed  to  he  sensitive  to  sr;JLl  changes 
in  the  electrical  parameters.  The  measured  equatorial  plane  pot  inns  are 
shown  in  figures  17(a)  to  17(f)  along  with  the  corresponding  ccr/.  ted 
patterns  fusing  expression  535  of  Walt2).  It  is  seen  that  the  agree  --amt 
between  theory  and  experiment  in  the  forward  direction  is  quite  satis¬ 
factory;  however,  at  the  higher  frequencies  (k  -  0.80,  0.85,  and  O.JG)  the 
radiation  level  measured  in  the  rear  was  considerably  higher  than  that  pre¬ 
dicted,  although  for  the  lower  frequencies  (k  ■  0.65,  G.70,  and  0.75) 
agreement  at  the  rear  is  at ill  satisfactory.  The  poor  agreement  in  the 
rear  direction  cay  also  he  attributed  to  the  critical  dependence  on  tbs 
parameters  of  H  and  W. 

AEiTTKKCE  0?  A  CGftTSD  AXIAL  SLOT 

From  the  work  of  Vait2  the  pertinent  fields  established  in  region  1 
(a  S  P  S  h)  sad  region  2  (b  $  p  S  ®)  of  figure  II-I  are  given  by  (336), 
(337),  (359),  (34o),  and  (5^2),  (3^3),  (3^5),  (346)  of  Wait,  respectively. 


or  the  sake  of  brevity,  these  expressions  viU  not  be  revr_tc>a  here  bu< 
it  is  noted  that  the  foUoving  difference  in  notation  is  used: 


Vait's  notation 

Rotation  of 
this  report 

h 

% 

Furthermore,  lB  noted  that  in  table  I  of  Volt2  (p.  128)  that  the 

coefficient  should  be  multiplied  by  u,  and  the  coefficient  is 

lacking  *  minus  sign. 

Here  again  the  tangential  electric  fields  on  the  cylindrical  sus-face 
p  «  a  ere  assumed  to  be 


20(a,0,s) 


off  slot 


|Eo  cos^j 


on  clot 


Ss(a,0,z)  .  0 


(20a) 

(20b) 


Froa  equation  (20a)  it  follows  that  the  transform  of  ^  is 


and  fren  (339)  of  Wait  the  transform  of  H,  is 


^(hja)  -  ^{b^fyua)  +  BaJgfua)] 


(22) 


''rca  Perceval's  theorem  the  external  (radiation)  adsitiance  in 

gp  a=4=  r*s  _  _ 

*c  “  p-jg  =  J^2*)2  ^  J  ^(hjtjE^afhjaJdh  (23) 

Using  the  nix  tangential  boundary  conditions  ^continuity  of  E0,  Ej, 
80,  Hj.  at  p  *  b,  and  continuity  of  E0  and  Ej  at  p  *»  a^  and  dcterai- 
nants  gives  expressions  for  Bq,  b-j,  and  Dq,  respectively,  vhere  is 
the  deterainant  forced  by  the  coefficients  a^,  etc.,  in  table  I,  p.  128- 
of  Wait. 

Explicitly  solving  for  Ba,  b^,  and  and  substituting  the  expres¬ 
sions  into  equation  (23)  then  gives 

-Jtf  ^ *  {->0  f?JV  } 


vhere: 


Ub  -  JnCua)^2)’^)  -  JB'(ub)^(2)(ua)  (25a) 

VB  -  Ja(ua)%v2){ub)  -  Ja(ub)%(2)(ua)  (25b) 

I*a  -  J,a(ub)BQ(2),(ua)  -  ^'(ua)^2)^)  (25c) 

Ta  -  Ja'fub)^2)'^)  -  JB’ (ua)Sjj^2) ' (ub)  (25d) 

is  e  partial  Chech  on  equation  (24),  consideration  of  the  special  case 
of  no  coating  (b  «  a)  or  an  air  coating  (h  =  1)  each  cause  equation  (24)  to 
reduce  to  equation  (15a)  and  (15b)  as  should  be 
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rationalisation  and  norsalixation  of  equation  (5”  then  gives 

('»*  -  ■  ia)as*2(">5^|(*J  *  r)(*  *  ^  5a^{J) 

v  ^  m  ■  £^z -  iw*c»r) _ . 

fc|Sil  *  *»  t1*^  p  -  -7*)  (i-  tO^/TT?  W’>H> 

-  i*2  -  7*  TA(i,'<V]  [[”  **  *A(i50p 

-  F7?  W2)'(x>]  -  kr>*(**  -  i>^A[*a(I)!I,>]i? 


where  a  CW -  y2,  X -  C  V  ^£r  _y£ 

-  [Ja(x)V(D  -  J=‘(g)^=(x)j 
Ta  -  jjn'(X)V(§)  - 

°n  -  ^'(X)  -  JE'(X}^)J 

va  -  JxD($)*n(x)  -  JB(X)^($)J 

X»3  **  -Jlfe 

Ua  -  -J?fe 

Va  »  -JVE 


It  iu  important  to  note  that  the  above  expression  vas  also  derived  inde¬ 
pendently  by  integrating  tie  Poynting  vector  over  a  spherical  surface  of 
radius  r  in  the  far  field.  As  a  partial  checi  on  equations  (27)  the  case 


of  ao  coeting  (b  «  i.r.,  V  =  1)  or  -^n  air  '  .  titiF  [M  »■  ).)  both  reduce 

equations  (27)  to  equation  (l**)  «  sbvJLd  be. 

It  is  noted  that  tbs  abort  reduction  of  y;  shows  thot  sal,,  the  inte¬ 
gration  fr os  0  £  y  ,  *  contributes  to  the  conductance  (which  nay  be  antic¬ 
ipated  it  one  considers  that  the  conductance  can  also  be-  obtained  using 
Pcyntihg'fi  vector  and  the  radiation  fields). 

However,  the  expression  for  b-  will  be  of  the  fora 


Mb&* 


pi  pH  r 

j  *  J  *  J 

Jq  Jy  Jj 


for  the  general  case  of  an  arbitrary  homogeneous  coating.  As  yet  the 
explicit  fora  for  the  integrands  of  equation  (26)  has  not  been  obtained. 

An  attempt  was  made  to  program  a  modified  form  of  equation 
(27)  in  order  to  .ospute  the  conductance  of  the  axial  slot  en 
the  13-inch-diaaeter  coated  cylinder.  However,  the  results 
seen  *d  to  be  in  error  by  at  least  6  percent .  Possible  errors 
in  this  computation  are  still  being  investigated  by  examining 
a  plot  of  the  integrand  of  equations  (27) . 
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TS':  C3JSDrlZAS-  QAP  AM3HA. 

The  a.'Stew.i  under  isc-oSsiOR,  shown  is  figure  15,  iy,  and  20,  is 
es.  jntially  j  cylindrical  gecser ary  version  o,  a  spi  erica!  slot  antenna 
d’s*'i*ibc.i  by  Musiahe  sod  Webster”1.  S-.-ch  uit  aatena.  b*u  interesting  fea- 
turec  chat  nake  it  extreme ly  useful  as  a  diagnocti"  tool  for  reentry  ple-tr- 
sheaths.  To  none  a  fev: 

1.  Its  cylindrical  structure  with  so  protruding  pert  allows  it  to  be 
an  integral  part  of  the  reentry  vehicle. 

2.  The  feed  systes  consisting  of  the  inner  conductor  of  a  coaxial 
cable  across  a  gap  such  smaller  than  a  wavelength  supports  a  uni for*  cur¬ 
rent;  as  a  result,  solutions  are  possible  without  assuring  an  aut-cure 
field  distribution. 

5-  In  spite  of  such  a  narrow  gap,  the  antenna  is  an  excellent  radiator 
when  its  circumference  is  approximately  equal  to  2  wavelengths.  At  YE?  fre¬ 
quencies,  the  physical  site  of  the  antenna  is  compatible  with  the  site  of 
nary  stall  reentry  vehicles. 

h.  An  additional  feature  of  the  narrow  gap  is  to  rule  out  the  existence 
of  axial  magnetic  fields  at  the  aperture  and  inside  tbs  gap.  Consequently, 
only  ari-juthal  "Ctj?  modes  are  excited  and  this  field  distribution  is  sain- 
tain.-d  *ven  In  the  presence  of  radially  nenhesogeneeus  plaasas. 

5.  The  antenna  is  a  resonant  structure  whose  radiation  and  ispedcuce 
p-qperties  :re  detersined  only  by  one  significant  arisstbal  mode.  Tress  the 
number  of  lobe 6  in  the  radiation  pattern,  it  is  possible  to  determine  the 
si^iificaat  nodal  ispedaoee. 
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6.  Flr.al.ly,  this  is  an  antenna  -which  is  amenable  to  simple  mathe¬ 
matical  analysis,  permitting  accurate  prediction  of  both  far  and  near  field 
performance . 

In  this  paper  the  antenna  input  impedance  both  with  and  without  plasma 
coating  is  discussed.  Experimental  checks  are  only  given  for  no  coating, 
and  are  correlated  to  the  radiation  patterns. 
fKPOT  IMPEDANCE  WITH  NO  PLASMA.  C  GATING 

Consider  the  cylindrical  gap  antenna  depicted  in  figure  18.  The  gap 
width  d  is  narrow  compared  to  the  Tree-space  wavelength  7^  (d  «  and 
is  formed  by  a  z-sdial  waveguide  short-circuited  at  its  center  by  a  spacer 
of  radius  p  =  pQ.  The  gap  is  fed  by  a  coaxial  line  located  at  the  rim  of 
the  guide  at  p  «  a.  The  coax  center  conductor  makes  electrical  contact 
with  the  guide  top  place  and  its  outer  conductor  Ls  shorted  out  against  the 
bottom  plate.  Because  the  gap  width  is  snail,  with  respect  to  the  -wave¬ 
length,  the  current  along  the  coax  inn-.--  co;»luctor  across  the  gap  can  he 
taken  to  he  constant  and  equal  to  I.  If  the  voltage  in  the  coax  is  Vc, 
the  antenna  input  impedance  presented  by  the  gap  to  the  coax  is 

Zin  =  (29) 

This  is  not  merely  a  definition  hut  also  a  measurable  quantity  related  to 
the  reflection  coefficient  T  in  the  coax  via  its  characteristic  imped¬ 
ance  v 

Zin  =  Zc  (30) 

In  order  to  calculate  Z^n,  refer  to  figure  21  which  shows  the  feed  region 
grossly  enlarged.  The  voltages  Vc  and  Va  in  the  coax  and  in  the  aper¬ 
ture,  respectively,  -e  related  by  Maxwell's  equations  to  the  magnetic  flux 


SI 


enclosed,  by  the  line  Integral  of  the  electric  field  yielding  these  voltages, 
that  is, 

v»  -  ‘  -&  If* '  “  <« 

for  a  perfectly  conducting  wire.  The  integration  ia  over  the  surface 
defined  by  the  loop  boundary  just  described.  The  right-hand  side  of  equa¬ 
tion  (31)  is  the  inductive  reactance  of  the  loop,  1 times  the  total 
current  I  flowing  in  the  wires  and  spreading  tnrough  the  radial  guide 
plates.  Equation  (29)  can  thus  be  rewritten  after  dividing  by  T 


_y 

-jS-  is  immediately  recognized  as  the  input  impedance  Zjn  being  sought. 

“V 

The  ratio  ~~  has  the  dimensions  of  an  impedance  and  is  defined  as  the 
antenna  aperture  impedance.  It  will  be  shown  that  the  aperture  impedance 
is  not  in  gereral  directly  measurable,  but  in  some  cases  it  can  be  calcu¬ 
lated  from  a  knowledge  of  the  fields  in  the  aperture.  The  quantity 
can  be  considered  to  be  solely  due  to  the  reactance  of  the  loop. 

Equation  (52)  nay  now  be  expressed  as  fellows'4 

^in  =  (^a  *  zaf)  +  (^fa  +  Zf)  (55) 

-V 

where  the  first  term  in  parenthesis  is  equal  to  and  tne  second  to 

izDgLj .  Each  parenthesis  consists  of  two  parts,  the  self -impedances  Zg, 

Zf,  and  the  mutual  Zej>  w  Zf&.  (The  f  and  a  subscript  denote,  respec¬ 
tively,  feed  and  aperture  impeiax.r  is. )  The  two  impedances,  aperture  and 

*Dr.  George  I.  Cohn  made  the  analysis  described  by  equations  (33) 


feed,  can  sometimes  be  calculated  and/or  measured  independently  of  each 
other  if  the  interaction  caused  by  the  natural  impedances  is  negligible. 
This  is  the  case  if  the  fields  can  be  divided  on  a  spatial  basis  into  sub- 
volumes  such  that  the  fields  in  any  one  subvolua-j  can  be  attributed  only  to 
sources  not  producing  fields  in  any  other  sub volume.  Thus,  if 

M  «  N  (3^) 

or 

|Z|tf|  «  |Zfl  (?**) 

-V 

then  Za  3  is  directly  calculable  fro®  the  fields  in  the  gap  -which  are 
obtained  by  solving  the  electromagnetic  boundary  value  problem-  Zf  is  the 
wire  inductive  reactance  in  the  absence  of  the  gap. 

The  flux  c  ceznon  to  two  sub  volumes  is  always  lesB  than  the  Eelf -linking 
flux;  therefore, 

M  2  j"a|  (35a) 

or 

|Zafj  3  J^fl 

One  way  to  insure  separability  of  impedances  as  given  by  equation  (33)  is 
to  demand 

|Za|  «  iZf|  (36a) 

or 

|*fJ«|Zft|  (#b) 

that  is,  either  the  wire  inductive  reactance  is  much  smaller  than  the  aper¬ 
ture  impedance  or  vice  versa.  It  will  be  shewn  that  the  feed  wire  induc¬ 
tive  reactance  can  be  made  to  fulfill  the  inequ/:>liiy  equation  (36b)  for  this 
gap  antenna  in  question  and  vanishes  as  the  gap  gpen  to  zero.  The  feed 
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impedance  Zf  and  the  aperture  inpedance  (in  tlje  absence  of  coupling 

between  them)  will  now  be  calculated,  keeping  in  mind  that  as  long  as  equa¬ 
tion  (36b)  is  satisfied,  the  antenna  input  impedance  Z^n  is  the  sum  of 
these  two  contributions. 

An  upper  bound  on  the  feed  wire  inductive  reactance  is  established  by 
calculating  the  Inductance  of  a  coaxial  line  of  length  d  with  inner  and 
outer  radii  11  and  rg.  Using  the  geometrical  approximation,  it  is 


easily  shown  that 


If  11-  A.  and  —2.  ~  10. 
30  ^ 


Zf  =  377A.  in  E2.  0 


|Zf|  ~  301  (38) 

Since  the  antenna  is  to  have  an  input  impedance  of  50G  resistive  in  order 
to  match  the  coax  characteristic  impedance,  equation  (5£>)  is  not  fulfilled 
and  the  input  impedance  is  not  separable  into  feed  and  aperture.  Changing 
the  outer  to  inner  coax  radii  ratio  in  equation  (37)  is  not  as  effective  as 
reducing  the  gap  width  since  the  logarithm  varies  3 lowly  with  the  ratio  in 
question.  However,  if  the  gap  width  is  reduced  by  l/2  or  l/3  corresponding 
to  d/?^  of  the  order  of  l/60  to  l/lOO,  equation  (36)  is  fulfilled  and  the 
inpedance  is  separable. 

The  solution  of  the  electrciaagnetic  fields  subject  to  the  boundary 
conditions  yields  the  aperture  admittance,  ^  »  -A-  for  each  peripheral 
mode  m,  whore  Vn  =  Eg,  d,  is  the  voltage  across  the  gap  for  the  mth  mode 
as  shown  in  figure  22.  The  total  aperture  impedance  is  in  turn  related  to 
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n 

1 

ii 


^  4'ii  B  0^  +  i(Bn1  +  B^) 


(59) 


where  ie  the  suoceptance  presented  at  the  feed  by  the  interior 

region,  and  Gm  and  Bne  is  the  admittance  contributed  by  the  exterior 
region.  For  a  narrow  gap,  with  uniform  feed  current,  the  resultant  modal 
admittance  normalized  with  respect  to  the  wave  impedance  in  the 
exterior  space  Y0  {( 120fi )  for  free  space]  has  been  found  to  be 

r*  • ♦  •(».*  ♦  *.') 

f  rl  tit  t3  i  j  rX  tin  e2  ,  1 

\Jo  \i.u){^irr?}'  *  S~?  A  ^pr-x)  ? 

. (*?[ r1  £(*E^>...yr?l*  ,u  ^ g» _ A-p .  r  't'xyP'-jl  «* 

“  liWM'.i  -  ^.wAk)] 


^  W 


{ho) 


x,  a0,  djj,  r0  being  respectively  normalized:  wave  number  k/ko, 
cylinder  radius  k^,  gap  width  k^d,  spacer  radius  koP0.  vas 
derived  by  calculating  the  fields  in  both  regions,  matching  them  at  the 
boundaries  between  antenna  and  exterior  region  taking  into  account  the 
discontinuity  in  current. 

For  each  mode  m,  there  exists  a  current  sheath  JOT  around  the  gap 
periphery  which  is  independent  of  z  (fig.  22).  The  electric  field 
across  the  gap  resulting  frem  this  current  density  is  also  independent  of 
z  at  the  interface.  It  follows  that  the  locally  induced  magnetic  field 
at  the  interface  does  not  vary  along  the  gap  width.  Pursuing  thin  reasoning 
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reasoning  further,  only  e  z-dependent  could  produce  an  E0  component 
and  this  one  in  turr  would  produce  an  Hz.  One  therefore  concludes  that, 
in  the  absence  of  z-dependent  field  components,  no  1L,  component  is  gener- 
-ted  at  the  gap  interface. 

EVALUATION  0?  APERTURE  BCEnANCE  Zq 

Equation  (40)  for  b^,  and  b^  was  used  to  calculate  the 

various  modal  impedances  for  a  frequency  scaled  up  model  of  the  cylindrical 
gap  antenna  under  study.  The  model  shown  in  figures  19  and  20  was  designed 
to  operate  around  1090  me  with  the  following  parameters: 

Gap  width,  d  =  inch,  -0.156 
Spacer  radius,  p0  =*  1  inch 
Gap  radius,  a  »  3  inches,  Bq  =  2.13 
Operating  wavelength,  Aq  =  8.8l  inches 

Graphs  showing  pertinent  calculations  of  gQ,  bQe,  and  'o^m  are  shown  in 
figures  23,  24,  and  25-  The  calculated  normalized  conductance  and  suscep- 
tance  for  the  above  parameters  as  well  as  the  unnormalized  corresponding 
impedances  are  shown  in  table  n  for  mode  numbers  ranging  from  m  -  0  to 
*-  -  9.  The  important  feature  of  these  results  is  the  fact  that  for  an 
antenna  circumference  approximately  two  wavelengths,  the  first  order  mode 
(m  =  1)  contributes  most  significantly  to  the  impedance.  The  modes  below 
(m  -  0)  have  a  reactance  that  tends  to  cancel  the  reactance  of  the  modes 
above  it  (m  5  2).  In  essence,  for  each  value  of  a0  the  antenna  is  a 
resonant  structure,  with  one  contributing  mode.  (For  a^  ~  1  it  can  be 
shown  that  Zq  is  the  resonant  impedance. )  The  resonant  mode  determines 
both  the  impedance  as  well  as  the  radiation  characteristics.  For  the 


specific  choice  of  parameterr.  above  a0  ■  2.13,  the  measured  radiation 
pattern  shovn  in  figure  2 6  exhibits  the  tvo  lobes  that  one  vould  expect 
fr<sn  the  a  *  1  cosinusoidal  aperture  distribution  of  figure  22.  There¬ 
fore,  the  size  of  the  antenna  establishes  only  one  specific  node  as  the 
contributing  one  and  this  node  in  turr  determines  uniquely  both  the 
input  impedance  and  the  radiation  fields. 

APERTURE  IMPEDANCE  IN  TEE  PRESENCE  CP  A  KOHHCMOGEHEOUS  PLASMA. 

The  plasm  aodel  under  consideration  is  represented  by  a  radially 
varying  dielectric  uniform,  in  the  0-  and  z -direction. 

The  approach  followed  is  the  one  developed  by  C.  T.  Swift®.  The 
plasma  is  subdivided  in  n -concentric  cylindrical  ..jeatha.  Each  sheath  is 
taken  to  have  a  uniform  plasma  and  collision  frequency  equal  to  its  aver¬ 
age  through  the  sheath. 

The  fields  in  the  radial  waveguide  are  matched  at  the  interface  to 
those  in  the  first  sheath.  The  field  transform:  In  each  sheath  are  no 
longer  forward  traveling  waves  expressed  in  terns  of  Eankel  functions. 
They  are  standing  waves  consisting  of  the  linear  superposition  of  Hankel 
and  Bessel  functions.  The  fields  in  Region  I  and  the  first-  sheath  of 
Region  H  are 


■  2a 


d  FaCV) 


"-‘--“ogrT 


(M) 


Region  X 
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Region  XI  (let  layer) 


-  ni 

Ezm  - 

%»B1  -  ’W//  -  T1*) 

Hen111  -  -*ol*ol 

S^111  -  U!^1 

+  iZol^olui^m1  ^ 

V111  -  +  Sol*ol 

The  foregoing  expressions  for  Region  I  were  obtained  after  caking  the 
substitution 

uo2  =*  *o2 

uo2amoEno(uoP)  =  *g~ 

In  Region  II,  forward  traveling  waves  are  replaced  by  the  standing 
waves,  l.e., 

cm(^)^i^2Hup)-^#ja1(h1up)  »  aa:L()c)Jia(u1p)  +  Aa1(k)Hm^2^(u1p)| 

da(k)Hn(9)(up)->/IIIi(k1up)  «  V^k^fuip)  +  Ea1(k)^,(2)(u1p)| 

J 

and  ySn  are  derivatives  with  respect  to  the  argument  ^  Uj_p  j . 


(42) 


(45) 


(44) 


k 


The  nodal  aperture  admittance  Ym  is  now  expressed  in  terms  of  these 


unknown  coefficients  by  demanding  that 

H^D1  0<z<d  (U5) 

Using  the  expressions  (Ul)  and  (U2),  it  follows  that 

pt*  I  n  op 

-iY^i  /  +An1Ha(2)  (u^Uje^dk  -  f  /  b^J^a) 

■<#  —CO 

+  BnV2)  (uia)^  dk  +  iY0  ^  ^  0  <  z  <  d  (*6) 

In  order  to  eliminate  the  dependence  on  z,  hoth  sides  of  equation  (h6)  are 
integrated  frets  0  to  d.  And,  use  of  the  relation 

1  f  e^dk  -  C0(k) 

U  v/q 

results  in 

-1*0*01  J ^  [•■VM  *  W2’1  (vi)»l'of<)“  -s  f__  [».ljn(“l») 

*  *  «.  ^  'f  ■  £  «s 

In  order  to  find  the  four  coefficients  %■*■,  Aq1,  b^1,  and  Bjj1 
icust  he  expressed  in  terms  of  VH.  This  is  done  as  follows. 


As  shown  by  8vift#  a  Transfer  Function  can  be  developed  vhich  relates 
via  a  4  x  4  natrix  the  unknown  coefficients,  a^1,  K?-,  b^,  and  B^ 

in  the  first  sheath  to  the  free-space  coefficients  c^dg,  beyond  the  last 
sheath.  The  transfer  function  depends  only  on  the  values  of  Jn(uiPi)  ant 
H^uiPi)  at  each  shea  a  and  the  corresponding  plasaa  parameters.  Hie 


matrljf  coefficients  axe  labeled  C 


Their  values  are  given  in  the  cited 


reference.  The  following  relation  holds  aaong  these  coefficients 


“a1 

C11  c12  C1J  cl4 

0 

A,1 

C21  c22  C2J  c24 

Cq 

^a1 

— 

C31  C32  c55 

0 

Ba1. 

c4l  c42  c4j  c44 

&3i 

m 


The  preceding  equations  reduce  the  number  of  unknowns  free  4  to  2.  Appli¬ 
cation  of  tive  two  tangential  boundary  conditions  at  the  antenna  surface, 
i.e.. 


4a  -  £  Ok) 


W, 


*The  field,  expressions  are  different  frosa  Swift  because  the  factor 
i8  absorbed  in  the  coefficients  aa1(k),  A^-^k),  bg^k),  and 
.  Furthermore,  the  fields  are  constructed  froa  while 

Swift  uses  negative  exponents.  The  above  expressions  can  be  converted  to 
Swift's  by  dividing  all  field  transforms  by  <ru0e^  and  changing  a  and 


k  into  -a  and  -k. 
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allows  the  unknown  coefficient  Cg  and  dg  and,  consequently,  Bq1,  Ag1, 
bg^,  and  Bj.1  to  he  expressed  in  terra  of  Ya.  The  derivation  now 
follows. 

Equation  (48)  can  be  written  as 


«»•  -  c12«n  +  CudJ 
Aa1  “  C22ca  +  C^dg 


t>a  -  CjgCa  +  c34dss 

Bjs1  “  C4'/°b  +  ClMa 

Therefore,  froa  equations  (47),  (/» 9),  and  (43),  after  using  equation  (50), 
it  follows  that 


(50) 


♦  SA(!I]  «.  ♦  (=l^r  ♦  «.  *  ^7  co*(‘) 

> 

♦  <*A(*5  *  tr«i  ^  ♦  <^§ S  »  *  <W*,‘  ‘  IToi  ;“j  jMfu-'.  ♦  «.  * 0 


(51) 

Solving  for  c-,  and  d-,  in  teres  of  Vn  and  expressing  Sg1,  Ag1,  bg1, 
and  Bg1  in  terns  of  Y_  via  equation  (49)  and  substituting  into  equa¬ 
tion  (4l)  results  in  the  following  expression  for  the  nodal  antenna  aper¬ 
ture  admittance  covered  with  plassa. 


*  *al*  J *  (e22cu  -  f  ♦  (*1^44  . 

-  ^7  /__  («:.<>  -  •  (^A.  •  «!•««») (i.*1’)*  *  fcrfu  -  Cj»S4j  .  ctfy. .  c^c^yy-.'21  c-£=  ■l|“  - 


(52) 
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where  A  the  determinant  of  the  !•  X  4  matrix  in  equation  (51)  is  given  by: 
A  *  (^12^ 54  -  -  CihCj^^J-IL^2) 

+  (<^54  -  c24c52)J='2=(2)  +  (c2^44  -  '  (55) 

It  is  understood  that  the  arguments  of  the  cylindrical  functions  ere 


uja  -  y(i£i  -  kZja.  The  quantity  k^  is  the  wave  number  in  the  first 
plasma  sheath  adjacent  to  the  antenna. 

As  a  partial  check  on  the  validity  of  expression  (52),  the  plasma  is 
removed,  in  vhich  case 

Cj2  **  *  Cgjj  b  Cxg  b  Cj4  fl  C]^2  *  o 


Cgg  *  C44  b  1 


*ol  “  *o 


The  admittance  n crs  reduces  to 


.La  r ^g)t c *r  dk  «*  r *(8)  c*c  ^ 
Y°^  JiT  0  0  *  ^  J^VJZT  0  °  T 


-  V  dk  a» 


vhere  C0*C0  *  sin  c2  i|L 

Bruation  (54)  is  equivalent  to  equation  (4o)  prior  to  breaking  it  into 
real  and  imaginary  parts. 
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CGHCUJDING  REMARKS 

Tbits  investigation  of  slotted  cylinders  has  led  to  the  following  con¬ 
clusions: 

(1)  For  the  case  of  resonant  coating  (a  «  top,  v  o),  the  admittance 
approach.. *  zero  as  the  thickness  of  the  coating  increases.  At  intermediate 
th' iknesses ,  the  susceptance  can  change  from  a  large  capacitative  to  a 
small  inductive  value  or  to  zero. 

(2)  For  all  practical  purposes,  the  admittance  of  identical  slots  on 
cylinders  and  ground  planes  is  equivalent  if  the  ctrcumference-to- 
wavelength  ratio  of  the  cylinder  is  greater  than  four  and  if  the  loss  angle 
of  the  plasma  lies  between  90°  and  l8o°,  and  the  magnitude  of  N  £  1. 

(3)  For  the  reentry  plasm  distribution  shewn  in  figure  7,  the  admit¬ 
tance  undergoes  pronounced  changes  when  the  propagating  frequency 
approaches  the  peak  plasm  frequency. 

(4)  The  admittance  of  apertures  on  large  coated  cylinders  (C  >  20) 
seems  to  be  sensitive  to  slight  changes  in  the  electrical  and  mechanical 
tolerances  of  the  coating.  As  such,  one  is  led  to  conclude  that  theoreti¬ 
cal  results  will  be  difficult  to  realize  experimentally  unless  electrically 
small  cylinders  are  used  (C  <  10). 

(5)  The  resonant  properties  of  the  cylindrical  gap  antenna  indicate 
that  it  may  be  useful  for  plasma  diagnostics.  The  presence  of  plasma  or 
any  change  of  antenna  parameters  brings  one  azimuthal  mode  into  prominence 
which  determines  both  the  ehaps  of  the  radiation  pattern  and  the  value  of 
the  antenna  input  impedance.  As  the  plasma  varies  in  density,  it  would  be 
expected  that  various  modes  would  ha  successively  excited;  therefore,  an' 
could  expect  the  pattern  and  the  input  impedance  to  change  accordingly. 
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APFEXDDC  I 

LIST  OF  SSKBOLS  US  ID  If?  THE  ANALYSIS  CCF  THE  AXIAL  SLOT 
radius  of  conducting  cylinder 


susceptance 

radial  distance  to  air-plasma  interface 
circumference  of  cylinder  in  wavelengths 
electric  field 

Fourier  transforms  of  electric  field 
amplitude  of  electric  field  at  aperture 
amplitude  of  incident  TE^  wave 
conductance 
magnetic  field 

Fourier  transforms  of  the  magnetic  field 
axial  mode  number 

length  of  slot  in  wavelengths  *  l/\ 

length  of  slot 

azimuthal  mode  number 

index  of  refraction 

power 

a  A 

Poynting  vector 
radial  mode  number 


r  222. 

Ty 


applied  potential  on  Blot 


74 


V  =  b/a 


w  width  of  slot 


Y 

7Ol 

p 

Pol 

r 

^mn 


«r 

n»S  . 

>«ol 

!*o 

v 

P/0>z 

0o 

CD 

®P 

Subscripts: 


admittance 

admittance  of  free  space 
waveguide  admittance  of  TB0]_  mode 
wave  number 

wave  number  for  the  SEqi  wave  vide  mode 
reflection  coefficient 
1 0 

Kronecker  delta  =  ( 

l1 

dielectric  coastant 
permittivity  of  free  space 
real  part  of  the  dielectric  constant 
transverse  components  within  a  rectangular  waveguide 
wavelength 
guide  wavelength 
permeability  of  free  space 
electron  collision  frequency 
cylindrical  coordinates 
angular  width  of  slot 
exciting  frequency 
plasma  frequency 

external 


a  n 
m  ■  n 


in  input 

l  per  unit  length 
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pla;*a  resonant 
vacuum 

vector  component  along  narrow  dimension  of  waveguide 
vector  components  along  the  three  principal  directions  in 
cylindrical  coordinates 
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APPENDIX  II 

ADMITTANCE  EXPRESSIONS  FOR  LONG  SLOTS  01,’  CYLINDERS 
The  pertinent  fields  in  regions  I  and  II  of  the  structure  shown  in 
figure  1  are: 

Hz1(p,0)  =  ^  Am[jja(pp)  +Bn£n<£%pj]e''5!I!?1  (II-1) 


H2E(p,0)  =•  Y1  CaHm(2)  (H-2> 


Uja-oo 


Y_  Aojja’Cpp)  +  BaBj/2) '  (pp)j  e“^  (n-3) 


E0D(p,0)  =  Yi  Ct'PjZ)' (fivQje-^  (H-4) 

*0  OT 


where  the  prime  denotes  differentiation  of  the  Bessel  and  Hanfcel  functions 
with  respect  to  the  argument. 

If  it  is  assumed  that  the  aperture  field  is  of  the  fora: 


then  the  houndary  conditions  at  p  =  a  and  p  =  b  give  the  following 
algebraic  relationships 


5) 


(II-6) 
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|  Anjja'(HCB)  +BffiHm(2)'(NCW)]  =  <^<2) ' (CW) 

^[jjKCW)  +  Bnft^NCW)]  =  C^2)(CV) 
Solving  the  above  equations  for  An  and  Bjj  gives: 


An  - 


NHm(g)(HCW)Hm(2),(CW)  - 

En(2)(CW)Tm  -  NHnteJ'tCV)!* 


^  (CW)  Jn'  (HCW)  -  HJm(HCW)^,^g| '  (CW) 

HHn(2)(HCW)Hn^2>'(CW)  -  { HCW) 

vhere 

Tm  =  jm' (ncw)%(2)  '  (rc)  -  jr'(Hc)Hn(2)'(KCw}| 

In  =  Jn(HCW)^2)'(HC)  -  Jm' ^0)^(2)  (new)  J 

Therefore: 

=  An[fn(KC)  +  BnBj2>(KC)] 


vhere 


=  -JVo 


SHy(g)'(CW)Vm  -  Hy(g)(CW)un 

y2)(Ch)^  -  NHn^’Ccw)!^/ 


%  «  JB(KC)%(2)'(RCW)  -  Jn'(KCW)^,(2)(Kc)| 
Vn  =  Je(RC)H_(2)(KCW)  -  Jm(RCW)^(2)(RC)  I 


(II-7) 

(II-S) 

(II-9) 

(11-10) 

(n-ii) 

(XI-12) 

(n-ij) 


Froa  Parsevai's  theoren,  the  external  admittance  per  unit  length  is  given 
by 


).  a 


.  S—  ,2mf~  MS* 

01  T0|2  .4.  |v„|  = 


(n-i4) 


And  the  substitution  of  equations  (II-5)  and  (H-I2)  into  equation  (H-l4) 
results  in  the  following  admittance  expression 

(n-15) 

2*a  Ba-B  p0\2  .  H^(2)  (cw)^ 


V  - 

Z. 


Tor  the  case  of  no  coating  (n>  n  1  or  b  =  a),  equation  (11-15) 


(U-J6) 


becosaes: 

r  t  .  -gj  f- 

%«'< 0  ’ 

It  has  been  shown1  that  at  plasma  resonance,  the  fields,  hence  the 
admittance,  are  indepent  of  the  azimuthal  coordinate,  0.  Therefore,  only 
the  m  =  0  mode  is  supported,  and  equation  ( 11-15)  reduces  to 

r,  ..£*  u»  „K(g)'( C»)V.  -  ■?(8W„'\  ( 

*~*°  \hJ2>(CW)T0  -  NH^W^j 
In  order  to  evaluate  the  above  equation,  the  following  expansions  are 


necessary: 


J0(x)  z  i  - 


Yo(x)  S  |  log  x 


J0'(x)  3-f 


V«  s  Jr 


(n-i8) 


to  give 


•\ 


To  -  ?['  ■  fl 

> 

Jc 

v0  S  |  log  w 

Therefore,  for  H  «*  0  and  V  >  1 


(H-19) 


r  JYq  _ Xytekw) _ 

2x11  [W2)  ’  (cvr)  +  H0(2){CVf)(w  -  lj] 


(n-20) 


To  analyse  the  inhomogeneous  plasm,  a  method  described  by  Svift^  la 
used.  Using  Swift's  notation,  the  axial  magnetic  field  with  the  plasm 
expressed  as  follows: 


(n-2i) 


where  the  ratio  **  tja(r)  +  jq^r)  is  numerically  evaluated  to  give 

Cj. 

solutions  of  tj!  and  Uj,  at  the  surface  of  the  conducting  cylinder, 
r  =  a.  The  details  of  the  technique  are  adequately  described  in  the  ref¬ 
erence  and  •will  not  be  repeated  here. 

At  r  =  a  the  boundary  condition  at  the  aperture  gives 


JSL 
2 jta 


sin' 


#•) 


C»[tn'(iO  +  Jum'(a)] 
__U= 


(US2) 


which  allovs  0^  aence  H^Ca,^)  to  be  cosqmted.  Therefore,  use  of 
Parseval's  theorem  gives  the  following  admittance  expression 
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REIATIOSSHIP  BETWEEN  INTERNAL  AND  EXTERNAL  AUCTIAKCE 
Computation  of  the  external  admittance  la  only  half  the  story  since 
the  external  admittance  as  such  is  not  a  measurable  quantity;  when  an 
experiment  is  performed,  the  input  admittance  is  referred  through  a  meas¬ 
urement  of  the  reflection  coefficient.  This  input  admittance  (or  reflec¬ 
tion  coefficient)  es  seen  in  the  feed  structure  muse  he  related  to  the 
external  admittance.  To  do  this,  however,  flow  conservation  iasaediately 
Inside  and  outside  the  slot  is  used. 

The  pertinent  fields  vhich  exist  in  the  vaveguide,  due  to  the  nave- 
guide  opening  onto  the  ground  plane,  are 

En  =»  E^e”J?olZcos^  jl  +  reJ2Pol“]  +  (h.t.)E  (m-la 


%  "  -Ooie"JPolZcos(f-)[l  -  reJ2polzJ  +  (h.t.)H  (m-lb 


Assuming  that  higher-order  terms  vanish,  i.e.,  (h.t.)g  »  (h.t.)g  -  0  the 

complex  power  flow  expressed  in  terms  of  the  fields  inside  the  slot  is 


pcq  -  |^|  2(1  +  r*)U  -  r)Ycl  (m-2 


But, 


4,0  v  1  +  r 


(m-3 


Therefore 
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But,  the  input  admittance  is  defined  as 


Y  a  5_.v  (l  -  r\ 

ln  S,,  ol\l  +  T) 


(m-5) 


Therefore,  the  following  relationship  exists  between  the  external  and  input 
admittance 


Xi5,  -  2| 


V*/ 


(m  -6) 


or,  in  normalized  fora, 

yin  -  ^  -  Sin  +  *>in  =  2(f)(^)yCy  (ro~7> 

Therefore  equation  HI-7  relates  a  calculated  normalised  external  admit¬ 
tance  v\  to  the  measurable  normalized  input  admittance  yln. 
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TABLE  II.-  TABLES  'T  MODAL  IMPEDANCE  OP  GAP  AKTHUfA 


m 

On* 

4ae  +  b*1 

Zg,  ohms 

0 

8.23 

27.1 

3.9  -  J12.7 

1 

7-06 

-.66 

53  +  J5-1 

2 

4.02 

48.8 

.6+  J7-7 

3 

1.11 

-107 

J3-5 

4 

.149 

-168 

J2.3 

5 

.011 

-229 

J2.0 

6 

-292 

J1.6 

7 

-355 

J1.4 

8 

-420 

i-9 

9 

-486 

J.8 

Waveguide  conductance  of  uncoated  axial  slot. 


Figure  15- -  Waveguide  conductance  c t  coated  axial  slot, 
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17.  ilL  RADIATION  PATTERN  OF  AN  ELECT  RIC  LINE  CURRENT 
Ei;c:X)SED  BY  AN  AXIALLY  SLOTTED  PLASMA  SHEATH 
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Department  of  Electric*.!  Engineering 
The  Radiation  Laboratory 
Ann  Arbor,  Michigan  48108 


ABSTRACT 

A  circular  uniform  plasma  sheath  of  complex  conductivity  is  enclosing  a 
unit  electric  line  source  at  the  origin.  The  plasma  sheath  has  am  infinite  axial 
slot  of  arbitrary  cross  section.  The  problem  has  been  reduced  m  a  Fredholm 
integral  equation  of  the  second  kind  in  which  the  unknown  funcUer.  is  the  electric 
current  >a  the  plasma.  An  approximate  solution  of  the  integral  equation  has  been 
obtained  in  the  case  of  a  locally  partially  transparent  plasma  sheatlu  On  the  basis 
of  this  solution  a  formula  for  the  radiation  pattern  is  derived.  As  an  example,  the 
radiation  patterns  are  computed  for  a  35°  wedge  slot  in  the  plasma  sheath  of 
three  different  thicknesses.  In  addition,  attenuation  calculations  have  been  per¬ 
formed  on  the  un-dotted  plasma  sheath. 
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I.  INTRODUCTION 

A  homogeneous  plasma  sheaJfe  of  inside  radii  a  and  outside  radii  b  is  taken 
to  enclose  a  unit  electric  an c  current  at  the  origin  of  the  coordieste  system  (r,  0, 2) 
as  shown  in  Fig.  1.  The  2-axis  coincides  with  the  line  source.  The  plasma  sheath 
may  cantata  an  axial  slot  of  area  Ag.  The  area  of  the  plasma  sheaih  we  have  indi¬ 
cated  by  Ap.  The  structure  is  uniform  in  the  2 -direction.  The  free  space  permit¬ 
tivity  and  the  permeability  nQ  is  taken  throughout,  including  the  plasma.  The 
electric  conduction  properties  of  the  plasma  (Rose  and  Clark,  1961)  are  accounted 
for  by  a  complex  conductivity 


<3  - 


9 

c  lT 

JLE_ 

l/*ju 


(1) 


where  u,  v  and  <jp  are  the  frequencies  of  Qe!d,  collision  and  picss&a  respectively. 

Since  the  primary  current  source  as  well  3s  the  plasma  snestk  does  not  de¬ 
pend  on  the  z-coordlnate,  it  follows  that  the  electric  ami  the  magnetic  fields  depend 
only  on  the  r  and  0  coordinates.  Furthermore,  the  electric  field  is  tangential  to  the 
Ase  source,  and  the  magnetic  field  is  funs  verse  to  it.  That  Is  to  say,  we  are 
dealing  only  with  the  electric  field  con*jonent  E^tr,  0)  and  the  magnetic  field  com¬ 
ponents  Hj.fr,  0)  and  H^(r,  0).  The  of  t  field  components  ue  zero.  If  we  let  the 
jplasina  DU  in  the  slot,  then  Ez(r,  0)-r- ^(r),  H^r,  0)  -*•  0,  and  K^tr,  0)  -*>  Hg(r) . 
Thus  we  have  only  one  competent  each  of  the  electric  and  the  magnetic  fields,  and 
they  depend  wily  on  the  radial  variable  r. 

In  this  report  we  are  concerned  with  Use  problem  of  finding  the  radiation  pat¬ 
tern  of  the  electric  line  source  enclosed  by  a  uniform  plasma  sheath  that  contains 
an  infinite  axial  slU.  7  hat  is  to  say,  -vc  wish  to  discover  the  0  variation  of  E.(r,  0). 
This  field  may  be  found  from  the  line  current  and  the  electric  current  induced  i. 
the  plasma.  The  latter  is  the  unknown  quantity  for  wnich  we  derive  ->  Fredholm 
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Integral  equation  of  the  second  kind  in  -Section  Bt.  We  may  solve  this  equation  only 
under  restrictions  which  make  the  solution  of  very  limited  practical  value.  Setting 
out  in  search  of  a  better  solution,  we  first  solve  the  unslotted  plasma  sheath  in  Sec¬ 
tion  HI.  Defining  the  sheath  current  with  a  slot  as  the  unslotted  sheath  current  plus 
a  correction  current,  we  are  able,  in  Section  IV,  to  derive  an  integral  equation  for 
the  correction  current.  We  use  the  unslotted  sheath  current  as  an  approximate  cur¬ 
rent  from  which  we  compute  the  radiation  pattern  for  the  slotted  plasma  sheath.  An 
approximate  solution  for  the  correction  current  serves  as  a  basis  for  checking  the 
reliability  of  the  radiation  pattern  calculations.  The  rationalized  MKS  system  of 
units  has  oeen  used  with  time  dependence  of  e^  understood.  In  section  V  we  pre¬ 
sent  numerical  results  and  conclusions  and  indicate  the  direction  of  the  future  work. 

H  REDUCING  THE  PROBLEM  TO  A  FREDHOLM  INTEGRAL  EQUATION  OF 
THE  SECOND  KIND 

The  unit  electric  line  source  induces  conduction  currents  iz{r',  0')  in  the 
plasma  sheath.  The  electric  field  (Harrington,  1961)  of  both  currents  is  given  by 

Ez(r,  jf>=  -  |  up0 [Hj>2)(k0r)+  j J  iz(r', {?-? |)dA-j  (2) 

AP 

(2)  _ 
where  H^  is  the  Hankel  function  of  the  second  kind,  k  and  }?-?'|  is  the 

distance  from  the  field  point  to  the  source  point.  We  assume  that  Ohm’s  law  holds 

in  the  plasma,  i.  e. 

iz(r,  5f)=CTE,(r,(3)  .  (3) 

We  are  neglecting  thereby  any  space  charge  affects  in  the  plasma.  Eliminating 
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E  (r,  0)  between  (2)  and  (3),  we  obtain  an  integral  equation  of  the  form 


z(v.t)+\w0a  Jj  yr'.jf'lH^tkJr-r-ljdA'-iu^aHf^r)  .  (4) 


We  recognize  that  this  is  a  Fredholm  integral  equation  of  the  second  kind  with  the 


kernel  (kJ  r-r  'I  ). 

ft  u 


We  put  the  integral  equation  In  a  more  convenient  form  by  defining 


Pouo 

s- 


4k_ 


(5) 


and 


i4i(p,5f)  =  izM) 

where  p  =  kQr,  and  hence 


(6) 


i(p»0)+T 


/f  w®!; 
S 


(2)(!p~p'l)dA'= 


(7) 


where  now  dA'=p'dp'd0' .  This  is  the  fundamental  equation  of  the  problem. 

The  integral  equation  may  be  solved  by  the  method  of  successive  approxi¬ 
mation  (Mikhiin,  1964}  provided  ihat 


1 


\r\<i 


(8) 


where 


!*  //  Jj  K2Wnf 


dAdA'  . 


(0) 


By  considering  unslotted  plasma  sheath  and  a  Hankel  asymptotic  form  for  wi 


rtain  from  (0;  that 
•'  B  —  2r  kgfb-a)  , 


(10 
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and  hence  from  (8)  that 


This  estimate  of  the  conditions  of  convergence  is  on  the  pessimistic  side.  Never¬ 
theless,  Eq.(ii)  gives  its  a  good  indication  that  the  successive  approximation  method 
of  solving  the  Fredholm  integral  equation  Lc  of  value  only  for  either  «ery  thin  plasma 
sheath  or  very  low  plasma  density. 

m  THE  CASE  OF  UNIFORM  PLASMA  SHEATH 


The  situation  is  as  in  Fig.  1,  with  the  slot  filled  in  by  the  plasma.  An  exact 
solution  to  this  problem  can  be  obtained  in  an  elementary  manner.  The  electric  Held 
may  be  expressed  in  the  form; 


Ez(r)=  -  \  upoH^)(k0r)+A1J0fK0r) , 

0  <  r  <  a. 

(12) 

'  A2Jo(kr)+B2No(kr). 

a<r <b , 

(lb) 

=  A3Ho>(kor) 

b  <  r , 

(14) 

where  JQ  and  NQ  are  Bessel  and  Neumann  functions  of  order  zero,  respectively,  and 

k  =  kjl  l-jo/(ue0) 

(15) 

Then  from  the  Maxwell's  equations  the  corresponding  magnetic  field  is  given  as 

0  <  r  <  a  , 

(16) 

a  <r  <b , 

(17) 

=j  — —  A„H^(k  r)  , 
up  3  1  o 

b  <  r , 

(18) 
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the  cylindrical  functions  bein'*  of  order  i. 

From  the  continuity  of  Ez(r)  and  Hz(r)  at  r=a  and  r=b,  we  obtain  four  algebraic 
equations.  By  solving  them  we  obtain  the  complex  amplitudes  of  the  waves: 

up 

V  (kcbg  + 
koNo(^)Hf)<koa)0^1{to)Ho,<^>b)-Mo{kb)Hi2><kob)]]  + 

kN^kamf  (k0a)  [k^ikblHf^^bJ-kJ^kblHf J  + 

kJi(ka)H^)(k0a([kNj(kb)H^\k0b)-k0N0(kbJH^(k0by|  . , 

(19) 

V  ^^Motfcbmf^W^OtbW^bg  ,  (20) 

o 


®2=  ~~~2  I 

1  2*aBk2  1 
o 

*  _  jt4i0 

A3  2  2  ,n 
j r  k^bD 


[kjjtkbjHf^w-v^kbjHf^bg 


D=  Ji(koa)Hl2)(kob)  [fodMlNoGcbJ-J^N^kaj]  +  - 

TT  l(koa)H!,2>(kob) [^(kMiyka)  J0(ka)Nl(kb)]  + 

Jo<koa>Hf  Kb)Do(kb)Ni(ka)'Ji(ka)No{klti  j  + 

(jr-  >2  J0(k0a’'Ho2>(k0b) [j^kajNjOdal-J^kblNjfkaj]  . 
The  current  Induced  in  the  uniform  plasma  sheath  i2(r)  is  given  by 
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i2(r)  srE2(r), 


a<  r  <  b 


-•  '  |I1 


•a 


_  ■  -*  ~S1 


1  s  ojA2J0(kr)fB2S0\Kr)|  . 

By  a  similar  definition  as  in  (G>  wo  introduce 

Up)  =  “  [A2Jo(kr)*B2Nc(kr)]  .  (25) 

ko 

The  electric  field  outside  the  plasma  shoe  nj  given  by  (14)  is  rewritten  in 


1  the  form 


The  first  factor  gives  the  cleetric  field  of  the  unit  eke  trie  line  source  in  the  free 
space,  the  second  factor  Introduces  the  effect  of  the  plasrua  sheath  on  the  radiation 
field.  vVe  define  the  latter  factor  as  the  radiation  pattern  fa.*  the  electric  line  source 
radiating  through  a  concentric  uniform  plasma  sheath,  i. 

PQ=  -  —  A-,  (27) 

o  UIIq  -i 

This  pattern  is  a  constant;  it  docs  not  depend  either  on  r  or  However,  it  is 
rather  difficult  to  calculate  because  of  the  involved  denominator  'a  A3.  However.,  in 
the  case  kt)a»  i  and  |ka| » 1,  P0  can  be  reduced  with  a  good  approximation  to  an 
elementary  form  by  using  the  Hankol  asymptotic  form  for  tha  cylindrical  functions. 
Taking  two  terms  in  the  iiankel  asymptotic  series,  we  obtain 
-jpUb-akkob-^kob)'1  \z\ 

?0  ~  - T - - -  (28) 

kojyd^a)  ^  j 


{where 


A 
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and 


and 


with 


dl=l  <  VTT (-‘-i;  )e->2«b-»>].u,{ko-  f . 

From  (1)  and  (15)  we  find  that  k»k'-Jk"  where 

***[{*♦  joap 

k"=ko[“  \  A^A2+B2'^ 

YX2 


A*  1-* 


l+‘ 


Y2  ' 


B< 


IV* 


XS 


ae. 


Y=-  . 
u 


(30) 

(31) 

(32) 

(33) 

(34) 


From  (28)  we  observe  that  the  amplitude  of  PG  la  determined  principally  by  the 
factor  fcxp[k"(b-a)J  as  would  also  be  the  corresponding  case  of  a  normally  incident 
plane  wave  radiation  through  a  plasma  slab. 

The  asymptotic  current  distribution  In  the  plasma  associated  with  (28)  is  given 

*  ,  rr  f 


i(r)=i 


!<i+ ^  + ib,a*  t?  5  ®',k<r"!0 


do^tskoa) 


(35) 


The  first  term  in  the  [  j  brackets  represents  outward  propagating  current  wave  and 
the  second  term  the  reflected  current  wave.  Again  the  exponentials  are  the  domi¬ 
nant  factors  in  this  formula. 
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^  -T  JJ  i(^2,<ip-/J'l)dA»  . 


However,  from  (38)  we  see  that  the  right-hand  side  of  (40)  reduces  to 

r  jf  i(pV2)({p-p«|)dA* 

As 

and  thus  .  * 

ic(p,0)«r  JJ Ic(p*.^)H;2)(|p-p*!)dA*=r  J J  i(p*)Hf(|p-p‘i)dA*  .(41) 


By  opening  a  slot  Ag  in  the  uniform  plasma  sheath,  we  change  the  sheath  current  by 
>c(p.0)  and  the  Integral  equation  (41)  defines  this  current.  Integral  equations  (?)and 
(41)  differ  only  in  their  right-hand  sides.  Thus  the  condition  (8)  on  the  convergence 
of  the  resolvent  applies  also  to  the  latter.  The  zeroth  order  approximation  to 
ic(p.0)  is  given  by 


hp.$)--7  J  j  i(p')H^\\p-p'\m'  . 


By  opening  the  slot  we  reduce  i(p)-0  in  the  slot  area  Ag  and  assuming  th»j  we  may 
neglect  ic(p,  0)  we  obtain  by  a  simple  transformation  from  (2)  and  (14) 


Vr*0" 


?>-  A3Ho  \  w0  J  j Kr'W^y  r-f’lWA' 


,  r>b  . 
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Using  the  asymptotic  fo*ms  for  the  Hankcl  functions,  we  obtain 

-4^0  Jdf  V  [V  j  j 1<r,)e  • 


r  »b 


(44) 


and  the  radiation  pattern  P  (jf).  consistent  with  the  earlier  definition  in  (27),  is  then 


Ps({&  *  PQ-  if  i(r*)e 


ik  r’cos^-^*) 


o 


dA*  . 


(45) 


It  is  clear  that  including  the  current  ij,  '(r,J?)  in  the  approximation  we  obtain  for 
the  radiation  pattern 


P,cW 


II 


f  jh  r*cos(jJ-0!) 


i(r')e 


1  f  f  {0)  jk  r'cos(jl-jl') 

dA'*  J  j  rV.ff'le  dA*  . 


(46) 


Equation  (45)  shows  us  that  opening  a  slot  amounts  to  adding  to  the  existing  cur¬ 
rent  its  negative,  i.  to  make  the  total  current  in  the  slot  area  zero.  For  over- 
dense  plasmas  ( P  { «  1  and  the  negative  current  in  the  slot  ar  ra  then  dominates  the 
radiation  pattern.  We  believe  fhat  (45)  gives  -  good  approximation  to  the  radiation 
pattern  whenever  i(r)  penetrates  into  the  sheath.  The  adding  of  i^(r,0  to  the 
radiation  pattern  PJfi  in  (46)  is  not  so  much  to  improve  the  approxi mntion  in  (45) 
os  to  provide  perhaps  s  means  of  feeling  out  its  bounds  of  validity.  Whenever 

P,(jf)  -  P^i  (47) 

one  caa  reasonably  say  that  (45)  gives  a  good  approximation  to  the  radiation  pattern. 

if  the  slot  is  a  wedge  slot  of  with*-  then  the  double  integral  in  (45)  may  be 
transformed  into  a  series,  and  wo  obtain  for  the  radiation  pattern 
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NT"  sin($J 
n=o  7  3 


<  1  I, 

n 


ft"  1,2,3,... 
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It  is  clear  that  the  pattern  is  an  even  function  in  j5.  The  variable  is  measured 
from  the  center  of  the  slot. 


V  NUMERICAL  RESULTS  AND  CONCLUSIONS 

In  Figs.  2,  3  and  4  we  present  the  radiation  through  a  uniform  plasma  sheath 
for  Y-0. 01,  0. 1  and  1,0  respectively.  With  Y-0. 01  the  plasma  is  only  slightly 
lessy.  With  Y«l.  0  the  plasma  is  very  lossy.  The  inside  radius  of  the  plasma  sheath 
has  been  kept  at  one  free  space  wavelength  X^  and  the  sheath  thickness  has  been 
taken  up  to  3A  .  The  plasma  parameter  X  is  varied  in  steps  tram  0.25  tc  i.  0. 

V 

In  Figs.  2  and  3  the  calculations  have  been  carried  oat  using  the  exact  formula 
(27)  for  all  X  values  .accept  X-2.  For  X-2  calc  lation  we  used  the  asymptotic  form 
/28)  with  dj  set  to  aero.  The  same  asy  nptotic  form  was  used  for  all  calculations 
In  Fig,  4,  The  asymptotic  form  for  PQ  gave  excellent  agreement  with  the  exact  re¬ 
sults  in  Figs.  2  and  3  where  it  was  applicable,  and  surprisingly  little  disagreement, 
even  for  cases  (e.  g.  Y-0. 01,  X=l)  where  the  arguments  of  the  cy  liir'-r  functions 
were  far  too  small  for  the  Hankel  asymptotic  expansion  to  be  valid. 

The  processes  that  inhibit  radiation  through  the  sheath. are  reflections,  ab¬ 
sorptive  attenuation,  and  reactive  attenuation  by  the  plasma.  The  attenuation  of  the 
radiation  for  0  <  X<  (/? approx. )  is  due  to  reflections  and  absorptive  attenuation,  and 
for  X> (& approx. )  due  mainly  to  both  reactive  and  absorptive  attenuation.  In 
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each  case  we  may  identify  conditions  where  only  one  process  is  primarily  causing 
the  attenuation,  Tl>us  for  X  <f?  and  Y=0.01  we  sec,  in  Fig.  2,  the  radiation  graphs 
oscillating  because  of  the  reflections.  As  Y  increases  to  0.1  in  Fig.  3,  the  ampli¬ 
tude  of  the  oscillations  decreases,  and  for  Y=i  in  Fig.  4,  the  oscillations  are  practi¬ 
cally  gone;  but  now  the  graphs  have  assumed  a  downward  slope.  This  is  caused  by 
the  absorptive  attc.  uation  in  the  plasma  taking  over  and  eliminating  the  reflections. 
Also,  we  see  that  we  get  quite  a  bit  more  attenuation  for  very  lossy  plasma  (Y=  1) 
than  for  very  low  loss  plasma  (Y=0.01).  For  the  case  of  X  >/?  the  opposite  is  true; 
low  loss  plasma  causes  more  attenuation  than  a  very  lossy  plasma.  This  is  due  to 
the  effectiveness  of  the  inactive  attenuation. 

In  Fig.  5,  we  show  a  36°  wedge  slot  Li  the  plasma  sheath,  and  in  Figs.  G  to  14 

2 

we  present  the  radiation  power  pattern  [Ps{0)|  calculated  from  (48)  for  the  sheath 
thickness  of  0.  oXq,  1.  OXq  and  2. 0  Aq,  respectively.  The  inside  sheath  radius  is 
kept  at  1.  OAq.  The  graphs  have  been  arranged  in  sets  of  three;  Y=l.  0,  0. 1  and 
0. 01,  and  X  is  a  parameter  in  each  set.  Since  the  power  patterns  are  even  functions 
in  0,  we  have  plotted  them  only  from  0°  to  180°.  By  the  constant  lines  we  show  the 
attenuation  in  the  power  radiated  before  the  plasma  slot  is  introduced,  i.  e.  the 
radiation  power  pattern  for  the  unslotted  plasma  sheath.  These  lines  are  obtained 
from  Figs.  2  to  4.  In  some  cases  the  attenuation  is  so  large  that  it  is  off  scale. 
These  values  we  will  give  in  a  note.  Asa  way  of  calibration  we  mention  that  the 
radiation  power  pattern  of  unit  amplitude  would  correspond  to  the  electric  line 
source  radiating  in  the  free  space. 

In  Figs.  6,  7  and  8  we  present  the  set  Y=l.  By  opening  a  slot  in  this  very 
lossy  plasma  we  have,  for  all  practical  purposes,  re-established  radiation  in  the 
forward  direction,  i. «.  0=0,  independent  of  the  sheath  thickness.  Outside  of  the 
front  lobe  there  is  no  significant  recovery  of  the  signal  strength.  Increased  slot 
width,  the  calculations  will  show,  narrows  the  front  lobe  with  some  additional  gain 


THE  UNIVERSITY  OF  MICHIGAN 
;>S23-4-T 


£3 


144 


THE  UNIVERSITY  OF  MICHIGAN 

5825-4-T 


ifl  the  signal  strength  for  0-0.  If  we  halve  the  slot  width,  t*  i  forward  lobe  is  con¬ 
siderably  reduced  with  no  benefit  for  the  rest  of  the  pattern. 

&  Figs,  9,  10  &hd  11  we  present  the  set  Y=0. 1.  By  opening  a  slot  ir,  tliis 
sheath  of  "medium"  lossy  Plasma  we  observe  that  for  X  <  0. 70?  the  formerly  unl- 
form  pattern  breaks  into  oscillations.  The  radiation  is  enhanced  slightly  in  the 
forward  dirtsetlen.  The  oscillating  character,  for  the  patt  ;ra  persists  as  X  ie  in~ 
creased,  i.  e.  a  dominant  forward  mode  is  not  formed  Li  the  range  of  the  validity  of 
s&prttdm&tton,  in  contrast  to  the  case  of  the  very  lossy  plasma  sheath.  Sections  oi 
graphs  for  which  the  approximations  are  expected  to  break  down  are  shown  in  the 
dotted  line.  For  X  >  0. 70?  the  slot  improves  the  radiation  very  substantially  as 
compared  to  the  chealh  without  a  slot. 

In  Fige.  12,  13  and  14  we  present  the  set  Y=D.  01.  i.e.  a  very  "low"  loss 
plasma  sheath.  Tbs  comments  for  this  one  are  substantially  the  same  as  for  the 
previous  case. 

The  sets  cf  data  have  been  presented  as  they  were  calculated.  A  dominating 
forward  lobe  formed  for  the  Y=1  set.  It  came  as  a  surprise  that  for  Y=0. 1  and  0.01 
seta  ao  aach  dominance  was  found.  However,  there  is  no  doubt  that  a  dominant  fbr- 
!  ward  lobe  would  develop  as  X  is  increased  beyond  $P.  Only  (48)  is  not  then  a  valid 
i  approximation  in  tne  case  of  our  thick  sheath. 

The  net  bod  we  have  presented  for  calculating  the  radiation  pattern  for  the 

slotted  piasm  sheath  makes  use  cf  the  electric  currents  in  the  plasma  before  the 

slot  Is  opened.  It  is  an  approximate  method  which  is  expected  to  give  the  radiation 

pattern  with  good  accuracy  only  if  the  electric  currents  sufficiently  penetrate  the 

plasma  sheath.  If  the  attenuation  of  the  unslotted  plasma  sheath  is  more  than  20  db, 
o 

i.e.  {Pc|  <  0.01,  then  we  are  approaching  the  bounds  of  validity  of  the  approxi¬ 
mation  for  the  pattern  as  a  whole.  However,  one  expects  that  the  forward  lobe  will 
be  given  to  a  good  order  of  approximation  even  tor  somewhat  larger  attenuation 
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than  20  db,  as  was  shown  by  calculations  performed  in  (46) 

If  the  attenuation  of  the  unslotted  plasma  sheath  is  such  that  the  radiation 
through  it  is  practically  insignificant,  we  then  find  that  the  electric  currents  pene¬ 
trate  the  plasma  only  a  small  fraction  of  the  free  space  wavelength.  If  we  open  a 
slot  in  the  plasma  under  these  conditions,  the  slot  in  fact  becomes  a  waveguide 
which  connects  the  source  with  the  outside  space.  For  suitable  slot  boundaries  one 
;  can  show  that  a  finite  number  of  modal  fields  can  propagate  in  the  slot  waveguide, 
j  These  modes  may  be  used  as  a  basis  for  calculating  the  radiation  pattern.  Work  is 
i  in  progress  to  exploit  this  waveguide  approach  for  these  cases  where  the  integral 
equation  method  becomes  too  cumbersome.  This  work  will  be  reported  as  Part  n 
of  this  study. 

The  two  significant  conclusions  indicated  by  this  siudy  are: 

1)  Opening  a  slot  in  the  plasma  sheath  re-establishes  radiation  in  the  forward 
direction  to  a  substantial  level. 

2)  The  plasma  slot  should  be  at  least  one-half  wavelength  wide  at  its  narrowest 
Gross  section  when  the  electric  field  is  polarized  tangential  to  the  slot  walls.  In¬ 
creasing  the  slot  width  beyond  one-half  wavelength  brings  only  diminishing  returns 
as  far  as  the  radiation  enhancement  in  the  forward  direction  is  concerned. 
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7.  THE  EFFECT  OF  A  FINITE  CONICAL  PLASMA  SHEATH 
ON  THE  RADIATION  FIELD  OF  A  DIFOLE* 

U.  C.  Pric'aore-Brown 

Aerospace  Corporation,  Los  Angeles,  California 


ABSTRACT 


The  plasma  sheath  that  forms  about  a  reentry  vehicle  is 
simulated  by  a  cone  of  finite  slant  height  and  negligible  thick¬ 
ness.  This  cone  can  support  surface  currents  that  arc  induced 
in  it  by  a  Hertzian  dipole  antonna  situated  on  its  axis.  The 
radiation  pattern  of  such  a  structure  at  wavelengths  long  com¬ 
pared  to  the  slant  height  is  calculated  using  a  Lobcdev-Kontorovich 
integral  transform  togothor  with  the  Wiener-Hopf  technique.  In 
the  opposite  limit,  the  radiation  pattern  is  obcained  by  a  pertur¬ 
bation  technique.  Plots  are  presented  showing  the  effects  of 
sheath  strength,  slant  height,  and  antonna  frequency  on  the 
patterns. 


This  research  was  supported  by  the  U.S.  Air  Force  under 
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I.  INTliOUUCTIOX 

1  2 

The  present  study  is  one  of  a  series  *  undertaken  to  gain  a  better 
understanding  of  the  role  of  diffraction  effects  or  antenna  pattern 
distortion  in  the  so-called  blackout  problem.  This  term  is  commonly 
used  to  describe  a  sharp  reduction  in  the  effective  transmission  of  infor¬ 
mation  by  radio  between  a  reentry  vehicle  and  a  ground  station.  This 
reduction  is  observed  as  the  vehicle  enter  he  earth's  atmosphere  and 
becomes  enveloped  in  a  sheath  of  partially  ionized  gas  or  plasma.  It  is 
due  partly  to  the  ability  of  this  plasma  sheath  to  reflect  electromagnetic 
waves  and  so  alter  the  radiation  impedance  of  the  antenna;  partly  to  the 
sheath's  ability  to  spread  waves  into  new  directions,  thus  causing  antenna 
pattern  distortion;  and  partly  to  the  fact  that  the  plasma  responds 
differently  to  different  frequency  components,  which  results  in  phase 
distortion.  There  will  also  be  an  additional  damping  of  the  signal  by 
ohmic  losses  in  the  sheath,  but  this  effect  should  be  small  in  the  case 
of  plasma  sheaths  that  are  very  thin  compared  to  a  wavelength. 

As  in  previous  studies,  we  shall  restrict  ourselves  to  a  point-dipole 
antenna  that  is  driven  by  a  constant-current  source  at  a  single  frequency. 
Thus  we  shall  be  mainly  concerned  with  antenna  pattern  distortion.  Two 
previous  studies  considered  the  case  of  a  magnetic  dipole  antenna1  and 
an  electric  dipole  antenna^  located  axisymmctrically  within  a  thin  conical 
plasma  sheath  of  infinite  extent.  The  plasma  was  represented  by  a  current 
sheot  whose  strength  was  calculated  on  the  basis  of  a  cold  electron  plasma. 
The  results  of  these  two  cases  showed  that  the  electric  dipole  antenna 
generates  a  strong  surface  wave  that  carries  most  of  the  radiated  onergy 
out  along  the  surface  of  the  conical  sheath,  whereas  no  such  wave  is 
produced  by  the  magnetic  dipole  antenna.  In  fact,  it  is  clear  that  the 
magnetic  dipole  induces  purely  azimuthal  currents  in  the  conical  sheath 
that  are  confined  to  the  vicinity  of  the  source  for  long  wavelengths. 


149 


whereas  the  currents  induced  by  the  electric  dipole  nm  up  and  down  along 
the  generators  of  the  cone.  These  facts  suggest  that  the  adaption  of  a 
mere  realistic  model,  consisting  of  a  cone  of  finite  slant  height,  may 
profoundly  modify  the  radiation  pattern  in  the  electric  dipole  case, 
particularly  if  the  wavelength  is  comparable  to,  or  longer  than,  the  cone’s 
slant  height.  The  present  work  is  concerned  witn  this  case. 

The  proolea  is  solved  by  Wiener-Hopf  techniques  using  a  Lebedev- 

Kontorovich  integral  transfora.  It  reduces  to  an  infinte  set  of  algebraic 

equations  which  can  be  solved  by  iteration  in  the  lev-frequency  limit.  This 

method  has  already  been  employed  in  the  problem  of  deteraining  the  scattering 

of  plane  acoustic  waves  by  a  rigid  disc3  and  the  impedance  of  a  biconical 
4 

antenna.  In  the  opposite  limit,  when  the  slant  height  is  long  compared 
to  a  wavelength  but  still  finite,  the  problem  is  solved  by  a  perturbation 
technique  based  on  the  known  infinite-cone  solution.  These  two  techniques 
are  used  to  calculate  radiation  patterns  that  show  the  Banner  in  which  the 
distortion  fades  and  the  radiation  resistance  gradually  resumes  its  free- 
spaco  value  as  the  slant-height  is  reducec  from  infinity  to  zero. 
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II.  ANALYSIS 


He  consider  tho  radiation  field  of  an  eleoontary  electric  Hertzian 
dipolo  that  is  located  on  the  axis  of  a  thin  partially-conducting 
conical  sheath  of  finite  slant  height  b.  He  adopt  spherical 
coordinates  centered  at  the  apex  of  tho  cone, which  then  becosas  the 
surface  6  «■  60,  0  <  r  (  b  (Fig.  1).  The  dipole  source 
is  at  T  •  a,  8  •  0  and  is  oriented  along  the  cone's  axis  so  that  the 
problea  is  axisyaaetric.  Since  then  all  currents  are  purely  radial, 
it  is  easily  verified  froa  Maxwoll's  equations  that  we  can  express 
the  fields  in  teras  of  a  single  scalar  potential  n(r,o)  such  that 


£  -  V  «  V  x  (rU) 

Using  HXS  units  and  suppressing  a  cosascn  tine  factor  »  we  then 
find  that  B(r,9)  satisfies  the  nave  equation 

(V2  ♦  k2)  n(r,8)  - 
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where  k  ■  u  uQc0<  The  nonvanishing  field  coeponents  are 


Er  "  ’  F  sIHS-  h  (sin8  w)]1 
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Now  the  current  density^  •>  J(r,0)  cr  is  Bade  up  of  the  antenna 
current  and  the  current  induced  in  the  sheath. 

J(r,8)  -  -iup  6(r  -  a)  ♦  o(r,6)Er  (6) 

Here  p  is  the  dipole  esossent  of  the  antenna  and  o(r,8)  tho  conductivity 
of  the  sheath;  o(r,0)  is  taken  to  be  zero  everywhere  except  at  tho 
sheath  itsolf, where  it  is  assumed  to  be  infinite  in  such  a  way  as 
to  a;Io^  induced  sheet  currents.  Thus,  we  write  it  as  follows 


o(r,8)  *  (iuc0r)  |u{b  -  r)«(0  -  8^  (7) 

o 

in  tenas  of  u  diwnsionless  sheath  paraaatcr  g,  the  step  function 
U(b  -  r)  (■  1  for  r  <  b,  G  for  r  >  b)  and  the  delta  function  6(8  -  0o). 

It  can  bs  shewn  that,  in  terss  of  plasaa  paroetera,2 


*d 

2r 


sinO 


o  ♦  ivcJ 


(*) 


where  snd  vc  sre  plasaa  and  electron  collision  frequencies,  respectively, 

and  d  is  tho  sheath  thictosss.  In  the  calculations  we  shall  assuas  that 
*) 

v  e  0  and  c.“  (d/r)  *  const., which  Bakes  z  a  positive  zeal  constant. 

V  P 

Tho  representation  of  a  thin  plasaa  slab  by  a  current  sheet  in  this  way 
was  first  used  by  tfsit.5 

Tho  solution  to  Eq.  (2)  can  be  represented  formally  as  the  solution 
to  an  integral  equation 
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Here  we  have  substituted  for  J  froa  Eqs.  (6)  and  (7)j  R  is  the  source-to¬ 
re  ceivcr  distance  shown  in  Fig.  l.d's  «  /sine  and  we  have 

suppressed  a  factor  (p/4*c0a)  for  convenience. 

A.  LONG  SAVE  LENGTH  LIMIT 

Although  the  integral  equation  (9)  can  bo  inverted  exactly  in  the  case 
of  an  infinite  cone  (b  ■  «)2,  it  does  not  appear  to  be  possible  to  do  so  in 
general.  In  this  part  we  present  an  analysis  which  yields  the  solution 
in  the  Unit  of  long  wavelengths,  1  »  b  >  a.  Part  B  will  be  devoted  to 
the  case  b  »  X  »  a. 


Wo  introduce  the  notation 


*(r)  -  (ra)1/2^n(r,0o)  ,  *Q(r)  »  (ra)1/2X^^ 
n0(r,e)  ■  »  (ra)-1/2^4/(r,e|a,0) 


(10) 


(11) 


(12) 


in  tons  of  which  2q.  (9)  bewau-s 
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which  con  be  deduced  fro a  the  well-known  expansion  for  tho  potential 
of  a  ring  source^ 

^tr,a|r',e*)  -  £  (2n»l)  PR(ccse)  Pn(cos8')  lnvl/2(7r<)  Kn+1/2(Yr>) 
n=*0 

(17) 


This  expansion  converges  for  nil  y  (real  or  imaginary).  Here  8<t  $> 
stand  respectively  for  tho  losser  and  the  greater  of  8  and  9',  and 
similarly  for  ir  and  r>.  To  establish  the  equivalence  of  Eqs.  (16)  and 
(17),  notice  that  as  |uj  •*  »  (u*±  0,  ±  1,  2,...) 


jsoc  unP^j^^ose^  Pu.i/2C-co»V  I  '  |w|“*  <r-p£-|lB(u)|Ce>-9<)] 


(18) 


(19) 


(see  Appendix  A,  Eqs,  (A-l),  (A-2)  and  (A-iJ)}  so  Eq  (16)  converges 
absolutely,  provided  we  assume  (temporarily)  that  y  •>  -ik  is  real. 
Furthermore  the  integration  path  can  bo  closed  by  a  semicircle  at 

t 

infinity  in  the  right-hand  plane  if  r  >  r.  Xf  this  is  not  the  case, 
one  can  appoal  to  the  relation  [sec  Eq.  (A-3)]  which  exists  between 
the  hyperbolic  Bessel  functions  to  show  that  for  any  even  function 
f(y)  [such  that  f(w)  ■  -f(-u)] 


(20) 


Consequently  r  and  r'  can  be  freely  interchanged  ir.  Eq.  (16)  provided  the 
integration  path  is  taken  along  the  imaginary  axis  (c  *  ).  Finally 
use  of  Cauchy's  theorem  together  with  Eq.  (A-ll)  leads  to  Eq.  (17). 


15 


From  Eq.  (16)  wo  deduce  also  tlie  integral  representations  for 
G(rjr')  and  *0(r) 


/+i“ 
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udu  ,.,2 
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(d2  -  l/4)Pu.l/2t-C0SVIu^)Ku(Yr’) 


(|c|  <  1/2) 


(22) 


llorc  wo  have  used  tho  fact  that  («C  ♦  H2  -  l/4)Py_1^2(cos8)  •  0. 

Lot  us  now  introduce  Eq.  (21)  into  Eq,  (IS)  and  reverse  the  order  of  the 
U  and  r  integration.  The  integration  over  r  can  then  be  carried  out 
immediately  by  virtue  of  Eq.  (A~8) ,  and  we  find 


m 

f  ^KvCYr)G(r|r') 
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f  ~  Kv(yr)V0(r) 


cosv* 


Pv-l/2<'C0SVKv^a) 


(24) 


Similarly, 


liero  the  integration  over  u  was  performed  by  noticing  that  because 
the  integrand  is  odd  in  u  the  integration  path  can  be  closed  on  the 
iasAginary  axis.  Introducing  those  expressions  into  Eq.  (IS)  we  find 


/ <Kr)  lyCr*}  ~  ♦  /  *(-)  ~ 

o  b 


/\(r)  Kv(Yr)  & 


D(v)  -  cosvi  ♦  g(v2-  t)  Pw-1/2Ccos6o)  Pv-1/2('c0*#o) 


This^is  an  intogral  equation  involving  two  unknown  functions  of  v, 
viz.,  f  #(r)Kv(Yr)y*1  dr  and  j£"$(r)  ^(yrjr*1  or.  In  general  we 
must  resort  to  complox  variable  techniques  to  solve  it,  but  first  let  us 
consider  tho  limiting  case  of  an  infinite  cone  (b  ■+■  »). 


2.  Infinite  Cone  Limit 

Denoting  this  case  by  the  subscript  <•  we  have  from  Eqs.  (24)  and  (2S) 


„  2 
jt.M  s„(Yrt 


Introducing  Eq.  (16)  into  Eq.  (13)  and  making  use  of  Eq.  (27)  we  find 


n.(r,e)  -  yr.e)  *  (ra)'1'2  i  /*  g  *_(r')  •  £  /^ggs*  P^Ccos «<) 

M  ^ 


*  Pu-l/2C*cose>)  V*r>  yv*’)_ 


n0(r.e)  ♦  (TzyUhsj1  Jg-  4^rVl/2Cc05®<)PM.l/2('cos9>J 


‘  Vl/S^V-V^  K„Clft) 


n  .  <  ,..,-1/2  Z'1',  r  1  1  1  P1,-1/2("M<)Pu.1/2('“"») 

o<  ■  ‘ lra)  Jj  “^  nsrl  — p-^rssg - 


*  Xu(y»)  s  Cta) 

*»  V 


ilere  0<l  6>  refer  to  0  and  6q.  Coaparison  of  the  first  integral  with  Eqs. 
(12)  and  (16)  raveals  that  it  exactly  cancels  no(r,8).  Thus  ws  are  left 
with 

nr,.,.  .~.W/2  .  L  f% d.U  .  Pu-l/2tCOSB<)  P„-l/2<-C9S«>J  . 

1  ;  T  JJGT  r”2'fci.sBo) - VyT)  V*a) 


This  solution  has  already  been  studied  in  sose  detail  in  a  previous 
paper.  There  it  was  shown  that  it  can  be  expressed  as  a  suocation  over 
the  roots  of  D(u)  »  0,  these  being  the  orders  of  those  Irgendre  functions 


PU-l/2^cos&J  £^at  C0BP0S«  the  natural  eigenfunctions  for  the  infinite 
cone,  Orto  result  of  that  -rtud/  that  we  nust  sake  use  of  in  the  p retail t 
paper  concerns  the  behavior  of  the  solution  in  the  vicinity  of  the  apex 
of  the  cone,  namely  the  fact  that  ^(r)  *  const  rc  as  r  -*  0  where 
c  «  |Ro{v0l|  and  *v9  is  that  pair  of  seros  of  Eq.  (26)  which  is  newest 
the  imaginary  axis.  If  wa  allow  a  small  amount  of  dissipation  in  the 
sheath  so  that  v£  >  0  and  consequently  Ira{g)  <  0,  then  it  was  shown" 
that  c  >  0,  and,  in  fact,  c  -*  0  uniformly  u  v^-  9.  Since  9(n)  is 
even  in  v,  this  means  that  thero  is  a  strip  |Re{y) |  <  e  containing  the 
imaginary  axis  where  D(w)  has  no  zeros.  Now  the  behavior  of  $(r) 
near  the  apex  must  be  like  that  of  from  which  we  conclude  that  the 
first  integral  in  Eq.  (25)  converges  within  this  strip.  Thus  Eq.  (25) 
is  to  be  considered  as  defined  within  the  strip  |Ra{v}|  <  c  ■  |Re{v0)|. 

#e  now  take  up  the  task  of  finding  a  solution  to  Eq.  (27)  in 
the  general  case  of  a  cone  of  finite  slant  height  b.  To  do  this  we 
regard  Eq.  (27)  as  a  Wiener-Hopf  problem  and  try  to  express  each  cf  the 
terms  appearing  in  it  as  sums  or  products  of  analytic  functions  of 
v  that  are  holoaorphic  »n  overlapping  half-spaces  of  the  complex 
v  plane  and  behave  suitably  at  infinity. 


S,  Kioner-Hopf  Decomposition 


Consider  first  the  second  term  in  Eq.  (25).  Since  Kv(yr)  is  an  even 
integral  function  of  v,  it  is  clear  that  the  same  is  true  of  this  term. 
Using  the  asymptotic  form  of  *v(yr)  given  in  Eq.  (A-5)  we  find  as  |v| —  » 


(29) 
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Thus  the  behavior  at  infinity  is  determined  by  that  of  t(r)  sear 

7 

r  ■  b,  and  this  is  given  by  the  edge  condition 

❖(r)  *  (r  -  b)“P  as  r-»b*  (30) 


where  p  {  1/2,  the  equality  holding  for  an  infinitely  conducting  cone. 
Thus  we  write 


/* 


«t)Ku(rr)  |i 


1.  ins  V 


£irb/2) 


(3D 


where  g* (v)  is  an  integral  function  of  v.  Pros  Eqs-  '29)  #a A 

(30)  wo  deduce  that  g* (v)  •—  0  as  |v|  -»  •  in  any  right-hand  hslf-piane, 

Ko{v}  » 


In  a  siailar  way  considor 
<-(r)Kw(Yr)  |I-  ^r(v)TC-v) 


« 

(Yb/2)v 


*(r)vI.w(Yr) 


dr 

r 


S'  (V) 


♦  JEjtSJL 

(Tb/2)‘W 


- j  *£r)vlv(yr) 

S'  (-v) 


£32) 


As  pointed  out  in  the  last  section,  this  expression  is  holonorphic  in 
|Ko(v}|  >  c  and  hence,  by  sycaetry  the  first  and  second  tents  on  the 
right  are  holoaorphic  respectively  in  the  overlapping  half-planes 
Ro(v)  >  -c  and  Re{v)  <  c  .  Furthermore  assisting  the  edge  condition 
•%(r)  *v  const  as  r  ->  b'  we  deduce  that 


£*{v)  »  0(v*X)  in  Re{v)  <  c 


£33) 
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It  rcnains  to  consider  the  coefficient  D(v)/cosvx.  He  should  liko 
to  bo  able  to  write  it  as  a' ratio  of  functions  holonorphic  respectively 
in  left  and  right  half-planes.  Now  it  has  been  demonstrated  that  for 
real  g  U(v)  has  simple  zeros  along  the  roal  axis  which  are  in  one-one 
correspondence  with  those  of  cosvs.  In  fact,  denoting  the  zeros  of  D(v) 
by  v#  it  was  shown  that 

±vs  ■  s  ♦  1/2  ♦  «  ,  s  -  1,2,...  (34) 


where  0  s  £  <  1.  In  addition  D(v)  has  a  pair  of  zeros  situated  near 
the  imaginary  axis  that  wo  denote  by  ±  vQ.  Thus  D(v)/cosvx  has  simple 
zoros  at  v  o  ±  vs,  s  ■  0,  1,  2,  ...,  and  simple  poles  »t  v  •  ±  (>  ♦  1/2), 
s  ••  1,2,...  .  (Tho  zoros  of  D(v)  at  v  ■  ±  1/2  are  cancelled  by  those 
of  cosvtj)  These  facts  suggest  the  possibility  of  writing  D(v)/cosv» 
in  terms  of  infinite  products 


U(v) 

■  ■  ■!  i  ■ 

COSVT 


D(0}^1  ♦  v/v^^Q*  (v)  •  (l  -  v/voj  /Cf*  (-v) 


(35) 


ri*  «  r  i  ♦  v/vs  1  *  r  v  4s 

M  •“jpnsTcrrw  "  I1'  Trrrrr *  rya'*  172% qr 


(36) 


It  is  clear  that  //  (v)  converges  and  is  non-zero  for  Re(v)  >  -3/2  and, 
furthermore, that  it  tends  to  zero  as  |vj  —  -  in  this  half-plane. 

To  determine  its  rate  of  fall-off  at  infinity  consider  the  asymptotic 
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no  of  D(v),  given  by  using  Eq,  (A10)  in  Eq.  (26) 

°(v)  *  iTsf no"  v  co*Cveft  "  */4)  cos(vt  -  vS  -  v/4) 


v  tends  to  infinity  along  rays  Baking  a  finite  angle  with  the  real 


Lm 


cgjv*  *  sine 


|v|  1 vl  -»  -,  arg(v)  *  0,« 


ron  this  we  can  deduce  that 


*(v)  •  0(jv|'1/2)  in  Re(v)  >  -3/2 


Ho  new  attempt  to  separate  Eq.  (25)  by  multiplying  by  (yb/2)v/r(v) 

>nd  making  use  of  Eqs.  (24),  (31),  (32), and  (35). 

(l  -  [«••(»)  •  ffer-  CYh/i}2*  «•-(-.)] 

•  [d(0)  U  .  ./.0)^(»)]-‘  m.  0 


The  first  tern  in  this  expression  is  holonorphic  and  Odvl’^2)  in 
Re(v)  <  c.  Tho  retaining  terns  will  only  tend  to  zero  in  the  right-hand 
piano  provided  a  >  b,  but  thoy  still  siss  being  holoaorphic  there 
bocauso  of  tho  simple  polos  o£$*(-v),  r(-v)  and  1/cos  vw  in  this 
region.  Thus  it  is  not  possible  to  transfors  Eq.  (25)  into  a  finite 
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algebraic  equation  which  can  be  solved  by  an  appeal  to  LiouvilXe's 
the ore*.  Instead  ono  oust  subtract  out  the  residues  at  these  poles  so 
as  to  obtain  art  algebraic  equation  with  an  infinite  nunbor  of  tens. 

The  restriction  a  >  b  according  to  which  the  source  is  excluded  fro* 
the  interior  of  the  conical  sheath  was  necessitated  by  the  fact  that  the 
last  ter*  in  Eq.  (40)  goes  as  |u|(b/a)vcxp[-|!*{v)0ojj  as  |v|—  *» 
in  the  right-hand  plena.  Fortunately,  however,  the  radial  cosponent 
of  the  electric  field,  Ey  ■  -r*1  j£-d,  satisfies  a  reciprocity  principle 
(see  Appendix  B).  Consequently  we  can  use  the  results  cf  the  present 
analysis  applied  to  the  case  r  <  b  <  a  to  obtain  as  well  the  fields 
in  the  reciprocal  case  a  <  b  <  r,  which  is  the  me  that  interests  us. 

Vo  rcaark  parenthetically  that  if  it  is  desired  to  have  a  fonalis* 
which  applies  to  all  cases  (including  the  regaining  one  a  <  b,  r  <  b), 
this  can  bo  done  at  the  cost  of  soae  extra  complication  by  subtracting 
the  quantity 

cS  f  ♦-«  *vCrr)  T-  *  jT ♦-(r)  MYr)  TT 

fro*  both  sides  of  Eq.  (25)  to  yield 

hS/* Vm  V”)  r  *  f V»e  T-  -fcS-  if 

o  b  J  b 

(41) 

Whore  VW  •  ♦(r)  -  ♦.(r) .  When  this  equation  is  sultiplied  through  by 
(tb/2)v/r(v)  the  last  ter*  will  dearly  go  to  zero  as  jv|  -*  -  in  the 
right-hand  plane  for  all  b. 


and  va  have  pvt  (1  -  v/vo)£t/*(-v)  «  L£v).  Also 


L(n  ♦  1/2)  -  v  ^  J1®  1/2  {[v-(n  ♦  l/2)]L(v)| 

Non  K*  (v)  ,  '4*(v)  are  clearly  "plus"  functions  (holoaorphio  and 
tending  to  zero  in  Rc{v)  r  -c)  while  2^J,  2^  and  2^  are  ‘‘ainus"  functions 
Thus  we  write  E<i.  '«!)  in  the  fora 

J(v)  -  L(v)  £v)  +  £“  ♦  2l"  ♦  r' 


*  -!i*(v)  -f*(v)/L£-v)  -  Q*(v) 


and  conclude  fwa  Uouvi  lie's  theorea  that  J(v)  •  0.  lienee  we  find 
«a  reintroducing  Eqs.  (44)  through  (46) 

»  f  .  ._,2n»l 

h(viT(w)  *E  a^*(-n  -  l/2)  S V  tl  *  Y/If 


To  solve  this  equation  we  assume  a  power  series  in  {yb/2}. 


S' M  «  ^  iin(v)  (Tb/2)n*i/2 


(52) 


n«l 


and  then  equate  successive  powers  of  yb/2  to  transfora  it  into  the 
following  infinite  systen  of  equations 


n/2-1 

lm  eom  .  y  _ 


T^I73T  £n-2a-lC*5‘1/2) 


(n-n/2  fa_ 

*  Z  ^  En-2at‘a) 

E»1 


n  ,  ^ 

V“(nVl>2)  n  * 


(S3) 


Here  the  integral  parts  of  the  upper  liaits  (n/3-1  and  {n/3-1/2  are 
to  be  understood.  Solving  the  systes  of  equations  by  iteration  and 
then  substituting  into  fcq.  (32)  we  find 


i 
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l«(v)  (vj 


C1 


W2)Z/2  *  vT§77  (tb/2)5^2 


<- 


C3 

^77T 


bi  ci 

tTOTF-lV 


(yb/2) 


7/2 


♦ 


C4 

7377 


(Yb/2)9/2 


(S4) 


Jtotico  that  (v)  is  essentially  the  L-K  transfoni  of  the  induced  current 
in  t.ie  shon'ch.  We  turn  next  to  the  task  of  expressing  the  radiation  field 
in  teras  of  g~  (v) 


4.  The  Fomal  Solution 

Hie  scalar  potential  n(r,8)  is  given  foraally  by  Eq.  (15).  If  in  this 
expression  no  introduce  Eq.  (16)  for^(r,e|r\eo)  and  make  uue  of  £q.  (32) 
we  obtain 


n(r,0)  *  "o(r*e)  *  (ra)”1/2  hr  f  C&  Pu.i/2(co*e<)pw.V2(’cose>)Ip<7r) 

-i» 


Sr  ♦  -rJ-~'L  £>*) 


(Yb/2) 


(yb/2) 


(55) 


167 


Making  use  of  Eqs.  (35),  (42),  and  (SO),  we  car.  rewrite  thia 
rosult  as  follows 

n(r.O)  -  n0(r,e)  ♦  (r^j"1/2  D(0 )/  PW.1  /2^co#0<)  Vl/2("6C!*®>)  * 

-i» 

•  I  (vr)  L(-y)  JM-  (k(u)  *  «*(„)]  (S6) 

“  (Yb/2)1*  LT*  J 


in  which  fora  it  i<  clear  that  the  y-ccntour  can  be  closed  in  the  right- 
hand  plant  providod  r  s  b.  The  only  poles  are  at  ji  -  n  ♦  1/2  due  to 
S"(y)  and  at  u  •  u$  where  D(ys)  •  0.  The  contribution  fro«  those  at 
n  ♦  1/2  is 

(r.)'1'2  fr  (-2wi»[0)  j  ^g-  W*  '„(“*»  WM  ' 

n«l 

•  «-“• 1/21  (Z'irfo  {»*“;‘«)  [tv‘"'1/2) 

On  saking  uso  of  Eqs.  (46)  and  (49)  it  is  easily  verified  that  this  reduces 
to  the  series  expansion  for  — Fo(r,0)  as  given  in  Eqs.  (12)  and  (17),  Thus 
the  contribution  fro*  these  poles  exactly  cancels  the  first  tens  on  the  right 
of  Eq.  (55).  As  a  result  the  total  solution  is  given  by  the  integral  in 
Eq,  (55)  evaluated  a»  a  sun  of  residues  at  the  poles  y  ■  yf.  Putting 
cosy*  -  D (y)  /  [d  (0)  L (y)  L  (»y)  ]  in  Eq.  (55)  we  have  then 

n(r,0)  ■  -(ra)"1/2  2g  0(0)  ^  5rCT  Py-l/2^COs9<*  Py-1/2*"CO*0>*  Vyr*  * 

y 


■  L(-y)  l(y)  6~M 

(Yb/2)“ 


(57) 


(r  t  b  <  e) 
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where  the  sisaation  is  over  the  roots  of  D(iO  «  0.  Notice  that  there  is 
no  contribution  froa  the  second  tens  in  tho  integral  in  Eq.  (55)  since 
i*(ns)  £"(-us)  *  o. 

As  a  partial  check  of  this  result  consider  the  limit  g  ■*  0.  In 
this  limit  it  is  easily  seen  from  Eq.  (26)  that  for  the  sth  root 

«s  s  ns  *  s  -  1/2  — ►  &  s(s  ♦  1)  | Ps(co»co)]2 
and  from  Eq.  ($1)  (since  an  •»  0) 

t(us)  £'(es)  —  «“*  c.(Yb/2)s+1/2 

where  c_  is  given  in  Eq.  (4!)).  Introducing  these  results  into  Eq.  (57) 
it  is  readily  verified  that  H(r,8)  —*■  no(r,e)  as  it  should. 

Me  now  have  a  solution  valid  for  the  receiver  inside  the  cone  and 
the  source  outside.  To  obtain  the  solution  in  the  reverse  case— source 
inside  and  receiver  outside— we  appeal  to  the  reciprocity  theorem 
(Appendix  B)  which  allows  us  formally  to  interchange  source  and  receiver 
coordinates  in  the  radial  co^xment  of  the  electric  field.  According 
to  Eq.  (5)  this  component  can  be  computed  froa  Eq.  (57)  as  follows 

5,  -  -r-l^B(r,e)  -  -  JL-  Craf3/22gD(0)£  Pu_1/2(cos9<). 

u 

*  pu-i/2(”cose>)  *„(**>  "  Ltp) 

"  1/2  *  p  (vb/2)p 


(rS  b<>) 


(58) 


and  wa  have  reinstated  the  factor 


(p/4*e0a)  that  was  suppressed  after  Eq.  (9).  In  this  expression  the 
receiver  coordinates  are  contained  in  the  ten  Py-l/2(cose)Ip(lfr) 
an  i'ii  ■■  >urco  coordinates  in  L(y)  £“(y)  and  hence  ultimately  in  the 
ccv.:fi«ussnts  cn  [see  Eqs.  (49}  and  (54)].  It  now  becomes  clear  that 
in  order  to  use  reciprocity  we  should  have  fc.aulated  the  problem 
for  a  ring  source  instead  of  a  point  source.  The  reader  will  readily 
verify  that  if  we  had  done  so  the  result  would  have  been  to  replace  the 
source  torn  Kn(ya)  in  the  coefficients  cn  in  Eq.  (49)  by  ?n(cosa)  Kn(r«), 
where  a  is  the  angle  subtended  by  the  ring  source.  (Notice  that  as 
the  ring  shrinks  down  to  a  point,  a  0  and  Pn(cosa)  -*•  1.)  Once  this 
is  done  we  can  formally  make  the  substitutions  r  ♦  a,  0  +  a  +  0  in 
the  recoiver  coordinates  in  Eq.  (58)  and  a  ♦  r,  a  ♦  B  in  the  source 
coordinates.  Finally  we  can  operate  on  the  result  by  -r  X”1  to 
obtain  the  scalar  potential 

n(r,6)  -  (ra)-i/22gD(0)  ]T  A(u)  X'1  [L(y)T(w)]  (S») 


where 


This  fora  is  now  valid  in  the  range  a  <  b  <  r.  Here  L(y)<£r (u)  is  given 
by  tho  right-hand  side  of  Eq.  (54)  so  X’1  eventually  operates  on  the 
cR  (which  are  now  functions  of  r  and  6).  Thus  if  we  introduce  the 
coefficients 


cn(r»9)  „  1  Pn(coa9)  W 

n(n  *  ^  r(n»l/2)  [J:+(n»l/2)/v0J^ (n>l/2) 


vs  800  that  the  ten  e’M]  in  Eq,  (59)  is  given  by  the 

right-hand  side  of  Eq.  (54)  with  cR  replaced  by^. 

S.  The  Far-Fiold  Radiation  Pattern 

To  calculate  the  field  at  large  distances  froa  the  cone  we  use 
the  asymptotic  form  Kn+1^2(yr)  *  (t/2tr)^2  c“Yr  for  the  radial 
functions  appearing  in  the  'c^.  This  allows  a  radial  factor  r”1  e*Yr 
to  be  taken  out  of  the  egression  for  H(r,f}  which  we  then  write 
formally  as  a  power  series  in  yb 

«o.«  ■  (jjSjt)  2  r»(,)  <*>" 


The  angular  coefficients  Fn(i)  can  be  found  by  equating  successive 
powers  of  yb  in 

£>„(»  w)” .  fcf2  «*»;£>> 

n«l  v 

Substituting  the  right-hand  side  of  Eq,  (54)  for  *^T3(1.(p#)  (us)  1 
with  v  replaced  by  and  cn  by  ^  given  in  Eq.  (61)  we  find  for  the 
first  few  coefficients 

V9)  -  •  (r  2  r^r  pit««V  P1(“S8) 
f2(*>  » -  i  (if2  rphr  2  r^rp2(«»»V  p2<“*9> 


(62) 
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1 


i 


1/2 


r3Cor  -  .  ^  {Zj  TXyfay  2  ~T777  P3(c«e0)  P3(COS8) 

“  fe(rj  rr-f/2)'  rfrif  2  pi<cos6> 

p4(si  “  -  i5r  (rj  rpW  S  ir^r  V^’V  p4 («■»•) 

*w(r)  rrr5737  rfnrj  5W  p2(co3o) 
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Tho  intensity,  as  given  by  the  radial  eswponont  of  the  Poynting 

flux. 


sr  *  7 

bocostos  in  tho  far  field 


Introducing  fcq.  (62)  for  It  wo  write  this  as  follows 


(64) 


in  terms  of  an  angular  distribution  factor 


2 

F(0)  .  _Ly  Y  F'(6){-ikb)n  (6S) 

(ka)2  i-J  n 

Ho re  we  have  reintroduced  k  ■  iy. 

In  the  absence  of  the  sheath,  i.e.,  as  g  ♦  0  it  is  readily  verified 

that 

5>nCe)CYb)n  —  (s/2ta)1/2  2  (2n»l)Pn(cos0)In+1/2(va) 

which  is  thn  expansion  of  the  plane  wave  expCyacosS)  in  spherical  waves.6 
Thus  it  follows  that  in  this  limit  F(6)  sin2e,  which  is  the  free  dipole 
distribution. 

6.  Nuaorical  Work 

To  compute  the  angular  distribution  factor  F(6)  given  in  Hq.  (65)  it 
is  necessary  to  ovaluato  the  coefficients  F^(S) ,  which  are  the  derivatives 
of  tho  expressions  given  in  Eq.  (63).  These  involve  the  infinite  products 
L(-p)  and  the  infinito  sums  S(p)  »  £MlO/(y-p)  for  various  integer  and 
half«intoger  p's.  The  first  30  roots  u  of  the  dispersion  equation  D(u)  »  0 
together  with  the  derivatives  D  (u#)  and  the  corresponding  Legendre 
functions  ^.1^2 ‘cos8'  were  obtained  on  an  electronic  coeputer  by  the  method 
outlined  in  the  report  of  Rtf.  2,  The  infinite  products  L(«p)  *  (1  ♦  p/u^/T (p) 
witheC*(p)  given  in  Eq.  (36)  were  then  computed  ss  follows 


In  tents  of  these  quantities  the  infinito  sum  S{p)  ■  ^A(u)/(y-p) 
wre  easily  computed  using  £q,(60)  for  A(u),  Notice  in  this  expression 


Ctb/2)u  y  V7  L 


f'a/2)2 
y ♦  T 


HS*-'  3 

.-.-4  ;§ 
S?Ji-  1 


'  f  I 
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Coaputed  values  of  S(p)  and  L(-p)  corresponding  to  the  case  ka  ■  1/3, 
kb  •  1  and  for  sheath  parancters  g  -  0.1,  1  and  10  are  given  in  Tabic  1. 

The  series  Eq.  (65)  unfortunately  docs  not  converge  for  larger  values  of  kb 
so  that  ita  uso  in  computing  is  limited  to  situations  where  the  wavelength 
is  large  compared  to  the  slart  height  of  the  cone.  For  such  leng  wavolengths 
the  radiation  pattern  turns  out  to  be  practically  indistinguishable  froa 
that  of  the  free  dipole,  although  its  level  can  be  considerably  different. 

In  fact,  for  the  cases  considorod,  the  radiated  power  levels  referred  to 
the  free  dipole  level  were  coaputed  to  be  14. S  db,  9.4  db,  and  -19  db  for 
t  ■  0.1,  1  and  10,  respectively. 

5.  LONG  COKE  LIMIT 

In  this  part  we  sake  use  of  the  known2  current  distribution  along  on 
infinite  cone  to  obtain  an  approxiaate  radiation  pattern  for  a  cone  of  finite 
length.  It  is  expected  on  physical  grounds  that  such  an  approximation  should 
bo  good  for  cones  whose  slant  height  is  long  coapared  to  a  wavelength. 

Let  us  return  to  the  fundamental  intogral  equation  uq,  (9)  and  rewrite 
it  in  an  obvious  way  as  follows 

*<*.*»  ■  MM>  -  f  fix'/  &  tfife’.y  (66) 

•'b  J  r-r 

*  Sr  M  * 


Here  is  the  solution  for  b  «  **.  It  is  clear  that  if  we  introduce 

()  ir.  place  of  a£u(r’,90)  under  the  integral  sign  in  Eq.  (66) 
tiic  result  will  be  to  yield  ar.  approximation  to  JI(r,8)  valid  for  large  b. 
Physically  we  are  replacing  the  actual  current  distribution  on  the  cone 
by  the  one  which  obtains  when  the  cono  is  infinitely  long.  This  approxi¬ 
mation  could  itself  be  introduced  into  the  right-hand  side  of  Eq.  (66)  to 
yield  a  second  approximation,  and  so  on,  but  we  shall  not  go  beyond  the 
first  one. 

nM(r,S)  can  be  obtained  from  Eq.  (28)  by  closing  the  contour  at  infinity 
and  evaluating  the  resulting  integral  as  a  sub  of  residues  at  the  zeros  of 
D(u)  (cf.  Kef.  2).  The  result  is 


n.(r,6)  -  2*ik£  -~- 
v  D  (v) 


Pv-l/2(c05e<)Pv-l/2(-C0se>> 


Pv-1/2{C0$V 


Here  we  have  used  Eqs.  (A-l)  and  (A-2)  to  express  the  modified 
Bessel  functions  in  terms  of  the  spherical  ones.  The  summation  in  Eq.  (67) 
is  over  the  roots  of  D(v)  *  0  whero  D(v)  is  defined  in  Eq.  (26).  Fro*  Eq.  (67) 
wo  obtain  the  expression 


;Cn.(r\V  ~2.it  £ 


which  is  to  be  substituted  for  ^n(r' ,80)  under  the  integral  sign  in 
Eq.  (06) . 


In  the  far  field  (r  ■•  •}  the  Green's  function  in  Eq.  (66)  can  be 
siaplificd 

r2*  exp{ikjr-r*  |)  ikr  /2*  «  ) 

I  " . ^ 1  ■•  —  -jr*'  /  cxpMkr  [sin8sin8  cos(6~8  )*cos8cos8  Jjjj’ 

J  i&£  I  J, i 


1S1 

■  ~ —  exp(-ikr'cos6coi8')  JQ(kr  sinBsine')  (69) 


Sisdlarly  the  llanksl  functions  appearing  in  Eqs.  (67)  and  (68)  can  be  replaced 
by  the  asynptotic  values 


hU>  (krl  ..  i*^1^  CL 
v-1/ 2l*r-  ~  1  To T 


This  mist  clearly  be  done  in  Eq.  (67),  and  in  Eq.  (68)  it  is  justified  since 
we  are  assisting  that  b  is  large.  With  these  ■edifications  we  introduce 
Eqs.  (67),  (68),  and  (69)  into  Eq.  (66)  to  obtain 
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The  far- field  radiation  pattens  bccoao*  {cf.  Eq.  (64)]. 
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where  F^(6)  is  the  derivative  mth  respect  to  6  of  the  expression  in 
braces  in  Eq.  (70). 
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Here  the  discontinuity  in  £^(8)  has  beer,  exhibited  explicitly  with  the 
help  o£  the  Wsbcr-Schafheitlin  intogra’o.  It  is  readily  verified 


that  ^v(0)  itsolf  is  continuous 


•  -T 

C  Jfl" 


cosvv 


7  SI 


2y 

'v 

sine 


-  0 


T9  shew  this  we  have  used  the  Wronsxian  relation  [Eq.  (A-12)]  and  the 
fact  that  v  satiriios  0{v)  >0. 

Note  that  £^(8)  do pends  only  on  8q  and  on  kb  and  is  independent  of 
the  sheath  pa roaster  g  and  the  source  location  a.  It  has  been  evaluated 
numerically  for  kb  »  20,  6q  *  10*,  and  the  results  are  given  in  Table  2. 

A  table  of  the  Legendre  functions  Py, ,j(cos6)  for  g  •  0.1,  1,  end  13  has 
bees  furnished  in  the  previous  report. 
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III.  RADIATION  PATTERNS 

Figure  2  is  a  plot  of  the  radiation  pattern  of  the  dipole  source  in  an 

infinite  conical  shoath  for  various  values  of  the  sheath  parameter  gc  The 

dipole  is  driven  at  a  frequency  corresponding  to  ka  ■  1/3,  and  the  cone  angle 

is  0  *  10*.  flic  dotted  curve  labelled  g  »  0  is  the  frce-space  pattern, 

°  2 

I  *  sin  6.  The  intensity  lovcls  are  given  in  decibels  referred  to  that  of 
the  froc-space  pattern  *t  DO*,  that  is,  to  unity.  Thus  the  quantity  plotted 
is  10  log1QF(9)  where  F(0)  is  given  in  Eq.  (6S)  or  Eq.  (71),  It  is  clear 
that  the  presence  of  the  sheath  leads  to  enhancod  power  output  (increased 
radiation  resistance)  for  low  valuos  of  g  and  to  a  surface  wave  which  is 
especially  pronounced  at  intermediate  values.  Those  effects  are  discussed 
at  length  in  the  previous  report,* 

Figure  3  shows  the  offect  of  increasing  frequency  on  the  radiation 
pattern  labelled  g  ■  0.1  of  the  previous  £i jure.  The  incipient  formation  of 
lobes  at  the  higher  frequency  is  evident.  The  radiation  pattern  is  fairly 
insensitive  to  frequency  with  the  level  proportional  to  (ka)'  for  (ka)  S 1/3, 

Figures  4  and  5  exhibit  the  effect  of  finite  slant  height  on  the  radiation 
patterns.  Thoy  have  been  drawn  for  shoath  parameters  g  ■  0.1  and  g  ■  10, 
respectively,  and  the  same  cone  angle,  10*.  The  curves  labelled  kb  *  1 
(corrospc  .ding  to  1  ■  2*b)  nave,  of  course,  been  computed  by  the  scheme 
described  in  Part  A  while  the  perturbation  scheme  of  Part  B  was  used  to  obtain 
the  curves  labollod  kb  ■  20.  In  this  last  case  the  slant-height  is  only  about 
3  wavelengths  long,  b  ■  (10/»)X,  so  the  numerical  values  say  be  somewhat  in 
error  as  they  ere  basod  on  the  first  iera  only  of  a  perturbation  expansion 
in  X/b, 

These  plots  indicate  that  truncating  the  cone  har  a  pronounced  effect 
on  the  radiation.  At  a  slant  height  of  threo  wavelengths  the  surface  wave 
is  sharply  reduced,  and  thoro  is  some  change  in  the  total  power  output. 
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At  a  slant  hoight  of  one-sixth  of  a  wavelength  the  surface  wave  has 
entirely  disappeared  although  the  total  power  output  is  still  markedly 
different  fron  that  corresponding  to  the  bare  dipole.  It  is  interesting  that 
the  power  output  remains  high  for  tho  weakly  reflecting  sheath  (g  *  0.1) 
but  it  is  greatly  depressed  by  tho  strongly  reflecting  one. 

The  tremendous  variation  in  intensity  radiated  at  different  angles 
necessitates  the  use  of  a  decibel  or  logarithmic  scale  to  display  the 
entire  pattern.  However,  the  distortions  this  introduces  are  misleading. 
Accordingly,  different  portions  of  the  radiation  patterns  of  Fig.  4  have 
been  reproduced  on  a  linear  scale  in  Fig.  6, 


ACKNOWLEDGMENT 


The  author  is  indebted  to  Dr.  R.  H.  Huddlestons  and  Professor  A,  Bafios 
for  valuable  discussions  of  this  problem. 


REFERENCES 


1.  D.  C.  Pridaore-Brown,  J»  Hath.  $  Phys.  43,  199-271  '1964), 

2.  A.  Bafios,  Jr.,  R.Ii.  Huddlestone,  O.C.  Pridaore-Brown  and 
E.C.  Taylor,  Aerospace  Technical  Report  TDR-269(S220-10)-2, 

(October  1964);  also  J.  Math.  5  Phys.  44,  (1965). 

3.  A.  Loitner  and.C.  P. -Wells,  IRE  Trans,  on  Antennas  and  Propagation, 
AP-4,  637-640  (1956);  also  B.  Noble,  Methods  Based  on  the  wiener- 
ilopf  Technique,  Pergaaon  Press,  London  ^1958]),  p.  213. 

4.  J.  A.  Meier  and  A.  Leitner,  Oesterroichisches  Ingenleur-Archlv, 

13.  141-151  (1959). 

5.  J.  R.  Wait,  Appl,  Sci.  Res.  8,  Section  B,  397-417  (1960). 

6.  P.  Morse  and  li.  Feshbach,  Methods  of  Theoretical  Physics,  McGraw-Hill, 

New  York  (1953),  p.  1466.  . 

7.  C.  J.  Bouwkaap,  Rep.  Pro2r.  Phys.  17.  35-100  (1SS4). 


I 


i 


182 


APPENDIX  A.  TABLE  OF  FORMULAE 


The  codified  or  hyperbolic  Bessel  functions  Iy(x),  Ky(x)  are 
related  to  the  cylindrical  and  spherical  Bessel  functions  by 

y»>  •  •>„«*>  •  (2r)I/J  w!i,)  <*•» 

*„(*>•  i"'*  J  (1i)  ■  l“*‘  (iji)1/2  hj»/2  (lit)  (»-2) 

and  to  each  other  by 

yx)  -  -  $  uT(u)r(-p)II.y(x)  -  Iw(x}]  (A-3) 

Note  that  they  are  integral  functions  of  y  and  that  Ky(x)  is  even  in  y. 
They  have  the  asys^totic  forms  s 

ly(x)  <w  (2«x)'1/2  ex  ;  Ky(x)  »  (i/2x)1/2  o'*  (A-4) 

for  x  »  |y |  and 


for  |e|  »  x. 
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(A- 9) 


and  having  the  asymptotic  behavior 


Vl/2(C059)  ' 


/  2  y\ 

l*|  u|  sine) 


cos  (v6-t/4) 


(A-IO) 


as  juf  ■  provided  0  t  6  <  i. 


When  y  ■  a  ♦  1/2 t  n  an  integer 


pn(-x)  -  c-r  pfl(s) 


(A-1X) 


and  otherwise  iy  -  n  ♦  1/2} 
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In  this  Appendix  we  verify  that  the  usual  reciprocity  principle  holds 
in  the  pwsenco  of  current  sheets. 

Consider  two  source  distributions  of  current  Jj (r)  and  (rl  and 
their  corresponding  fields  Ej,  llj,  Hj  ?n  a  region  V  enclosed  by  surfaces 
3.  Then  the  relation 

/%  x  ax  -  tj  X  H2)  •  JC  ds  .y'fEj  *  J2  -  E2  •  jp  dV 

follows  directly  froa  Maxwell's  equations  on  the  assusptiot;  that  the  sources 
have  the  saas  httsonic  tiae  dependence.  If  now  nf-  take  the  current  sources 
to  bo  radial  and  to  have  the  fora 

ii  -  ^  «Ct  -  It) 

£  “  £2J 

and  if  we  consider  the  surfaces  S  to  consist  of  the  inner  and  outer  surfaces 
of  the  conical  sheath  together  with  the  surface  at  infinity,  then  this 
equation  reduces  to 

“  Elr^2J  "  E2r&l} 


/  [W  “  Wl  ^ 

am.  «-  J 


m 


Hhcro  tlso  integration  is  now  taken  over  just  one  surface  of  the  cone,  and 
IH^]  denotes  tho  juap  in  across  this  surface.  Icposition  of  the  juap 
boundary  conditions  across  this  surf  ice  then  causes  the  integral  to  vanish. 
Thus  wa  conclude  that  the  radial  coeyonent  of  the  electric  field  is  syeactric 
in  source  and  receiver  coordinates. 
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VI.  PROPAGATION  OF  TRANSIENT  SIGNALS  IN  PLASMAS 
C.  T.  Case  and  R.  E.  Haskell 

Air  Force  Cambridge  Research  Laboratories 
Bedford,  Massachusetts 

INTRODUCTION 

Most  analyses  of  wave  propagation  in  plasmas  are  based  on  the 
assumption  of  a  monochromatic  wave.  Although  no  real  signals  are 
truly  monochromatic,  this  assumption  is  often  a  very  good  one  for 
many  applications.  However,  with  the  advent  of  very  narrow  pulses 
and  extremely  fast  rise-times,  the  frequency  spectra  of  certain  signals 
arc  broad  enough  for  a  large  number  of  frequency  components  to  inter¬ 
act  significantly  with  the  plasma.  1 

The  dispersive  nature  of  a  losslees,  isotropic  plasma  may  be 
seen  from  the  simple  dispersion  equation 

k  =  -  -  n2  (1) 

c 

•where  k  is  the  wave  number,  w  is  the  angular  frequency,  c  is  the  free- 
space  velocity  of  light,  and  II  is  the  plasma  frequency.  This  dispersion 
relation  is  plotted  in  Figure  1,  where  it  is  observed  that  no  propagation 
occurs  fo»  frequencies  below  the  plasma  frequency.  As  the  frequency 
increases  above  the  plasma  frequency  the  phase  velocity  decreases 


T 


from  inanity  toward  c  while  the  group  velocity  increases  from  zero 
toward  c. 


The  response  of  a  nonmonoehromatic  signal  propagating  in  a  { 

3 

dispersive  media  can  in  principle  be  determined  by  a  superposition  j 

of  all  of  the  Fourier  frequency  components,  each  of  which  propagates  j 

according  to  Figure  1.  If  the  frequency  spectrum  is  sufficiently  narrow  1 

i 

and  peaked  about  a  carrier  frequency  iu  ,  then  the  signal  may  be  treat- 

o  { 

ed  as  a  wave  packet  and  the  wave  number  k  given  by  (1)  may  be  expanded  | 

in  a  Taylor  series  about  <so«  This  technique  has  been  widely  used  to 

j 

determine  the  overall  degradation  of  pulsed  signals  propagating  in  dis-  | 

2-5 

persive  media.  However,  this  anal/sis  is  not  valid  if  the  frequency 

a 

spectrum  is  very  oracd  and  thus  it  can  not  be  used  for  very  short  pulses, 
for  describing  the  initial  arrival  of  the  signal,  or  for  analyzing  with  no 
carrier  frequency.  In  this  paper  these  special  cases  will  be  discussed 
by  solving  the  initial  value  problem  of  waves  propagating  in  plasmas 
using  Laplace  transform  techniques 

FORMULATION  OF  THE  PROBLEM 


Maxwell’s  equations  which  describe  the  propagation  of  electromag¬ 
netic  waves  in  an  isotropic  plasma  are 


curl  t  =  -  n 


curl#  =  J+  tg 


a* 


3£ 


(2) 


Z  at 


Z{a,  z)  =  E{s,  0)e 


(6 


where  E(s,  0)  is  the  Laplace  transform  of  £(t,0).  The  time  response 
C(t,  r)  is  then  obtained  by  taking  the  inverse  of  {6}.  That  is. 


qt.  z)  = 


St 

£{s,  0)e 

v 

s 


5n3 

5+S 


ds. 


(7 


If  collisions  can  be  neglected  then  (7)  reduces  to 


t%  z)  = 


—  S 

Zri  J 


st 


E(s.  0)e 


f  Vsa 


+  D 


ds. 


{« 


We  shall  investigate  various  solutions  of  (8)  in  the  following  two 
sections.  Equation  {7}  will  then  be  investigated  in  order  to  determine 
the  roie  which  collisions  play  m  the  propagation  of  transient  signals  it 
plasmas.  Finally,  reflections  from  a  plasma-vacuum  interface  wiU 


be  considered. 
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NUMERICAL  SOLUTIONS  FOR  SHORT  DISPERSION  LENGTHS 

There  are  very  few  initial  transforms  E{s,  0)  for  which  the 
inverse  transform  (8)  is  known.  One  exception  is  a  unit  step  mag¬ 
netic  field  at  z  =  0  for  which  H{s,0)  =  1/  s.  From  Maxwell's  equations 


so  that 


E(s,  z) 


H(s,  z) 


E(s,0)  = 


(9) 


Equation  (8)  may  then  be  written  as 


I  0  0  <  t  <  z/  c 

£U.  r  1 


(10) 


where  J  (x)  is  the  zero  order  Bessel  function  of  the  first  kind.  Note 
o  ' 

that  no  signal  arrives  prior  to  the  time  t  3  z/c. 

The  solution  (10)  is  a  damped  oscillatory  function  wnose  frequency 
approaches  the  plasma  frequency  for  large  times.  The  argument  of  the 
Bessel  function  in  (10)  may  be  written  as 

$  3  Ii  \J t2  -  z'/c‘  j  (11) 


If  we  define  the  "instantaneous  frequency"  of  the  oecillation  by 
Uj  -  d$/ut,  then  from  (11)  we  find  that 


n 


t.  i 

V 


I  -  z2/c*  t2 


(12) 


Thus  for  early  times  (t£  z/c)  the  instantaneous  frequency  is  very  high 
and  approaches  the  plasma  frequency  for  very  long  times.  The  instantaneous 
frequency  also  increases  for  larger  plasma  frequencies  and  longer  dispersion 
lengths. 

For  other  forms  of  E(s,  0),  the  inverse  transform  (8)  can  not  be 

found  from  tables  and  alternate  methods  for  evaluating  the  integral  in  (8) 

must  be  found.  For  the  particular  case  of  a  step  sine  wave  input  the  solution 

6-8 

can  be  written  in  terms  of  a  sum  of  Bessel  functions  ,  which  can  be 

9-11 


expressed  in  terms  of  tabulated  Lommel  functions, 


Integral  solutions  which  may  be  evaluated  numerically  for  short 
dispersion  lengths  may  bn  obtained  either  by  contour  integration  or  by 
using  the  convolution  theorem.**'  ^  In  this  section  we  shall  present 
results  obtained  by  use  of  the  convolution  theorem. 

The  time  response  given  by  (8)  may  be  written 


where  U(t  -  z/c)  is  the  unit  step  function,  we  can,  by  using  the  convolution 
theorem,  write  (13)  as 


If  a  unit  step  E  field  is  turned  on  at  z  =  O.then  (14)  becomes 


{t.  z}  =  1 

- 


*  ji  nr  **/'z 


dt*  U(t  -  z/c) 


Equation  (15)  has  been  evaluated  numerically  and  the  results  are  shown 
in  Figure  2.  Note  that  the  argument  of  J  ^  in  (15}  is  die  same  as  (11)  and 
that  the  instantaneous  frequency”  shown  in  Figure  2  is  in  accord  with  the 
discussion  following  (12). 

If  a  unit  step-carrier  sine  wave  is  turned  on  at  z»  0  then  (14)  can  be 
written  in  normalized  form  as 


£,  (Tf  1r\)  *  Xm  j  « 


,i(T  -  n) 


T  -iU 


-  Pqe 


f  •"‘Jl  PV» 

1  — ~rr*~T~ 


du  U(T  -  q) 


/  u  -  q 


where 


r  »  M  t 
o 

q  a  wez/c 

P  s  E/>0  . 


205 


20? 


If  we  substitute  (18)  into  (14)  and  use  (17)  we  obtain 


The  results  of  numerical  integration  of  (19)  are  shown  in  Figure  6  where  the 
envelopes  of  the  oscillations  are  plotted  for  different  values  of  a  . 

APPROXIMATE  SOLUTIONS  FOR  LONG  DISPERSION  LENGTHS 
For  very  long  dispersion  lengths  (i.  e.,  large  values  of  q)  the 
numerical  integration  method  described  in  the  previous  section  loses  its 
usefulness.  This  is  due  to  the  fact  that  under  these  conditions  it  is  necessary 
to  take  the  difference  of  two  very  large  numbers  in  order  to  o'otai  .  a  small 
number  leading  to  possibly  large  errors.  However,  it  is  just  in  this  region 
of  very  large  n  that  approximate  solutions  to  the  transient  propagation 
problem  work  very  well.  We  shall  consider  approximate  solutions  of  the 
integral  in  (8)  for  the  propagation  of  a  unit  step-carrier  sine  wave. 

The  arrival  of  the  signal  wavefront  at  the  time  t  8  s/c  (t  *  ij) 
has  been  analysed  by  Sommerfdd.  U*  15  If  the  integrand  in  (8)  is  expanded 
for  high  frequencies,then  one  finds  that  under  the  condition 


T  .  «<<  £JL  , 


10.  15 


the  solution  to  (8)  may  be  written  in  the  form 


£  P  .  n)  ®  —  /T  Jj  (/^  I  U(u) 

p  n  l  / 


where 


u  a  2P  t>(t  -  tl)  » 


This  solution  is  plotted  in  Figure  7.  Note  from  the  condition  {20)  that  thie 
solution  is  valid  for  times  immediately  following  the  arrival  of  the  wavefront. 

After  this  initial  time  period  it  is  possible  to  use  the  saddlepoint 
method  of  integration^'  ^  to  evaluate  the  integral  iu  (8).  This  method  is 
based  on  the  fact  that  if  the  contour  of  integration  can  be  deformed  through 
the  saddlepoints  along  the  lines  of  steepest  descent,  then  the  major  contribution 
to  the  integral  will  result  from  integrating  in  the  vicinity  of  the  saddlepoints. 

If  the  second  derivative  with  respect  to  s  of  the  exponent  in  (8)  is  sufficiently 
large  then  the  saddlepoint  solution  will  be  a  good  approximation  to  the  actual 
solution. 

In  the  integral  in  (8)  the  location  of  the  saddlepoints,  found  by 
equating  to  zero  the  first  derivative  with  respect  to  s  of  the  exponent,  are 
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given  by 


{23} 


As  time  increases  after  the  arrival  of  the  wavefront  the  two  saddle  point! 
move  in  from  +oo  along  the  imaginary  axis  approaching  +  i  H  for  long  times. 
These  saddlepoirts  and  the  paths  of  steepest  descent  through  them  are 
shown  in  Figure  8. 

The  time  interval  from  just  after  t  =  r/c,  until  the  saddlepointc  approach 
the  vicinity  of  the  poles. is  called  the  precursor  region,  -o  this  region  the 
saddlepoint  solution  is  given  by*** 


£  (r,  *l)  =  j  £  {  cos 


(24) 


where 


,wl  i  w2-i  ! 


(25) 


and 


14  * 


'fy7 

p^2u-fV/4 


(T2  -  tl2)  (1  -  P2) 

.V 


(26) 


A  universal  curve  representing  the  envelope  of  the  precursor  is  plotted 
from  (25)  in  Figure  9.  Note  that  the  amplitude  of  the  precursor  increases 
with  increasing  t  for  w  <  1  and  decreases  for  increasing  r  for  w  >  1. 

The  amplitude  becomes  infinite  atw  =  1, which  occurs  when  ~  ~  9  ^1  -  P 
and  corresponds  to  a  time  of  arrival  at  the  group  velocity  (i.  e. ,  t  3  */w  }  * 

o  o 

By  setting  a  =  its  in  (23)  we  see  that  this  also  corresponds  to  the  time  at 
which  the  saddlepoint  crosses  the  pole.  The  standard  saddlepoint  solution 
given  by  (24)  breaks  down  in  this  region. 

This  time  interval  during  which  the  saddlepoint  crosses  a  pole  includes 
the  arrival  of  the  main  signal  at  the  group  velocity.  An  approximate 
solution  which  is  valid  in  this  time  interval  fo:  large  t]  may  be  obtained  by 


expanding  the  integrand  in  (8)  about  the  pole  at  s  =  +i«o.  The  corresponding 

,  .  .  .  .  1°.  <6 

solution  is  given  by 


£  it,  V)  3  A  sin  (t 


+  0  ) 

o' 


(27) 


Somme r£eld  solution 
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6  =  tan 
o 


-1 


C(v)  -  S(v) 


1  +  C(v)+S(v) 


and 


C(v) 


[  °°*  X1 


dx 


S(v) 


sin  —  x^j  dx 

A  *  * 


are  the  Fresnel  cosine  and  sine  integrals  respectively. 

The  amplitude  function  A  given  by  (28)  is  plotted  as  a  function  of  v 
in  Figure  10  and  will  be  referred  to  as  a  master  Fresnel  envelope.  The 
time  history  of  the  envelope  (as  a  function  of  t)  for  a  particular  value  of  P 
and  tj  is  determined  by  using  (29)  which  acts  like  a  "stretching  factor" 
and  simply  changes  the  abscissa  scale  in  Figure  10.  Note  carefully  that 
the  master  Fresnel  envelope  has  a  value  of  0.5  at  y=  0,  which  corresponds 
to 

In  order  to  determine  the  time  scale  at  which  the  ringing  of  the  master 
Fresnel  envelope  occurs  we  note  that  on  the  v  scale  the  ringing  occurs  for 
changes  of  v  of  the  order  Av^  1.  Thus  from  (29)  we  see  that  the  ringing 
will  occur  on  a  time  scale  A  t  ~  1/ iw  ,  where 


'"tF7  and  is  the  time  of  arrival  at  the  group  velocity  a  zjv 


mm 


The  solution  (27)  is  based  on  an  expansion  of  the  exponent  in 
(8}  in  a  Taylor  series  about  s  *  hiw^.  Only  the  first  three  terms  in  this 
expansion  are  kept  and  this  is  justified  as  long  as 


This  condition  is  based  on  the  requirement  that  a  large  number  of  cycles 

occur  between  the  arrival  of  the  wavefront  at  tf  *  zjc  and  the  time  of  arrival 

at  the  group  velocity  t  *  zfv  .  We  note  from  (30)  that  for  times  greater 
8  8 

than  tg  the  solution  (27)  will  be.  a  good  approxim  itioa  for  very  large  values 

of  tj,  depending  on  the  value  of  P.  For  values  of  P  greater  than  .866  the 

coefficient  of  »1  in  (30)  decreases  from  .25  toward  r.ero  as  P  approaches  unity. 

The  condition  (30)  therefore  becomes  more  and  more  difficult  to  satisfy. 

This  is  due  in  part  to  the  fact  that  the  pole  comes  very  close  to  the  branch 

point  and  the  expansion  fails  in  this  region.  This  problem  can  be  partially 

9  10 

overcome  by  a  suitable  conformal  transformation  ’  which  will  improve 
the  range  c*’  validity  of  the  solution  for  values  of  P  very  close  to  one  (P  >  .99). 

The  results  of  this  section  arc  summarized  in  Figure  11  where  the 
total  envelope  response  to  a  unit  step-carrier  sine  wave  input  is  plotted  for 
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the  particular  case  of  P=  0.  8  and  n  -  100.  This  total  solution  is  made  up 
of  three  separate  parts.  The  first  is  the  Sommerfeld  solution  for  the 
arrival  of  the  signal  wavefront.  This  is  then  joined  to  the  precursor  region 
which  is  evaluated  by  the  standard  saddlepoint  method.  As  the  saddlepoint 
approaches  the  pole  in  the  complex  s -plane  this  solution  is  then  joined  to 
the  master  Fresnel  envelope  solution.  Thij  total  solution  then  approaches 
the  steady  state  value  of  1  for  large  times. 

COLLISIONAL  DAMPING  OF  TRANSIENT  SIGNALS  IN  PLASMAS 

In  order  to  determine  the  effect  of  collisions  on  the  envelope  response 
found  in  the  previous  section. we  consider  the  integral  in  (7).  For  this  case 
the  trajectories  of  the  saddlepoints  move  off  of  the  imaginary  axis  in  the 
complex  s -plane.  The  arrival  of  the  main  signal  occurs  when  the  saddlepoint 
moves  into  the  vicinity  of  the  pole.  Since  we  are  interested  primarily  in  this 
part  of  the  solution, we  shall  consider  a  modification  of  the  pole  expansion 
of  the  previous  section. 

If  a  step-carrier  sine  wave  of  unit  amplitude  is  turned  on  at  t  »  0 
then  the  solution  given  by  (7)  may  be  written 


where 


W(s,  t)  «  st  - 


An  approximate  solution  of  (31)  may  be  obtained  for  the  main  signal 
build-up  by  expanding  W{s,  t)  in  a  Taylor  series  about  The  solution 
may  be  written  in  the  form^ 


£  {t,  z)  »  -  Im 


|eW^0)  £i+erf  0{-i^| 


where 


=  re-i+ 


f 2W"  H“0) 


and  where  <}>  may  be  shown  to  be  a  positive  angle  between  0  and 
The  envelope  of  (33)  is  given  by 


|  e  (t.  x)l  =  j  eWH“°)  |  |  J-  jfl  +  erf  9  (-i*0)] 


Wf-iu  ) 

The  exponential  function,  e  °  ,  is  just  the  attenuation  which  a  steady 
state  wave  of  frequency  «0  suffers  due  to  collisions.  The  remaining  factor 


can  be  written  as 


r.  9}  =  — 1-|  £l  +  erf  9  (-i“0)j 


M 


saeecgarssBiB 


where  6  *  z/z  -  <*>  and  we  have  used  the  relations 


erf  8  *  -  F(r) 

V» 


where  z  •  ie  *  re*0  and 


*  u{r,  0)  +  iv(r,  0}  * 


{36) 
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The  function  F(z}  given  by  (36)  is  tabulated  for  complex  arguments. 

The  behavior  of  the  function  B(r,  0)  given  by  {35}  is  shown  in 
Figure  12.  The  trajectory  (r,  0}  as  a  function  of  time  depends  on  the 
collision  frequency  v .  When  v  *  0,  0  *  45°  and  the  envelope  solution  is 
that  given' by  (23)  with  r  equal  to  the  argument  (29).  If  v  ^  O.the  angle  6 
is  not  constant  as  time  increases,  but  rather  decreases  from  90°  approaching 
some  angle  between  45°  and  9C°  asymptotically  for  large  times. 

It  is  therefore  seen  that  the  effect  of  collisions  is  two-fold.  First, 
the  entire  envelope  is  attenuated  an  amount  equal  to  the  attenuation  of 
a  cw  wave  with  an  angular  frequency  Second, the  “ringing"  of  the 
envelope  is  inhibited  and  quickly  damps  out. 


REFLECTIONS  FROM  A  PLASM  A- VACUUM  INTERFACE 
An  electromagnetic  signal  which  in  reflected  from  or  transmitted 
through  a  plasma-vacuum  interface  will  be  ui3persed  as  a  result  of  the 
discontinuity  in  the  medium  at  the  interface.  Let  a  lossless,  isotropic 
platana  occupy  the  half-space  a  >  Orand  free -apace  occupy  the  half-space 
a  <  0.  A  wave  incident  on  the  boundary  from  th}  free  space  side  will 
give  rise  to  a  reflected  wave  in  free  space  and  a  transmitted  wave  in  the 
plasma.  Applying  the  boundary  conditions  by  equating  the  tangential 
components  of  the  electric  and  magnetic  fir'  ds  on  either  side  of  the 
boundary, one  can  readily  obtain  the  complex  reflection  coefficient 


R  * 


./  2  .  „2 

s  -  Vs  ▼ n 

TTftTF 
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(37) 


If  the  incident  wave  is  a  step-carrier  sine  wave,  then  the  reflected 
signal  at  'he  boundary  is  given  by 
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ds 


(38) 
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Using  the  convolution  theorem  together  with  the  expressior ;  from  tables 


we  obtain 


£rw 


-21m 


JeiUo{t-  t'j 


J,(nt'} 


dt 


(39) 


Equation  (39)  has  been  evaluated  numerically  and  the  results  are  shown 

i-  Figure  13.  Reflection  from  a  perfect  conductor  is  indicated  by  the  dashed 

curve  in  Figure  13a.  The  major  feature  of  these  curves  is  the  apparent  shift 

or  "  delay"  of  the  initial  cycle  of  the  response  which  increases  as  the  plasma 

18 

frequency  decreases.  This  "delay  time”  has  been  discussed  by  Chen  and 

19.  20  z. 

Knop  f  who  find  it  to  be  related  to  the  time  t^  a  2/tL 

In  order  to  obtain  a  physical  meaning  for  this  time  delay  t^j  let 

us  write  the  reflection  coefficient  R  given  by  (37)  in  the  form 


If  (40)  is  substituted  into  (38}fthen  the  saddlepoints  of  the  resulting  exponential 
function  are  given  by 


V±  spj-  -  n2  (41) 

Thus  at  t  *  0  the  saddlepoints  are  at  +ra  on  the  real  axis.  As  time  increases 
they  move  inward  along  the  real  axis.  At  t  3  2.A  they  merge  at  the  origin  to 


Fig..  I '  Reflected  signal  at 

plnrirta-vacuam  boundary 

i 

form  a  third  order  saddlepoint.  For  t  >  2/  II  they  move  along  the 
imaginary  axis  and  approach  the  points  ±  ill  for  large  times.  For  times 
t  <2/11  the  saddlepoints  are  located  on  the  real  axis  and  the  solution  is 
nonoscillatory  in  nature,  increasing  monbtonically  in  the  negative 
direction.  The  oscillatory  nature  of  the  soluHon  occurs  for  t  >2/0.. 

At  t=  2/0  a  transition  occu-s  between  these  two  types  of  solutions.  The 
saddlepoints  moving  in  along  the  real  axis,  therefore,  represent  the 
"delay  which  occurs  before  the  oscillatory  nature  of  the  solution  can  be¬ 
gin.  Unfortunately,  the  saddlepcint  solutions  are  not  good  approximations 
for  this  problem  since  the  "heights"  of  the  saddlepoints  are  not  sufficient 
for  the  integral  t;<  be  evaluated  except  in  the  neighbor  hood  of  the  saddle¬ 
points.  However,  the  motion  of  the  saddlepoints  lends  insight  into  the 
nature  of  the  "delay  time"  as  described  above. 

CONCLUDING  REMARKS 

Two  types  of  transient  solutions  occur  depending  upon  the  distance 
the  wave  travels  through  the  plasma.  Far  short  dispersion  lengths  {less 
than  a  few  free  space  wavelengths),  exact  solutions  must  be  used.  This 
usually  means  some  type  of  numerical  integration  procedure  is  required. 
For  long  dispersion  lengths  (a  very  large  number  of  free  space  wavelengths) 
approximate  solutions  may  be  used  with  good  accuracy.  The  wavefront  of 
the  signal  arrives  with  the  velocity  of  light  c.  The  signal  then  builds  up 
slowly  in  the  precursor  region  until  the  main  signal  arrives  with  the  group 
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velocity.  As  the  propagation  distance  in  the  plasma  increases.  tin 
precursor  region  becomes  longer  and  the  amplitude  of  thn  signal  in  this 
region  becomes  smaller. 

The  nature  of  the  dispersed  signal  can  be  divided  into  two  categories. 
First  the  envelope  of  the  signal  becomes  distorted,  oscillating  about  the 
steady  state  value.  Second,  the  carrier  frequency  becomes  phase  modulated. 
Where  the  saddlepoint  method  of  integration  can  be  used, the  "instantaneous 
frequency"  of  the  dispersed  signal  is  equal  to  the  frequency  corresponding 
to  the  location  of  the  saddlepoint  in  the  s -plane  at  that  instant  cf  time. 
Whereas  the  envelope  of  the  distorted  signal  is  a  critical  function  of  the 
envelope  of  the  input  signal,  the  phase  modulation  characteristic  of  the 
dispersed  signal  is  primarily  a  function  only  of  the  nature  o.  the  dispersive 
plasma  medium. 


i 

i 


REFERENCES 


1.  H.  J.  Schmitt,  Dispersion  of  Pulsed  Electromagnetic  Waves  in  a 
Plasma,  IEEE  Trans,  on  Microwave  Theory  and  Techniques. 

Vol.  MTT-13,  472,  (July  1965). 

2.  G.  I.  Cohn,  Electromagnetic  Transients  in  Wt.veguides,  Proc.  Natl. 
Electronics  Conf.  Vol.  8,  284,  (1952). 

3.  R.  Elliott,  Pulse  Waveform  Degradation  due  to  Dispersion  in  Waveguide, 
IRE  Trans,  cn  MTT,  Vol.  MTT-5,  254  (Oct  1957). 

4.  M.  P.  Forrer,  Analysis  of  Millimicrosecond  RF  Pulse  Transmission, 
Proc.  IRE,  46,  1830  (Nov  1958). 

5.  C.  M.  Knop  and  G.  I.  Cchn,  Comments  on  Pulse  Waveform  Degradation 
due  to  Dispersion  in  Waveguide,  IEE  Trans,  on  MTT,  Vol.  MTT-11, 

p.  587  (1963). 

6.  C.  M.  Knop,  Pulsed  Electromagnetic  Wave  Propagation  in  Dispersive 
Media,  IEE  Trans,  on  AP  Vol.  AP-12,  No.  4  (July  1964). 

7.  A.  Rubinowicz,  Propagation  of  a  Cut-off  Train  of  De  Broglie  Waves, 

Acta  Physica  Polonica,  Vol.  X,  79,  (1950). 

8.  R.  C.ajewski,  On  Transients  in  Waveguides,  Bull,  ds  L’ Academic 
Polonaise  Des  Sciences,  Vol  3,  29  (1955). 


9.  M.  Cerrillo,  Transieut  Phenomena  in  Waveguides,  TR  No.  33; 
Research  Laboratory  of  Electronics,  MIT  (Jan.  8,  1948). 


230 


10.  R.  E.  Haskell  and  C.  T.  Case,  Transient  Signal  Propagation  in 
Lossless,  Isotropic  Plasmas,  AFCRL  Tech  Rept. ,  to  be  published. 

11.  E.  N.  Dekanosidze,  Tables  of  Lomir...  s  Functions  of  Two  Variables, 
Pergamon  Press,  New  York  (I960). 

12.  C.  T.  Case,  On  Transient  Wave  Propagation  in  a  Plasma,  Proc.  IEEE 
(Corr. }  53,  730  {July  1965). 

13.  -  Erdelyi,  Magnus,  Oberhettix^er,  and  Tricomi,  Tables  of  Integral 

Transforms,  Vol.  1,  McGraw-Hill  Book  Co. ,  Inc.,  New  York,  1954. 

14.  A,  Sommerfeld,  Ann.  Physik,  44,  177-202  (1914). 

15.  L.  Brillouin,  Wave  Propagation  and  Group  Velocity,  Academic  Press, 
New  York  and  London  (I960). 

16.  M.  V.  CerriUo,  An  Elementary  Introduction  to  the  Theory  of  the  Saddle  - 
point  Method  of  Integration,  Research  Laboratory  of  Electronics, 

MIT,  TR-55:2a,  {May  3,  1950). 

r*  x2 

17.  K.  A.  Karpov,  Tables  of  the  Functions  FJs}*  \  e  dx  in  the  Complex 

Domain,  The  Macmillan  Co.,  New  York,  1964, 

18.  S.  N.  C.  Chen,  The  Transient  Phenomenon  in  an  Isotropic  Plasma  without 
Collision  Loss,  Proc.  lE'SE  (Corr. )  51,  1045  (Jul.  1963). 

19.  C.  M.  Knop;  Comments  on  the  Transient  Phenomenon  in  an  Isotropic 
Plasma  without  Collision  Loss,  Proc.  IEEE  (Co.?. )  52,  99,  Jan.  1964. 

20.  C.  M.  Kaop,  Further  Comments  on  the  Transient  Phenomena  in  an 
lyotropic  Plasma  without  Collision  Loss,  Prwc.  IEEE  (Corr.)  53,  751, 
(July  1965). 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  RID 

(Security  classification  of  tide,  lody  of  abstract  and  indexing  annotation  must  be  entered  vdi^n  th'  overall  report  is  classified / 


1.  ORIGINATING  ACTIVITY  (Corporate  author)  2a.  REPORT  SECURITY  CUA»lP!CATlON 

Hq  AFCRL,  OAR  (CRD)  Unclassified 

United  States  Air  Force  rer  croup - 

Bedford,  Massachusetts  01730 


r.  report  title  proceedings  of  the  Third  Symposium  on  tiie  Plasma  Sheath— Plasma 
Electromagnetics  of  Hypersonic  Flight:  Volume  I.  Radiation  Characteristics  of 
Plasma -Covered  Antennas 


*.  DESCRIPTIVE  NO  TCI  (Type  of  report  and  inclusive  dates ) 

Proceedings  of  Symposium 


5.  AUTHORS  (Last  name,  first  name,  'initial) 

Rotman,  W.,  Moore,  H.,  Papa,  R.,  Lennon,  J. 


«.  REPORT  OATE 


i a .  CONTRACT  OR  GRANT  HO. 


l  project  ano  task  ko.  4642,  all  tasks 
C.  OOO  ELEMENT  62405394 
i  OOD  SUBELCMCMT  681000 


7a  TOTAL.  NO.  Or  PAGES 

237 


ta  ORIGINATOR'S  REPORT  HUMBER'S; 

AFCRL-67-0289  (Volume  I) 
Special  Reports  No.  64  (I) 


tb.  or Sl  (Aaj  ctirr  amber,  that  may  be 

AFCRL-67  -0280  (Volume  I) 


IX  ABSTRACT 


~~S>This  volume  is  a  collection  of  the  unclassified  papers  presented  at  the  Third 
Symposium  on  the  Plasma  Sheath— Plasma  Electromagnetics  of  Hypersonic  Flight. 
This  symposium  consisted  of  the  review  of  progress  in  reentry  communication 
studies  during  the  three  year  period  since  the  prior  meeting.  The  program  of  this 
symposium  on  plasma  electromagnetics  of  hypersonic  flight  involves  a  wide  range 
of  scientific  disciplines,  including  electromagnetics,  aerodynamics,  aerothermo- 
chemistry,  plasma  dynamics,  electronics,  and  high -temperature  phenomena.  The 
papers  were  selected  to  explore  as  many  of  these  facets  of  research,  including  the 
results  of  1r/horatory,  flight,  and  system  tests,  as  time  permitted.  •  } 

\  ■  f - 


Security  Cisasifirsuon 


Plasma 

Electro  ma  gnetlc  s 
Hypersonic  -Fliglii 
Reentry 

Communications 


instructions 


t.  OtlCMATlNC  ACTIVITY:  Eater  Ike  urn;  sad  aWress 
• I  tie  contractor,  nkoaincM.  gras  tec.  Dcnutaeit  of 
Defeaae  activity  or  other  orgsai ratio*  (corporate  amtkor) 
fwaiag  tie  raped. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Eater  tie  over¬ 
all  narky  claaaiffcatiort  el  tie  report.  Indicate  whether 
** Restricted  Data”  U  iaelodod.  MarLlag  la  to  ie  ia  accord- 
*a ec  wlti  appropriate  aoccity  regri^loa*. 

2b.  CROUP:  Aaloaale  dowagrsdiag  la  specified  ia  DoD 
Directive  3300.10  aad  Armed  Forcesladaatrisl  Haul. 

Ee*«r  tie  group  aatiirr.  .'.iso,  vdiea  applicable,  shew  that 
tjunal  enrkiaga  have  been  aaed  for  Grosp  3  aad  Creap  4 
a*  astir*  i  zed. 

1  REPORT  TITLE:  Eater  tie  complete  report  title  la  all 
capital  letters.  Titles  ia  all  caaca  aieold  be  aaclaullied. 

U  e  acaaiaelal  title  cacao:  be  selected  withost  class  iilcs- 
tlce,  skew  title  elasaificatioa  ia  all  capitals  ia  paresthesia 
lame  likely  (oliowiag  tie  title. 


S.  AlTTHOSiS):  Eater  tie  aeaieU)  of  avtborla)  as  shows  A 
or  ia  lie  report.  Eater  last  same,  first  «£rae,  middle  iaitial. 
II  military,  show  nub  aad  bruch  s <  service.  The  ease  of 
the  priacipsi  astber  ia  as  obsolete  taiakasm  reqeirement. 

4.  REPORT  DATE:  Eater  tie  date  of  the  report  u  day, 
■oath,  year,  or  aoaeh,  year.  If  Here  than  oae  date  appears 
ea  tie  report,  seo  date  si  pabUckioa. 

7a.  TOTAL  HUH  BEK  OF  PACE?-.  Tie  lota!  page  cooat 
shea  Id  follow  anal  paglaatioe  procedarea,  i.e.,  eater  the 
a  amber  of  pages  coataiatag  iaforaetioa. 

74.  NUMBER  OT  REFERENCES:  Eater  the  total  somber  of 

refers  aces  cited  ia  tie  report. 

8s.  CONTRACT  OR  C8ANT  MMBEB:  If  appropriate,  eater 
the  applicable  aameec  el  the  cosrract  or  great  seder  which 
lie  report  was  writ  tea. 

84.8c.  4  84.  PSOJECT  HUMBER:  Eater  the  appropriate 
Military  department  idemliiicatiot,  asri  ax  project  s saber, 
smbpreyect  a  saber,  system  « sabers,  tasi  aaaber,  etc, 

9a.  ORIGINATOR’S  REPORT  NUMBEKS):  Eater  tbe  off N 
ciai  report  aaahrr  by  which  tie  decs  neat  will  be  idestifi e„ 
aad  ceatraliedby  ibe  origiaettag  activity.  This  avaber  neat 
be  aaiqse  to  tii*  report. 

94.  CITHER  REPORT  NUMBER®:  II  tbe  report  has  beets 
aaaigacd  say  other  report  avabers  (eitktr  4r  the  originator 
or  1 7  ike  tpoatorX  also  eater  this  aaaberis). 


10.  AVAILABILTrY/LMTATION  NOTICES:  Enter  any  iimi- 
tatloaa  on  further  die  semination  of  the  report,  other  than  those 
imposed  by  secviiy  classification,  using  standard  statements 
Bach  as: 

(I)  *'Qaslifred  requesters  au.y  obtaia  copies  of  this 
report  from  DDC" 

(S  "FereUa  aaeoaarerienl  sad  diasesd nation  of  this 
report  by  DOC  is  not  authorized.'’ 

<3)  "U.  5.  Government  agencies  may  obtain  copies  of 
this  report  directly  Iron  DOC.  Other  qualified  DDC 
nsers  shall  request  through 

f* 

(4)  **U»  5*  military  agencies  may  obtain  copies  of  tbia 
report  directly  from  DDC*  Other  qualified  users 
stall  request  through 

»♦ 

(5)  **AII  distribution  of  this  report  is  controlled*  Quali¬ 
fied  DDC  users  shsll  request  through 

M 

II  the  report  hss  been  furnished  to  the  Office  of  Technical 
Services,  Oepstfment  of  Commerce,  for  ssle  to  the  public,  indi- 
cste  this  fact  sad  eater  the  price,  if  known, 

11.  SUPPLEMENTARY  NOTES:  U»e  for  addltioftsl  explana¬ 
tory  notes* 

11  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  fpo>- 
f*g  for)  the  research  and  development*  Include  address* 

11  ABSTRACT:  Eater  sa  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  tbe  report,  even 
though  it  may  also  appear  elsewhere  in  the  body  of  the  tech¬ 
nical  report*  U  additional  space  i?  required,  a  continuation 
sheet  *5*11  be  attached* 

It  is  highly -desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragraph  of  the  abstract  shall 
end  with  aa  indication  of  tbe  military  security  classification 
of  the  information  in  tbe  paragraph,  represented  as  (TS)t  (SX 
(CX  or  (U) L 

There  b  no  limitation  on  the  dearth  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  ISO  to-  225  words. 

Id.  KEY  WORDS:  Kev  wotda  ere  technically  meaningful  terms 
or  short  phrases  that  cherseterixe  s  report  and  may  he  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  bo  security  classification  is  required*  Identi¬ 
fiers,  such  as  equipment  model  design  at  ton,  trade  name,  mill- 


context*  The  assignment  of  links,  rules,  and  weights  is 
optional* 


Unclassified 


unty  utsimeatton 


